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FOREWORD 

Affixed vibration study is an important or even a major 
influence in many branches of modem engineering practice where 
fluctuating stresses originated by periodically varying com¬ 
ponents of the motive forces may seriously jeopardise structural 
stability. 

The general term “ vibration ” is given to the effect produced 
by periodic forces, and this may range from physiologically 
unpleasant tremors to disturbances of sufficient magnitude 
to cause mechanical failure of structural components. 

Recognition of the ill effects that are apt to appear when 
vibration problems inherent in engineering structures are neg¬ 
lected has led to intensive studies of fatigue phenomena in 
materials, and other relevant subjects. This has resulted in 
more widespread appreciation of the importance of avoiding 
sharp corners and other abrupt structural discontinuities from 
the point of view of minimising local concentrations of stress; 
and the necessity for adjusting the dynamical characteristics of 
the system so that important zones of resonant vibration do not 
occur in the operating speed ranges. 

The foregoing remarks are particularly true of sea, land, and 
air transport services, where modem requirements demand a 
clear understanding of attendant vibration problems if opera¬ 
tional troubles of one kind or another are to be avoided. In 
this connection it should be borne well in mind that modem 
requirements are not necessarily fulfilled by designs which 
merely provide for freedom from mechanical breakdown. 

Most engineers now and then experience the satisfaction of 
handling a machine in which the different design factors have 
been so carefully balanced that the resultant blend is a product 
of outstanding merit. Even to-day, however, such examples are 
comparatively rare and in many cases are confined to individual 
specimens from a production batch, so that the achievement 
appear^ to be largely the result of a happy accident which some¬ 
how or other provides a product which is neither too sluggish 
and lifeless nor too lively and over-sensitive. This suggests that 
an important contributory factor towards the achievement of • 
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these ideal results is accidental realisation of the correct solution 
of the vibration problems associated with each case. Vibration 
study is therefore likely to exert considerable influence on 
engineering development, and already several industrial con¬ 
cerns making products which are particularly susceptible to 
vibratory influences have established special vibration engineering 
departments. 

In the past the tendency has been to disregard vibration 
during initial design and consequently the work of the vibration 
specialist has been largely confined to the correction of faults 
which have all too often appeared in the product after it has 
been put into service. Only occasionally has this procedure 
proved to be entirely satisfactory, because it is almost impossible 
to cure a really difficult case without resorting to a drastic 
re-design of the whole system. There is nearly always some 
increase of weight accompanied by a reduction of operational 
efficiency due either to the necessity of imposing restrictions 
on operating speeds or introducing energy-absorbing damping 
devices. In other words, failure to take the vibration problem 
into account during initial design nearly always results in 
inefficiency of one kind or another. 

It is to be hoped, therefore, that in the future many more 
engineers will become interested in vibration study so that 
design work will proceed along lines which are fundamentally 
sound. At the same time, there are certain aspects of vibration 
study which must always remain the special province of the 
vibration engineer ; in particular, the task of accurately recording 
and analysing structural vibration is a matter which requires 
the skill and experience of trained personnel. 

Although a considerable volume of literature has accumulated 
on this important subject, there is a lack of information on the 
fundamental principles underlying present-day vibration study 
presented in such a way as to be readily assimilated by the busy 
engineer. The author of the present work is therefore to be 
congratulated on having produced a book which should be 
welcomed by everyone wishing to acquire a good working know¬ 
ledge of these principles. In particular, the author is to be 
commended on dealing in a straightforward manner with 
mechanical vibrations as such without resort to electrical or 
other analogies. 

There is no doubt that the student or engineer who studies 
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Mr. Manley’s work will find little difficulty in extending the 
principles so ably set forth to meet the needs of his special 
problems. This is largely because the work is based on the 
results of actual first-hand experience in the Vibration Depart¬ 
ment of a leading industrial establishment. 

W. KER WILSON. 

London, 1942. 



AUTHOR’S PREFACE TO SECOND EDITION 


This second edition represents a compromise which, it is hoped, 
partially meets conflicting requirements. Several reviewer, and 
others, have suggested that in some places an amplification of 
the descriptive matter would be advantageous. On the other 
hand, the reception of the original, and the fact that another 
edition was called for inside two years (although it could not be 
provided at that time), have prompted an avoidance of any 
alterations which would add unduly to the delays incidental to 
publication. 

In the text, additional material has been confined almost 
entirely to the insertion of “ Introductory ” and “ Summary ” 
sections in each chapter ; it is hoped that these will aid the reader 
in gaining an appreciation of the background to the subject and 
of its mode of development. Two further appendices have been 
added ; in the first of these are given some details of the applica¬ 
tion of the theory to the solution of the simpler practical problems, 
while the second includes tables of some of the quantities com¬ 
monly used in vibration calculations. Slight extensions have 
been made in the analytical exercises at the ends of the earlier 
chapters, and sundry small errors have been corrected. 

The Bibliography has been brought more up to date. Since 
revision, J. Morris’s book, The Escalator Method in Engineering 
Vibration Problems, has appeared (Chapman and Hall, 1947), and 
provides an admirable guide to the calculation of the vibration 
characteristics of very complicated systems. 

The author’s thanks are due to many colleagues whose sug¬ 
gestions have been adopted, or are kept in mind in those instances 
where it has not proved possible to give affect to them in this 
present revision. 

Newcastle-upon-Tyne, 

September, 1947. 

R. G. MANLEY. 
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AUTHOR’S PREFACE 

The need for a really introductory book on vibration theory has 
for some time been felt, both by technical staff engaged in vibra¬ 
tion research and by others whose work brings them into contact 
with the problems of the subject. The available textbooks are, 
of course, excellent and comprehensive, but experience has shown 
that it is difficult for beginners to extract from the somewhat 
lengthy standard texts the precise information they are seeking. 

The present work is designed to serve as an introduction to 
the subject, so that the existing textbooks can subsequently be 
approached with some degree of confidence. In accordance with 
this plan, attention has been confined to the basic theory, a 
knowledge of which is essential to a thorough understanding of 
the physical phenomena ; instrumentation and testing technique 
fall outside the scope of the book. No previous knowledge of 
differential equations is assumed, this troublesome part of the 
theory being dealt with by the powerful and easy method of 
operators; indeed, the only mathematical equipment is that 
which should be at the command of anyone who has gone but 
a little way 'beyond the Matriculation syllabus. A series of 
notes in the appendices serves as a bridge between the school 
work and those parts of analysis which are required in the text. 

The most important advance, in recent years, in vibration 
theory has been the development of the concept of effective 
inertia, with its collateral “ mechanical impedance ” (herein 
termed dynamic stiffness for reasons stated in the text); these 
methods of attack have been utilised with conspicuous success 
in the treatment of torsional vibration problems, particularly in 
aircraft applications. An entirely new approach to the idea of 
effective inertia is described in the text—an approach which, 
besides being direct, has been “ tried out ” with students and 
is evidently easily grasped. 

Particular emphasis is given to resonance, which is the really 
significant phenomenon for the engineer; the usual assumption 
concerning the practical equivalence of natural and resonant 
frequencies in lightly-damped systems is justified by an argu¬ 
ment based on the method of partial inertias which is described 
in the chapter on effective inertia. 

ix 
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A few exercises are included at the end of each chapter; 
these are for the most part analytical in nature, and serve slightly 
to extend the work covered in the chapter. Some numerical 
examples are collected at the end of the book. 

In the Bibliography are listed four standard works, only. 
These are books which should be readily available to all who 
have any working connection with the subject, and between 
them they cover the entire field of vibration theory and practice. 

The volume embodies the subject-matter of a course of 
lectures delivered by the author to members of the Vibration 
Department of the de Havilland Aircraft Company ; in this and 
other connections the author’s thanks are due to many members 
of this Company, particularly: to Dr. W. Ker Wilson for his 
interest and encouragement, and for the many helpful suggestions 
he has made ; to Mr. R. N. Hadwin for the facilities so readily 
granted for the compilation of the book, especially access to his 
library of technical works; and to Miss M. K. B. Harwood, 
both for much numerical work embodied in the text and exer¬ 
cises, and also for ungrudging co-operation in the tedious task 
of checking the typescript and proofs. 

R. G. MANLEY. 

Edowabe, Janmry 1942. 
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The symbols listed below are those that occur most frequently in the 
text. All symbols are defined where they first occur, and in many cases 
the meaning is reiterated in subsequent references. 

The figures in parentheses denote the pages on which the symbols are 
first used or defined. 
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CHAPTER I 

SYSTEMS HAYING ONE DEGREE OF FREEDOM 

(Undamped Motion) 


Introductory. 

This first chapter is concerned with the vibration of a simple 
system when disturbed from its position of rest. The system 
consists essentially of a concentrated mass and an elastic agent 
which opposes motion of the mass by a force proportional to the 
displacement from the position of rest. 

The most commonly encountered example of such a system 
in engineering can be regarded as equivalent to a mass suspended 
on a light coiled spring, as depicted in Fig. 1, page 2. Many 
industrial machines and power-plant installations are mounted 
on flexible suspensions of various types, some practical features 
of which are described in Appendix IV. [Some equivalent 
systems illustrated diagrammatically in Fig. 3, page 11, are also 
of practical importance. The cantilever spring in Fig. 3a merely 
replaces the coiled spring of Fig. 1, so that the systems are essenti¬ 
ally the same; a very useful model can easily be built with a 
light steel strip for the spring, and can be employed for investigat¬ 
ing practically the properties of simple vibrations. The pendulum 
at (6) in the diagram is the basis of many time control or indicator 
mechanisms; the ordinary pendulum clock depends upon the 
vibration characteristics of the pendulum for its accuracy as a 
time-keeper. The torsional pendulum at (d) replaces the simple 
pendulum in watches and similar mechanisms, and also corres¬ 
ponds dynamically to the suspension used in some galvanometers.] 

The important physical properties of such a system are that it 
vibrates, when deflected and then released, at a definite frequency 
(i.e. with a definite number of vibrations per minute) dependent 
solely upon the mass of the suspended body and the stiffness of 
the spring ; and that the motion is of a special type. The fre¬ 
quency, which does not depend upon the size of the vibration 
provided this is sufficiently small to ensure that the elastic law 
is not violated, is termed the “ natural frequency ” of the system ; 
the motion is in the form of a sine-wave (although this is modified 
in practice by the damping effects of friction forces). These 

l . b 
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properties are the result of the nature of the restoring force, which 
is-proportional to the deflection, and can be established mathe¬ 
matically on the assumption of such a force. For a thorough 
understanding of the subject it is essential to follow out the 
mathematical study, which does not however entail more than 
an elementary acquaintance with the calculus. Fortunately, the 
method of operators, introduced in Section 2, enables the differ¬ 
ential equations arising in the study of vibration to be solved 
without difficulty. 

1. A simple vibrating system. 

The system shown diagrammatically in Fig. 1 is an example of 
the simplest type of mechanical vibrating system that can exist. 
It consists simply of a mass m fixed to the lower end of a light 
coiled spring, the upper end of which is attached to a rigid 
support. If the mass is pulled downwards and then released, 
practical experience shows that it will oscillate in the vertical 
direction, the amplitude of oscillation diminishing in time until 
finally the system appears to come to rest. This decay of the 
motion is due to dissipation of energy by the action of damping 
forces, the nature and effect of which are described in Chapter II. 
It will be convenient to consider first the motion of the system 
in the absence of such forces. 



Fig. 1.—Simple spring-mass system: (a) static position; (6) at time L 
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The spring is therefore considered to be idealised, so that the 
force required to extend or compress it along its axis is propor¬ 
tional to the extension or compression only. Fig. la shows the 
position of the system in static equilibrium, the spring being 
extended a distance d beyond its natural length by the action 
of the gravity force mg. The spring-constant of proportionality, 
or stiffness, k, is then given by the equation: 

kd = mg . . . (1.1) 

The mass is then pulled downwards a further distance x, and 
released with zero velocity. After a time t from the instant 
of release the displacement of the mass below its equilibrium 
position is x, as shown in Fig. 16. The real forces acting on 
the spring are : 

(i) the gravity force mg downwards, and 

(ii) the spring force k(d+x) upwards, as the reaction of the 

spring is proportional to its extension beyond the 
natural length and tends to return the spring to its 
natural length. 

These two forces have a resultant kx upwards, as the part kd of 
the spring force is balanced by the gravity force mg, as shown by 

(1.1) . The downward acceleration of the mass is given by the 
equation of motion: 

k 

mx = — kx, or x = — x . . (1.2) 

m 

(The dot notation for deferential coefficients with respect to time will be 
employed wherever possible.) 

(1.2) is a linear differential equation of the second order in one 
dependent variable, i.e. there are no second or higher powers of 
any variable quantity, the highest differential coefficient involved 
is the second (two dots) and there is only one dependent vari¬ 
able (x). On account of this last property the system is said to 
have one degree of freedom. In order to obtain an equation 
giving the displacement x directly in terms of time, it is necessary 
to integrate the equation of motion twice, and at each indefinite 
integration an arbitrary constant must be inserted for generality. 
The general solution must therefore contain two arbitrary con¬ 
stants, the values of which can be chosen to make the solution 
correspond to particular conditions. It has been postulated that 
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at the instant of release (t = 0) the displacement is # 0 and the 


velocity is zero. With the substitution 

* 


2 k 
m 

• 

• (1.3) 

the three equations to be satisfied are : 

equation of motion, x = — co 2 x 

(«n 


, initial conditions, x = x a at t = 0 


• (1-4) 

o 

II 

Is 

o 

II 

•a 

(c)J 


2. The differential operator D. 

It can easily be shown by substitution that 

the function 

x = x 0 cos cot . 

. . 

• (2-1) 


is a solution of (1.4). Differential equations are frequently solved 
by guesswork—indeed, guessing is very often the only way of 
arriving at a solution—but there is a perfectly general method 
of solving equations of the type of (1.4a). 

Write D for — so that Dx stands for x and D 2 x for x. The 
dt, 

symbol D represents the operation of differentiating with respect 
to time, and is called the differential operator . (1.4a) can be 

written, with this notation, as 

D 2 x — — co 2 x .... ( 2 . 2 ) 

Proceeding in an experimental manner, the common factor x is 
“ dropped ” by analogy with ordinary algebraic equations, leaving 

D 2 = — a) 2 . . . . (2.3) 

(2.3) is to be interpreted as stating that the operation of double 
differentiation with respect to time, when performed on the 
appropriate function, is equivalent to the operation of multiply¬ 
ing the same function by — co 2 . It remains to determine the 
appropriate function. Performance of the square-root operation 
on (2.3) leads to : 

D = ± ico, i.e. ^ = ± icox . . (2.4) 

dt 

where i 2 = — 1 


Separating the variables x and t, and taking first the positive sign, 


dx 


ico . dt 
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whence Iog e x = icot+ a constant 

= to)<+log e G x (say). 

This last result can be written as: 

x = 

Direct substitution shows that this function does in fact satisfy 
(1.4a). Taking the negative sign in (2.4) leads to another 
solution: 

x = C t e~ M 

and substitution again shows that the combined function 

x = C l e ia,t +C t e- i0,t . . . (2.5) 

satisfies the equation of motion (1.4a), and it contains the two 
arbitrary constants necessary for generality. Before determining 
the values of these constants so that the initial conditions (1.46, c) 
are satisfied, it is advisable to restate (2.5) in more suitable terms. 
Elementary mathematical analysis (see Appendix I, section 33) 
shows that 

e <9 = cos 0-M.sin 0 
e - * 9 = cos 0 — i .sin 0 

Putting (at for 0, (2.5) can therefore be written: 



x = A . sin cat +B . cos (at . 

. (2.6) 

where 

A = t(C\ - C t ) 



B = 0,+C, 


and 

2 k 

Ct) 2 = — 



m' 

The form x = A.sin (ot+B .cos cot is the general solution of the 
equation of motion (1.4a). From the initial conditions (1.46, c) 
the values of A and B are determined as 

A = 0, B = x 0 

and the particular solution for these initial conditions is therefore 


x = x „. cos (at 

where o>* = - ’ ' * ‘ (2,7) 

m 

The solution can be obtained by more direct methods in this 
particular case, but the method of operators described above is 
a very powerful one and can easily be extended to more com¬ 
plicated linear differential equations ; it is for this reason that 
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the method has been described in some detail Using exactly 
the same procedure it can be shown that the solution of the 
equation - 

x Dx = ax 

is x = Ce at 

and that the solution of the equation 

D 2 x+PDx-\-Qx = 0, 

where P and Q are independent of x and t, is 

x = C^+C^, 

where a, /? are the roots in D of the equation D 2 -\-PD+Q = 0. 
(See Exercise 1 at the end of this chapter.) 


3. General solution and interpretation. 

The function x — A.sin cot+B .cos cot has been shown to be 
the general solution of the equation of motion (1.4a), and it has 
already been noted that the constants A and B are to be deter¬ 
mined by specific conditions of the motion. The conditions most 
commonly given are the displacement and velocity at the instant 
t = 0. At this instant let 


* _ X -° 1 at t = 0 . . . (3.1) 

X — Xq j 

Then x„ = B 

and x 0 = Aco, 

and the solution is : 

x = — sin <w<+a; 0 .cos cot . . . (3.2) 


The general solution (2.6) can be put into a more convenient 
form by using the relation 

sin(0+^) = sin 0.cos ^-fcos 0.sin^. 


Dividing (2.6) by VA 2 +B 2 gives 
x A 


VA 2 +B 2 VA 2 +B 2 


sin cot- 


VA 2 +B 2 


cos cot, 


A 

VJJ+B 2 

B 

VIZ+B* 


= cos <f> 
— Bin<f> 


and if 
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the result is obtained : 

x = a.sin (a>£+<£)\ 
where a = \/A 2 -\-B 2 I 

tan <f> = j ) 


(3.3) 


Note. —In determining <f> from (3.3), the sign of A or of B enables the 
correct choice to be made from the two alternative values of </> correspond¬ 
ing to the value of tan <f>. Thus if B is positive, <f> lies within the limits 
0 and 180°, whereas if B is negative, <f> lies within the limits 180° and 360°. 

In this form of the solution there are again two constants 
(a and </>) to be determined. If the conditions at the instant t = 0 
are x = x 0 , x — x 0 , then 

x 0 = a. sin <f> 
x 0 = ao) . cos <f> 


whence 



(3.4) 


The angle <f> is termed the phase-angle and the factor a is termed 
the amplitude , a useful interpretation of both terms being readily 
available. In Fig. 2 the point P moves in the counter-clockwise 



Fig. 2.—Generation of sine-wave by rotating vector. 


direction on a circle, radius a and centre O, so that the line OP 
rotates about 0 with a constant angular velocity m. If QA is 
a reference line, the projection AQ of the rotating vector OP on 
a line perpendicular to OA is a. sin 0, where 0 is the angle POA. 
If P is at P 0 when t = 0, and the angle P 9 0A is then 
AQ = a.sin (o)t+<f>). Thus AQ represents x, and the right- 
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hand part of the diagram shows the graph of the displacement 
x plotted against time t. 

A fundamental property of the sine function is that if the 
argument is increased by 2ji radians the value of the function 
is unaltered, i.e. 

sin (0+2Xn) — sin 0, 

where X is any integer; hence after successive intervals of time 
— the motion repeats itself exactly. 


Displacement x = a.sin (cot+<f>) 
Velocity x = aco. cos (cot -| -<j>) 
Acceleration x = — aco 2 .sin (cot+<f>) 


and 



(sin) 

(cos) 


(cot-\-<f>). 


The time-interval — is termed the period of vibration, and 
co 

the reciprocal of the period is the frequency, or number of cycles 
(i.e. complete repetitions of the motion) in unit time. As the 
unit of time in every system of measurements is the second, the 
frequency is given as a number of cycles per second, commonly 
abbreviated to c.p.s. In practical work it is more convenient 
to state the frequency in cycles per minute, usually written c.p.m. 
The whole motion of the mass in Fig. 16, under the initial con¬ 
ditions specified, can be described concisely thus : the displace¬ 
ment of the mass below the static position is a sine-function of 
time, the maximum displacements being ±* 0 and the frequency 
of vibration being 



i II 

2n\J m 


C.P.S. 


60 Ik 

2ji\l m 


c.p.m. = 



c.p.m. (3.5) 


This frequency of free vibration of the system, with no damping 
or external applied forces, is termed the natural frequency of the 
system. A table giving corresponding values of F and co 2 = k/m 
is included in Appendix V, page 145. 


4. Units and general remarks. 


In making use of the frequency formula 


1 k 
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due consideration must be given to the units of measurement. 
In scientific units the spring constant k is in poundals per foot 
deflection, or in dynes per centimetre, and the mass m is in 
pounds or grams. In the engineers’ system of units & is in 
pounds per inch and m is in slugs, one slug being that mass to 
which unit force (one pound weight) imparts unit acceleration 
(one inch per second per second). Thus one slug equals g pounds, 
where g is the gravitational constant in inch units = 386 inches 
per second per second; so that one slug equals 386 lbs. (see 
Appendix II, section 37). Any of these three systems may be 
used so long as care is taken to ensure that all the quantities 
involved are given in consistent units. 

Checking the units of the frequency formula (3.5), 
k = lbs.ins. -1 


and thus 


lbs. 

m = — 

g 


lbs.ins. -1 secs. 2 


‘'-Vi 


lbs.ins. -1 
lbs.ins. - ^ecs. 2 


= secs. 


-l 


That this result is correct is evident from the fact that the 
frequency is the reciprocal of a time-period. 

As a numerical example, suppose that the mass m is 200 lbs., 
and that the static deflection of the spring under the weight of 
200 lbs. is 0-5 inches. Then 


k = 2 x 200 lbs./in. 
200 . 

m ~ 386 ShgS ' 


and 


6-28 


.V772 = 4-42 C.P.S. 
= 265 C.p.m. 


This example suggests an alternative frequency formula in terms 
of the gravity deflection. If <5 S is the static deflection of the 
spring under the gravity-load of the mass, then 


k- 1 ^ 

S’ 


and 


F = 


1 

2n 



This formula is sometimes quoted as F 


• (4.1) 


188 

vs 


where F is the 
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frequency'in cycles per minute, and <5„ is the static deflection in 
inches. A table giving corresponding values of natural frequency 
and static-deflection is included in Appendix V, page 146. 

The amplitude and frequency are usually of far greater 
practical importance than the phase-angle. The two forms 

A.sin (ot+B. cos <ot 
and a .sin (cot+<f>) 

both contain two constants which have to be determined by 
initial conditions, but in the first form both A and B must be 
found before the amplitude can be calculated; the form 

x — a.sin {wt+<f>) 

is therefore normally to be preferred. There are two conventions 
with regard to the sign of the amplitude and the range of the 
phase-angle': (i) if the amplitude is considered as essentially 
positive, the phase-angle may have any value 0 < <f> <[ 2ji ; 
(ii) if the amplitude is allowed to take a negative value, the 
phase-angle is restricted to the range 0 < ^ < ?r, for 

a.sin (arf+a+jr) = — a.sin (cut+a.). 

The choice between these alternatives depends upon the nature 
of the problem. 

It is well to note here what assumptions have been made in 
the preceding analysis. First, the effect of damping forces has 
been neglected ; these forces are taken into account in Chapter II. 
Secondly, the spring is assumed to have asymmetrical linear 
characteristic, i.e. the forces required to extend and to compress 
it by equal distances are equal and proportional to the distance. 
This assumption is justified for small vibrations in very many 
practical applications ; the analysis of unsymmetrical and non¬ 
linear spring characteristics is beyond the scope of this present 
work, but an excellent discussion is given by Den Hartog (see 
reference 1 in the Bibliography at the back of the book). The 
physical conditions postulated for the system analysed above do 
not arise in practice, but the type of resultant motion (i.e. sinu¬ 
soidal) does occur. It is shown in Chapter VI that any cyclical 
variation can be expressed as the sum of a number of sinusoidal 
variations of different frequencies, and for this reason the 
properties of sinusoidal motion are very important. 
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5. Simple harmonic motion, and equivalent systems. 

The type of motion in which the restoring force is proportional 
to the displacement from a mean position, and which therefore 
is expressed by an equation similar to (3.3), is termed simple 
harmonic motion. Prom the equation of motion (1.2) it is found 
that the acceleration for a unit displacement ip numerically equal 
to k/m ; the frequency equation (3.5) can therefore be written : 


or 


frequency _ / — (acceleration for unit" 

(C.P.S.) 2n\J displacement) 

frequency _ ^ ^ / — (acceleration for unit ’ 
(c.p.m.) V displacement) 


(5.1) 



Fig. 3.—Systems equivalent to Fig. 1: (a) light cantilever ; (b) simple 
pendulum ; (c) compound pendulum ; (d) torsional pendulum. 


and this property is typical of simple harmonic motion. If it 
is known that the restoring force on a body is proportional to 
the displacement from a fixed position, formula (5.1) can be 
applied at once to determine the frequency of vibration. The 
system of Pig. 3a, which consists of a mass m on the end of 
a light cantilever spring, is fundamentally the same as that of 
Pig. 1, for transverse vibrations; the spiing-constant of the' 
cantilever can be calculated from the physical constants (see 
Appendix II, section 38). Two other examples of systems to 
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which the formula (5.1) can be applied are (i) the pendulum, 
Figs. 36, c and (ii) the torsional pendulum, Fig. 3d. 

Rotatory systems. 

It is shown in Appendix II, section 37, that if a rigid body 
rotates about a fixed axis the torque T and angular acceleration 
B about that axis are linearly related by the equation 

T = IB, 

where I is the moment of inertia of the body about that axis; 
and if the body is considered as being composed of infinitesimal 
particles typified by the particle of mass 8m at a radius r from 
the axis, I = Ur 2 .8m, the summation being extended over the 
whole body. There is a general analogy between rectilinear and 
angular motion, the corresponding quantities being listed in 
Table I together with the British scientific and engineers’ 
units. 


Table I 


| Linear 

Angular 

Action 

Brit. Scientific 

Engineers’. 

force 

poundals 

lbs. 

torque 
poundals-ft. 
lbs.ins. 

Response 

Brit. Scientific 

Engineers’. 

acceleration 
ft. /sec. 2 
ins./sec. 2 

angular acceleration 

J radians/sec. 2 

Inertia 

Brit. Scientific 

Engineers’. 

mass 

lbs. 

lbs.ins.-^ec. 2 

moment of inertia 
lbs.ft. 2 
lbs.ins.sec. 2 


Fig. 36 illustrates a simple pendulum, consisting of a particle 
of mass m hung on a light inextensible string of length l. If 
the string is deflected in a vertical plane through the point of 
support A, and the deflection from the vertical position is 0 
radians, the restoring torque about A is that due to the weight 
mg resolved normal to the string, i.e. mgl. sin 6. As this is a 
restoring torque tending to decrease 0, the equation of motion 
analogous to (1.2) is 

mgl. sin 0 — — IB = — mltf 
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For small displacements sin 0 = 0, so that the equation of motion 
becomes 

mglO = — ml 2 6 

or (5.2) 

b 

As g and l are constant the motion is simple harmonic and the 
frequency is given by 

F — c.p.s., where co 2 — ~ . . (6.3) 

2 71 l 

Similarly, in the compound pendulum illustrated in Fig. 3c, if 
m is the mass of the pendulum, I a its moment of inertia about 
an axis through the support A and normal to the plane of motion, 
and h is the distance between the point of support and the centre 
of gravity of the pendulum, then 

mghd = — I A S, or 0 = — r H^h$ . . ( 5 . 4 ) 

for small values of 0 . 

Hence the frequency of vibration is given by 

F = c.p.s., where a> 2 = . . (5.5) 

2 n I A 

Torsional systems. 

Fig. 3d illustrates a torsional pendulum. A rigid flywheel A 
is attached to one end of a light shaft B, the other end of which 
is clamped. The gravity load of the weight of A is taken on 
a bearing which is supposed frictionless, and the flywheel A is 
free to rotate against the torsional stiffness of the shaft. Let 
this torsional stiffness be C, so that the torque that must be 
applied at A to twist it through an angle 6 is CO, and let the 
moment of inertia of A about the axis X — X be J. Then the 
equation of motion analogous to ( 1 . 2 ) is 

JO = - CO 

or 0 = -jO ... . (5.6) 

J 

Hence the frequency of vibration is given by 

F — C.P.S., where co 2 = ^ . . (5.7) 

2 71 J 

Note .—The quantity to, which is the angular velocity of the radius 
vector OP in Fig. 2, page 7, is called the “ phase velocity of the vibration. 
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The units of measurement are radians per second, so that the quantity 
has the dimension of the reciprocal of time; the connection between it 
and the frequency F is given by 

ft) CO I 

F = 2n = £28 C • P • S • . (5-8) 

= 9*55 (o c.p.M.J 

It will be apparent that the mathematical study is more easily carried out 
in terms of the phase velocity than in terms of the frequency, and in fact 
the main equations obtained in the study are equations in a> or co* ; owing 
to the close connection between co and F expressed in (5.8), such equations 
are commonly termed “ frequency equations.” 

The value of co is also the ratio of velocity amplitude to displacement 
amplitude, or similarly of acceleration amplitude to velocity amplitude, 
in a sine wave vibration. Tables of values of w* are given in Appendix V, 
page 145. 

6. Two-mass systems with one degree of vibrational 
freedom. 

The torsional system depicted diagrammatically in Fig. 4 is 
an example of a two-mass system having only one degree of 


(J|) FRICTIONLESS (J,) 



Fig. 4.—Two-mass system with one natural frequency. 


freedom so far as vibrations are concerned. Let the angular 
displacements of the inertias and J t from the static position 
be 0 lt 0, in the same sense, and let the torsional stiffness of the 
shaft be C. Then, as the torque T transmitted by the shaft is 
proportional to the twist between the ends, 


JA = - G(d 1 - e t )\ 

i = — G(0t — 0i)J 


( 6 . 1 ) 


These equations can be written in the operator form: 


(J 1 D»+C)0 1 - <70, (o)\ 
(JtD'+oe, = CO i (6)/ 


( 6 . 2 ) 


Equations (6.2) are simultaneous operational equations which 


may tentatively be treated as algebraic equations, so long as 
the result obtained in this way is checked in the original equations. 
The method of treatment is as follows : operate on both sides 
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of (6.2a) with the operator (J i D i +C) which occurs on the left- 
hand side of (6.26). Thus : 

(J t D*+C)(J l D*+C)B 1 = C(JtD 2 +C)6 t 

The right-hand side of this last equation is, by reason of the 
equation (6.26), equal to Hence, 

(J t D 2 +C)(J 1 D i +C')0 1 = C 2 6 1 

Dropping the variable 0 X and expanding the brackets, 

•J 1 J 1 D 4 +(Ji+«/,)£D s = 0 

i.e. D* = 0 or - W£±+£A . . . (6.3) 

J iJ t 

Comparing the solutions (6.3) with equation (2.3) it is evident 
that the two frequencies of vibration are 

F = c.P.s., where a>* = 0 or — —f~Y • (6.4) 

The zero frequency solution corresponds to a steady rotation 
0i = 0* = At+B (for which clearly D 2 — 0 ), while the finite 
solution gives the frequency of free vibration of the system. 
That there is only one degree of vibrational freedom is shown 
by the fact that the amplitudes and phase-angles of the motion 
of the two inertias are related by ( 6 . 1 ). Let the motion of the 
two inertias be 

0 X == -4|.sin ((at-\-<f> x ) 

0 , = 4 s .sin (a >t+<f> t ) 

then 0 , = — a>*4 1 .gin (a>t+(f> 1 ) 

and 0 , = — a> 2 .4 t .sm (cot 

But from ( 6 . 1 ) J X S X = — JJS t 

hence- = — JtA t sm(cot+<l> t ) . (6.5) 

If the amplitudes A x and A t are regarded as positive (see the 
note on the sign convention in section 4) then sin (a)t+<f> t ) = — 1 
whenever sin (oit+<f> 1 ) — 1 . Hence <f> x and <f> t are separated by 
an odd integral multiple of n, for if 

cot+<f> 1 = 2K x n+ji/2, then sin (cot+<f> x ) = 1 
and = (2K t +\)n-\-n/2, K x and K t being integral. 

Thus 

= [2(^t — ^i)+l]7r. 

Furthermore, (6.5) gives the result J x A t = J t A t . If the motion 
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of Ji is -A 1 .sin then that of J t is 

0, = ^i^i.sin (a)t+<f>+n ). . . (6.6) 

J i 

- (S)P( X * - Jr *>+ 1 *'- 

The vibrational motion of the two masses is therefore such that 
they move in opposite directions, and it is important to note 
that there are only two arbitrary constants A 1 and </> in the formula 
(6.6). The initial conditions of the motion can be given as the 
displacements and velocities of the two inertias at the instant 
t = 0, and it is left as an exercise for the reader to include 
appropriate terms in the solution to express the steady rotation 
of the system, and generally to determine the four necessary 
arbitrary constants to make the solution fit the general initial 
conditions. (See Exercise 2 at the end of this chapter.) 

The equivalent linear system would consist of two masses 
connected by a spring and supported on a smooth table, the 
motion being horizontal along the line joining the two masses. 
The simplified system of Fig. 1 is the result of making one such 
mass infinite and rotating the system from a horizontal to a 
vertical position. 

Summary. 

Whenever the physical properties of a mechanical system are 
such that the force acting on a mass tending to restore it to its 
position of equilibrium is proportional to the deflection from that 
position, the mass can execute a sine-wave vibration at a constant 
frequency which depends solely upon the mass and the intensity 
of the restoring force. The formula for the frequency is given 
in (5.1) on page 11. 

The amplitude and phase-angle (relative to the chosen time- 
datum t — 0) depend upon specified initial or other conditions; 
for example, if the initial displacement and velocity are given, 
these constants can be calculated and the vibration is then com¬ 
pletely determined. 

EXERCISES I 

1. Find the general solution of the equation 
x+Px+Qx = 0, 
where P and Q are independent of x and t. 

(See end of section 2.) 
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2. Determine the motion of the two flywheels in Fig. 4 if the initial 

displacements and velocities are G lf <9 2 , & 2 ; and find the conditions 

(i) for steady motion only, with no vibration, and (ii) for vibrational 
motion only, with no steady rotation. 

(Method : show that the steady rotation formula 6 X = 0 t — At-\-B is 
a solution of the equations of motion, and determine these two constants 
and the two constants contained in the formula for vibratory motion, 
section 6, in terms of the initial conditions. Note that in the steady 
rotation there is no relative motion between the flywheels.) 

3. Find the general solution, in terms of exponential functions, of the 
fourth-order equation 


d*y 
dx * 


= fcfy 


(This equation occurs in the theory of beam vibration, see Chapter V, 
section 25. Proceed as in section 2 ; the general solution contains four 
arbitrary constants.) 

4. Assuming that the vibration is in the form of a sine-wave, determine 
the natural frequency of the system in Fig. 1 by considering the energy 
in the system. The kinetic energy is a maximum, and the potential 
energy zero, in the positions where the displacement is zero ; the kinetic 
energy is zero, and the potential energy a maximum, where the displace¬ 
ment is a maximum in either direction. Make use of the fact that the 
total energy content is a constant in the absence of dissipating forces. 
(See page 92.) 

5. If a sinusoidal vertical force is applied to the mass of Fig. 1, calculate 
the amplitude of the force required to maintain unit amplitude of vibration 
of the mass ; hence find the natural frequency, as that frequency at which 
the required force has zero amplitude. 

(The extension of this process is very useful in the study of more 
complex systems ; see page 73.) 

6. Calculate the inertia which must be added to the flywheel in the 
torsional pendulum (Fig. 3d) so that any given frequency may be the 
natural frequency of the system. 

(This inertia, which is a function of the given frequency, is termed 
the “ tuning inertia ” ; it has the same magnitude as, but the opposite 
sign to, the “ effective inertia,” which is widely used in the analysis of 
more complex systems; see page 60.) 


c 



CHAPTER n 


DAMPING AND FORCED VIBRATION 
{One Degree of Freedom — continued) 

Introductory. 

As already noted, the sine-wave motion discussed in the preceding 
chapter is in practice modified by the effects of energy dissipation. 
These effects must be taken into account when studying the 
vibratory response to applied forces, otherwise the over-simplified 
theory leads to results which are not obtained in practice. 

In the first section of this chapter it is shown that a force 
opposing motion, and proportional to the velocity of the motion, 
does a definite amount of work during each cycle of the vibration. 
Such a force, usually termed a “ viscous damping force,” is 
normally assumed to be present in all vibrating systems, and is 
found to account quite well for many of the actual observed 
properties of such systems. The behaviour of the simple system 
of Fig. 1, modified by the inclusion of a “ dashpot ” exerting 
such a force on the mass, is studied in Section 8. 

Generally speaking, the free vibrations of structures are not 
of importance in themselves ; what is of much greater practical 
importance is the response to imposed vibratory forces. For 
various reasons, the imposed forces are at first supposed to be 
sinusoidal. The vibratory force due to an unbalanced rotor is 
of this form, and the practical result of Fourier’s theorem (see 
Chapter VI) is that any type of periodic force can be represented 
as the sum of a number of sinusoidal forces; finally, the response 
to a sinusoidal force is a sinusoidal displacement. Consequently 
the force acting in the system of Fig. 9, page 25, is assumed to 
be of this type. 

The mathematical study of the response to the applied force 
displays the extremely important phenomenon of resonance, which 
may truly be said to be the key to the whole theory and practice 
of vibration treatment. 

As examples of the application of the theory, the use of 
flexible suspensions to reduce the transmission of vibratory 
forces, and the response characteristics of a simple seismic 
vibrograph, are briefly described in Section 12. In the second 
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part of this section the use of complex numbers is illustrated; 
the introduction of these quantities greatly simplifies the mathe¬ 
matical analysis of damped vibrations. The question of insulation 
is further considered in Appendix IV, page 134. 

7. Dissipation of energy. 

It has already been noted that the amplitude of free vibra¬ 
tion of a practical system diminishes in time so that the motion 
decays to zero, and that during this process there is a loss of 
energy. It is now necessary to determine the origin of this 
energy loss. If a system is constructed to resemble that in Fig. 1, 
and it is set in vibratory motion by an initial displacement, the 
motion is observed to be approximately sinusoidal during suc¬ 
cessive small intervals of time ; in fact, it differs from sinusoidal 
motion only in the diminution of amplitude. In order to discover 
the cause of the loss of energy, the motion will be assumed sinu¬ 
soidal and the work done during one cycle of the motion by an 
impressed force will be calculated. 

Let the displacement be given by 

x =X .sin mt .... (7.1) 

and let the force / be a linear function of the displacement, 
velocity and acceleration, say 

f = A+Px+Qx+Rx . . . (7.2) 

This force is applied to the mass, and the work done by the 
force is given by 

W = J/.d* 

dx 

But = x, and so dx = x.dt. Thus 

W = jfx.dt. 

If attention is confined to one cycle of the motion the limits of 
integration differ by the period of vibration. The frequency 
being m/2n, the period is 2n/co, and tl^e work done by the force 
during one cycle is 
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_ 

Thus W = \ (A+Px+Qx+Rx)x.dt. 

Jo 

Giving x its sinusoidal form (7.1) the integral becomes 

2 n 

W = I m [^4+(-P — -Bco 2 )A - .sin cot+QcoX .cos <o<]g>X. cos cot.i 
Jo 

2 n 2 n 

= | AcoX .cos cot.dt+ | \(oX\P — Rod 2 ) sin 2a>t.dt 
Jo Jo 


Pix 

Jo 


“ = o)X 2 Qji 


2 co 2 Q( 1+ cos 2 cot)dt 


. (7.3) 


It is evident that the only part of the force that has done work 
over the whole cycle is the part Qx = QXco cos cot ; the other 
parts of the force each contribute nothing to the work done 
over the cycle. This result suggests that an explanation of the 
dissipation of energy during the motion of the practical system 
may be that work is being done against a damping force propor¬ 
tional to the velocity of the mass. The theoretical results 
obtained by analysing the motion on such an assumption are in 
good agreement with practical results. 


8. Damped motion of single-mass system. 

The theoretical system of Fig. 1 is now extended, as shown 
in Fig. 5, by the addition of a “ dashpot,” the function of which 
is to exert a force cx proportional to the velocity of the mass 
and opposed to the direction of motion. The symbolic spring 
and dashpot do not necessarily represent separate parts of the 
system, but rather different properties of the same part, i.e. the 
practical spring. The symbolic spring represents that property 
of the spring whereby it opposes motion of the mass with a force 
proportional to the displacement, and is merely an energy¬ 
storing device Incapable of dissipation; the constant of pro¬ 
portionality, k , is sometimes referred to as the ncm-dissipative 
stiffness. The symbolic dashpot represents that property of the 
spring (and surrounding air) whereby motion is opposed with 
a force cx proportional to the velocity, and is merely a dissipative 
device incapable of storing energy ; the constant c is called the 
damping coefficient or resistance . 
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Apart from air-resistance, the energy is dissipated by way 
of internal friction in the spring; in the system of Fig. 5 the 
dashpot represents both the resistance of the spring and that of 
the air surrounding the system. The force exerted by the dash- 



Fig. 5.—System of Fig. 1 modified by inclusion of a dashpot. 

pot on the mass is — cx, and the equation of motion is therefore 

mx = — lex — cx . . . (8.1) 

As before, let to 2 = k/m, and further put 2 y = c/m. (8.1) then 
becomes 

x+2yx-\-a) 2 x = 0, 

and the operational equation obtained by using the operator D 
and dropping the variable x is 

D 2 +2yD+co 2 = 0 . . (8.2) 

Proceeding on the lines of the previous work (Section 2) this 
equation is solved for D, thus : 

D — - y ± Vy 2 - <o 2 

and in this solution a and take the place of ± ion in (2.4). 
The form of (2.5) suggests the solution 

x = C^+C^ 1 
i.e. 

x=Cj. exp (-y+Vy 2 —(o 2 )t-{-C i exp (—y—Vy 2 —(o 2 )t .. (8.3) 

where the notation exp Q is used in place of e® for convenience. 
Direct substitution shows that the function (8.3) does in fact 
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satisfy the equation of motion (8.1), and as it contains the neces¬ 
sary two arbitrary constants it is the required general solution. 

The physical significance of the solution (8.3) depends upon 
the nature of the expression Vy 2 — co 2 . Three cases need be 
considered, according to whether y is greater than, equal to, or 
less than co. 

(i) If y > co, both roots a and /3 are real and negative, for 
Vy 2 — co 2 is then real and numerically less than y ; the graph 
of the displacement plotted against time is therefore the sum of 



Fig. 6.—Displacement: time curve for damped motion (more than critical 

damping). 


two decreasing exponential curves (see Appendix I, section 31). 
The shape of the graph depends upon the values of G x and C t , 
which in turn depend on the initial conditions of the motion. 
Fig. 6 illustrates some typical shapes ; in this diagram the 
function x/C x is plotted against time for the particular values 
a = — 1, p = — 3, and for various values of the ratio G t /C 1 . 
(ii) If y = co, both roots a and are equal to y. (8.3) becomes 

* - (Ct+CJe-*, 

but this function involves only one arbitrary constant (Ci+C t ) 
and cannot be the most general solution. To find the general 
solution the following procedure is adopted: 
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Let x = 2 .exp (— yt), where 2 is to be determined. Then 
Dx = — ye^.z+er^.Dz 
and D 2 x = y*eT vt .z — 2ye~ rt .Dz+e~ rt .D*z. 

Substitution of these values in ( 8 . 2 ) leads to the result 
(to 2 — y 2 )e-*.z+e- , *.D 2 2 = 0 
i.e. D 2 z = 0, as y = to. 

Double integration with respect to time gives : 
z = C a t-\-C t , 

and x = ( C i t+C i )e~ vt . . . (8.4) 

The displacement-graph plotted against time may have either 
of the forms indicated in Fig. 7, depending on the initial con¬ 
ditions. (See Exercise 1 at the end of this chapter.) 


Fig. 



7.—Displacement: time curve for damped motion (critical damping). 


The condition y = co is known as the condition of critical 
damping ; the corresponding critical value of the damping 
coefficient (c = c 0 ) is a convenient standard of comparison for 
damping coefficients and is easily found: 



and c = c„ = 2 ray = 2 ma> 

whence c 0 = 2Vmk ..... (8.5) 

The physical significance of the critical-damping condition is 
that in this condition the system is brought practically to rest ;, 
in a shorter time than if the damping is either greater or less 
than the critical value (see Exercise 4 at the end of the chapter). 
For this reason, many indicating instruments have their suspen¬ 
sions critically damped. 

(iii) If y < to, let q = Vo> 2 — y 2 , so that q is real. 

Then x = e-*(C l e* t +C t e- i ° t ), 
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and as in Section 2 of Chapter I this expression may be rewritten 
in trigonometric form: 

x = e~ vl (A . sin qt +B . cos qt)' 
where A = i(C x — C t ) 

B = C x +C t 


or alternatively, 


( 8 . 6 ) 


x = e y *.o.sin ( qt+<f>) 
where a = VA 2 +B 2 

tan <f> = B/A 

The graph of the displacement plotted against time has the 
form indicated in Fig. 8, the curve being a damped sine-wave 



Fig. 8.—Displacement: time curve for damped motion (less than critical 

damping). 


lying between a decreasing exponential curve AB and its mirror- 
image A'B’ in the time-axis. 

Case (iii), when the damping coefficient c is less than the 
critical value e 0 , is the most important case in practice. A useful 
property of the curve (Fig. 8) is the manner in which the wave 
decreases in amplitude. The points P x , P t are successive maxima 
at times t x and t t , so that t t = t x +2n/q ; the amplitudes x x and x t 
at these times are very nearly equal to the values of the expon¬ 
ential function a.exp (— yt) at P x and P t , i.e. 


x x 

Xt 


= o.e -y< *j 


approximately. 


Let the ratio x t /x x be [i. Then fi — exp [— y(t t — $0] 


i.e. n = exp (— 2 yn/q) = exp (— jic/mq) 

This ratio is independent of the initial amplitude and the time 
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h, and is thus the same for all pairs of consecutive maxima. 
If A — — log« n 


then A = — . . . . (8.7) 

mq 

This quantity is termed the logarithmic decrement, and provides 
a ready means of determining the damping coefficient c of a 
system, as all the quantities in (8.7) except c can be evaluated 
experimentally from the physical constants of the system and 
a displacement curve or vibration record. 

There is no real cycle of the motion, as the amplitude of 
vibration is constantly decreasing, but the maxima and minima 
of the displacement curve (Fig. 8) each recur after successive 
equal intervals of time 2n/q, and the quantity q/2n is termed the 
“ frequency ” of vibration. The larger the value of (y = c/2 m) 
is made, the lower the frequency 




2n 2 ji 


becomes. 

As a practical result it is interesting to note that after a>/y 
“ cycles ” the amplitude of the motion is only 0-19 per cent, of 
its original value, so that in most practical examples the motion 
is effectively damped-out after this number of cycles (see 
Exercise 5 at the end of the chapter). 


9. Forced motion. 

So far the analysis has been concerned with free vibrations, 
i.e. vibrations maintained solely by the energy stored in the 
system at the commencement of the motion. If the vibration 


Fig. 9.— System of Fig. 5 
modified by inclusion of an 
applied force. 


of a damped system is to be maintained at constant amplitude 
there must be an energy input. Fig. 9 represents the system 
of Fig. 5 acted on by a harmonic force P = P 0 sin pt. The 
equation of motion is 

mx — P 0 sin pt — lex — cx 





P 0 sin pt 


l ™ I 
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i.e. x+2yx-\-(o 2 x — — sin pt 1 

to \ . (9.1) 

where 2 y = c/m, to* = k/m J 

In operator form this becomes 

[D 2 -\-2yD+a) 2 ~\x = — sin pt. 

m 

Operating on both sides of the equation with [D 2 — 2 yD-f<u 2 ] 


[(D 2 +«o 2 ) 2 — 4y 2 D 2 ]a; = —[D 2 — 2yD+oo 2 ] sin pt 

m 


and this last equation can tentatively be written as : 

P. [D 2 -2 yD+(u 2 ] . ^ 

* to ■ [(D 2 +to 2 ) 2 - 4y 2 D 2 ] Sm P 

or x = —[D 2 — 2yD+a> 2 ][(D 2 +a ) 2 ) 2 — 4y 2 D 2 ) -1 sin pt ( 9 . 2 ) 

TO 

If the second bracket, with its negative index, is now expanded 
by the Binomial Theorem, it is seen that every term in the 
expansion contains powers of D 2 . But 

D 2 sin pt = — p 2 sin pt 
D* sin pt = p* sin pt, etc., 

so that wherever the operator D 2 appears in (9.2) the multiplier 
— p 2 can be substituted. (It is. to be noted that the validity 
of this procedure is to be checked by ascertaining whether the 
derived solution satisfies the original equation.) 9.2 now becomes 

_ (P„/to)[co 2 — p 2 — 2 yZ>] sin pt 
(a) 2 — p 2 ) 2 -f4y 2 p 2 

and performing the operation indicated by the square bracket 
the result is obtained: 

* = -*•> 8in ^ - 2 yp.cosp<] 

which can be written as: 


x = 6 .sin (pt — y>) 


b 

tan y> 


P,/to 


[(co 2 - p 2 ) 2 +4y 2 p 2 ]* 

_2yp 


ft) 1 


P * 


(9.3) 
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Using the equations 2 y = c/m, and c 0 = 2Vmk, ( 9 . 3 ) can be 



This solution contains no arbitrary constants and cannot 
therefore be the general solution. It is termed the particular 
integral of the equation (9.1). To determine the general solution 
the procedure is as follows : 

Let x = u+v, where v stands for the solution (9.3). Substitu¬ 
ting in the operator form of (9.1), 

[Z> 2 + 2 yZ)+ft>*](tt+«;) = — sin pi 

m 

P 

i.e. [Z) 2 + 2 yZ)+a> 2 ]i 4 + [-D 2 + 2 y 2 )+cy 2 ]v = — sin^tf 

m 

But v satisfies (9.1) and hence 

p 

[D a + 2 yD+ft) 2 ]v = — sin pi 
m 

thus [D 2 + 2 yZ)+co 2 ]tt = 0 . 

This last equation is exactly similar to that discussed in Section 8 , 
and if y < co the solution is given by ( 8 . 6 ), i.e. 

u = e _v< .a.sin (qt+<f>) . . . ( 9 . 6 ) 

The function u contains two arbitrary constants (a and <f>) and 
is itself the solution of the equation obtained by putting zero 
in place of the right-hand side of (9.1). Such a function is 
termed a complementary function, and it is a general result in 
the theory of differential equations that the complete solution 
of an equation such as (9.1) is the sum of the complementary 
function and a particular integral determined as in (9.3). (See 
Exercise 2 at the end of this chapter.) 

The complete solution x is given by 

x = u+v = 6.sin (pt — y)+e -y< .o.sin (qt+</>) . (9.6) 
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The term due to the complementary function u is a decreasing 
time-function, as in Fig. 8 ; after a time this part of the motion 
is practically damped out, and effectively the displacement 
function consists merely of the particular integral v. Such a 
motion is termed “ steady-state forced vibration,” and is the 
only part of the complete motion that is at all important. 

10. Resonance. 

The frequency of the steady-state motion is the same as that 
of the impressed force, namely p/2n. It is of interest to investi¬ 
gate the relation between the amplitude of steady-state motion 
and the frequency. From (9.3) the amplitude is given by 

6 = -PoA» 

[(co 2 - p 2 ) 2 +4y 2 P 2 ]* 

If P is independent of the frequency, b has a stationary value 
whenever 

B = (co 2 — p 2 ) 2 +4y 2 p 2 

has a stationary value (maximum or minimum). Differentiating 
B with respect to p, which is proportional to the frequency of 
the impressed force, 

= 8 y 2 p-f 2 (<o 2 — p 2 )(— 2p) 

= 4[2y 2 p-f-p(p 2 — co 2 )]. 

Hence b has a stationary value when p — 0 

or 2 y 2 -fp 2 — co 2 = 0 

i.e. p = V(o 2 — 2 y 2 . . . ( 10 . 1 ) 

It can be shown that when p = 0 , b is a minimum, and for the 
value ( 10 . 1 ) 6 is a maximum ; thus the amplitude of the steady- 
state motion has a maximum at the frequency 

= jr-V'co 2 ~ 2y* 

2 71 271 r 

and the graph of displacement-amplitude plotted against fre¬ 
quency has the general form indicated in Fig. 10 . 

The condition of maximum displacement-amplitude for a given 
impressed force is termed resonance, and the corresponding fre¬ 
quency is termed the resonant frequency. The practical impor¬ 
tance of resonance cannot be emphasised too strongly. Indeed, 
the treatment of vibration problems in engineering is almost 
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entirely a matter of detecting resonant frequencies and avoiding 
resonant conditions, for at resonance the vibration-amplitude of 
the system is maximum and the resulting stresses in the system 
similarly have maximum values. The normal procedure in 
dealing with practical vibration problems is to determine the 
source of vibration (i.e. the origin of the input force) and the 
resonant frequencies, either theoretically or practically, and then 
to arrange matters so that the machine does not run in a resonant 
condition. There are usually four ways in which this may be 
done, any or all of which may be possible in a particular case : 
(i) eliminating the source of vibration, (ii) absorbing the input 



energy in some way, (iii) altering the speed-range of the machine 
to avoid resonance, and (iv) removing the resonance from the 
operating-range of speed by suitable alteration of the dynamic 
constants of the system. The precise technique involved is a 
matter of engineering detail and as such is outside the scope of 
this work, but one type of energy-absorber is discussed in 
Chapter III, section 15. 

A measure of the intensity of the resonance effect is afforded 
by the quantity known as the dynamic magnifier . The dynamic 
magnifier at a point in a vibrating system is defined for a par¬ 
ticular frequency as the ratio of two deflections : (i) the dynamic 
amplitude at that frequency (as given, for example, by b in (9.3)) 
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and (ii) the deflection caused by the application of a static force 
equal in magnitude to the amplitude of the forcing input. In 
the system of Fig. 9 the dynamic forcing input is P, sin pt, and 
the static deflection under the action of a force P„ is P 0 /k at 
the mass. The dynamic magnifier at the mass is thus given 
from (9.3) by 


as 



co 2 = k/m 


a) 2 

p a ) a +4yW ' 


or 


M = 




+4- 




c 0 a eo a 


using the relations 2 y = c/m, c 0 = 2 Vmk. 


(10.2a) 

(10.26) 


Using the equation (10.1) it can easily be shown that at resonance 
this expression becomes : 



where c„ = 2Vmk as before. 

The greater the ratio of the damping coefficient c to its critical 
value c 0 , the less becomes the value of the dynamic magnifier 
at resonance. Fig. 11a represents the dynamic magnifier curves 
for the one-mass system for two values of the damping coefficient: 
10% and 20% of the critical value. The magnification of 
amplitude due to the resonant condition is seen to be very 
much greater for small damping forces than it is for large damp¬ 
ing forces. If there is no damping, the dynamic magnifier at 
resonance has an infinite value, but this is of course a purely 
theoretical consideration; the curve for zero damping is drawn 
dotted in Fig. 11a, and it is evident that for a lightly damped 
system the magnifier is not very different from that for zero 
damping, except in the neighbourhood of the resonant frequency. 
A table of values of the dynamic magnifier in the absence of 
damping forces is given in Appendix 'V, page 146. 

The general shape of the dynamic magnifier curve Fig. 11a 
is typical for single-mass systems with a forcing input of the 
form P = P„ sin pt, where P 0 is constant; this type of forcing 
is commonly termed “ constant excitation.” A type of excitation 
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that is of very frequent occurrence in engineering is that due to 
centrifugal action on unbalanced rotors; the reactions at the 
bearings supporting the shaft are then of the form P=Qp 2 sin pt, 
as the centrifugal force acting on the unbalanced rotor is propor¬ 
tional to the square of the rotational speed. The dynamic 
magnifier curve for this type of forcing, commonly termed 




tion ; ( d) inertia excitation. 

“ inertia excitation,” has the form indicated in Fig. 116, from 
which it can be seen that, for inertia excitation, the magnifier 
has zero value at zero frequency and approaches unit value at 
very high frequencies. (Note ,—In this case the relevant static 
deflection to be used in the definition of the dynamic magnifier 
—see the paragraph immediately preceding equations 10.2—is 
the deflection Q/k under the action of a static force Q, not Qp % 
which has zero value at zero frequency.) 
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Although the free vibration, given by the complementary 
function (9.5), is damped-out after a short time owing to the 
decreasing exponential factor, its amplitude may be quite con¬ 
siderable at the beginning of the motion, the amplitude of the 
total motion (9.6) being thereby increased temporarily to an 
extent which can be troublesome in practice; this effect is 
particularly significant in cases where the exciting force is started 
and stopped frequently, when the “ starting transient ” (i.e. the 
complementary function) may seriously increase the risk of 
fractures being caused by excessive vibration amplitudes. Many 
instruments are critically-damped to avoid this vicious condition. 
An elementary treatment of transients has been given elsewhere 
by the author, see references 13 and 14 in the Bibliography, 
page 151. (See Exercises 6-9 at the end of the chapter.) 

11. Resonance and natural frequency. 

In the present section the resonance characteristics of a 
single-mass system under the action of a constant excitation 
force are discussed. Analogous results are easily obtainable for 
the case of inertia excitation, the most important practical 
characteristic of which is concerned with the problem of insulation, 
see Section 12, Example I. 

The resonant frequency of the single-mass system of Fig. 9 
is, from (10.1), 

r ~-i- ■ ■ (111 > 

and this expression can be rewritten in terms of the critical 
damping coefficient (c„ = 2 Vmk) as 



It can be seen from (11.2) that if the damping force is small so 
that the square of the ratio c/c t can be neglected in comparison 
with unity, the resonant frequency is equal to cu/2 n. This 
frequeitcy is the natural frequency for zero damping (Chapter I), 
and if the system is lightly damped no great error is involved 
in neglecting the damping force altogether and making use of 
the natural frequency of the resulting system instead of the true 
resonant frequency. 

It is shown in Chapter IV, section 18, that this result can 
be extended to multi-mass systems, for such systems vibrating 
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at any frequency can be considered as built-up of partial single- 
degree-of-freedom systems, plus an auxiliary mass known as the 
effective inertia of the complete system ; and the resonant frequency 
of each such partial system is approximately equal to the natural 
frequency for zero damping, so long as the damping forces are 
small. Some examples of the error involved are listed in Table II 
for various values of the ratio c/c 0 . 


Table II 


Damping coefif. 

Natural frequency 


Critical value 

for zero damping 
Resonant frequency 

Error 

c/c # 

<o/p 

0/ 

/o 

0 

100 

0 

01 

101 

1 

1 0*2 

104 

4 

j 0-3 

! iio 

10 

! 0-4 

1-21 

21 

0-5 

1-41 

i 

41 

The error quoted is that involved in using the natural 

frequency co/2ti instead of the resonant frequency p/2n. 


The practical importance of this result is that it is very much 
easier to neglect damping forces altogether, and calculate the 
natural frequencies of the resulting modified system, than it is 
to calculate the resonant frequencies ; Table II shows that the 
error involved in this procedure is negligible in the case of lightly 
damped systems. 

One further useful property of the resonant condition of 
lightly damped systems is the value of the phase-angle y> in the 
expression (9.3) or (9.4). When p = a>, i.e. at the natural fre¬ 
quency for zero damping, the angle ip has the value n/2, as the 



Fro. 12._90° phase-angle between force and displacement waves at resonance. 

D 
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denominator of the expression for the tangent is then zero, and 
the tangent has an infinite value. Thus in a lightly damped 
system the phase-angle at resonance is approximately jr/2. This 
angle is the angle by which the displacement-wave b. sin {pt — ip) 
lags behind the force-wave P 0 sin pt, as illustrated in Fig. 12 . 
Fig. 13 shows representative curves of the phase-angle for various 
values of the ratio c/c„, plotted against the non-dimensional ratio 
p/a) which is proportional to the frequency p/2n. For very small 
values of the ratio c/c 0 the curve approximates to the discontinu¬ 
ous line for c = 0, the characteristic of which is that for excitation 
frequencies below the natural frequency the displacement is in 



phase with the force, whereas for frequencies above the natural 
frequency the two waves are 180° out of phase so that the rotating 
vectors in Fig. 12 are separated by 180° or n ra dians 

A fairly simple relation between the dynamic ma.gnifW M 
and the phase-angle <f> can be obtained from equations (9.4). 
This relation is 


where 


M = 


cos <f> 

1 - X 2 


X = p/a) = 


freq uency of vibration ’ 
natural frequency . 


(11.3) 


12. Practical examples. 

In order to illustrate the use of the foregoing analysis, two 
practical examples are discussed in this section. The first involves 
the principles of vibration insulation, whereby the force trans¬ 
mitted from a vibrating machine to its support is reduced by 
means of an elastic mounting; the second is the case of a seismic 
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vibrograph for recording vibration-amplitudes, and in the treat¬ 
ment of this problem the complex-number method of calculation 
is introduced. 


I. Insulation (see also Appendix IV, page 134). 

A machine, which may be considered as a rigid mass of 
1,930 lbs., is acted on by an alternating vertical force, whose 
amplitude is 400 lbs., at a frequency of 3,000 c.p.m. The damping 
forces in the system can be made to have 20% of the critical 
value. If the alternating load is taken equally by four spring 
feet of identical design, what must the stiffness of each mounting 
be in order that the total alternating force transmitted to the 
base is only 100 lbs. maximum ? 

In essentials the system is that of Fig. 9 inverted, the four 
identical spring mounts being considered as combined together. 
In Fig. 9 the total alternating reaction on the support is kx+cx=f r 
and if x = b . sin (pt — y>) the force is 

f r = b[k. sin (pt — ip) +cp. cos (pt — y>)] 

= bVk 2 +c 2 p 2 .sai (pt- f-e) 

where the phase-angle e is not important. The maximum value 
of the reaction force is 


fmax. = bVk*+C 2 p 2 


which can be written in a form similar to that used in (9.4) as :. 


where 

and 


fmax. — bm(o 2 Vl + 4A/f = 2 bmoi 2 V\-\-Xp 



p = (pj = (0-2) 2 = 0-04 


i.e. f max , = 2bnuo 2 V 0-25+0-04A. 

Substituting for 6 from (9.4), 

2P 0 VO-25+004A 


fmax. — 


[(1 - A) a +0-16A]* 


and when P = 400 lbs., f max . must be 100 lbs. 
Hence 8V0-25+0-04A = V(1 - A) a +0-16A 

Solving for A, 

= A — 6-66 



i.e. 


£ = 2-58 

(O 
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But it is given that p/2n — 3,000/60 = 50 o.p.s. 
Hence co = 50 x 2n 4- 2-58 = 122 


Also 

Hence 


x> — ^ /—, and m = 1,930 4- 386 = 5 slugs. 
V m 

Jc = 5 x 122 2 = 7-45 x 10 4 lbs./in. 


The total stiffness of the four similar mounts is thus 7-45 x 10 4 
lbs./in., and that of each mount is one-quarter of this value, 
i.e. 1*86 x 10 4 lbs./in. 

This example shows the property of partial insulation possessed 
by a spring mounting ; by choosing the mounting so that the 
working frequency is 2-58 times the natural frequency, the force 
transmitted to the base is reduced to a quarter of the applied 
force. 

In connection with this insulation effect, it is a matter of 
some practical importance to determine the precise result of 
increasing the damping forces in the systems. By the method 
given above the general formula for the transmitted force due 
to an applied force P = P„ sin pt is obtained as : 


f _ P 0 Vl-\-4:^p 

^ max. — / - - ■ -- • 

V(1 — A) 2 + 4A/« 


• ( 12 . 1 ) 


By differentiating this expression with respect to p, keeping A 
constant, it is found that for A less than 2, i.e. for a forcing 
frequency less than V2 times the undamped natural frequency, 
an increase in the damping force causes a decrease in the trans¬ 
mitted force, whereas for A greater than 2, i.e. for a forcing 
frequency greater than V 2 times the undamped natural frequency, 
an increase in the damping force causes an increase in the trans¬ 
mitted force. Representative curves are shown in Fig. 11c. 

If Qp 2 is substituted for P 0 in (12.1), so that the expression 
corresponds to the case of inertia-excitation, the increase in 
transmitted force due to an increase in the damping force, when 
the forcing frequency is greater than V 2 times the undamp ed 
natural frequency, is more violent than for the case of constant 
excitation ; representative curves for this type of excitation are 
shown in Fig. lid. 

The choice between great or small values of the damping 
coefficient in any particular practical case depends to a large 
extent on whether there is any likelihood of the system being 
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vibrated accidentally in a resonant condition. For high fre¬ 
quencies of excitation, large damping forces are a disadvantage, 
but in the neighbourhood of resonance the effect of decreasing 
the damping forces may easily be disastrous, owing to the much 
greater effect of magnification at resonance. 

II. Seismic vibrograph. 

Fig. 14 illustrates the seismic element of a simple type of 
vibrograph. The base A is attached to the structure whose 
vibration-amplitude it is required to measure, and means are 
included for recording the amplitude of the relative displacement 
between A and the mass m. It is required to draw the “ response 
curve ” of the instrument, i.e. a curve showing how the ratio 
amplitude of relative motion 
amplitude of motion of base 

varies with the frequency. The natural frequency of the instru¬ 
ment is 600 cycles per minute, and the damping can be made 
to be between 0-1 and 0-7 of the critical value. 



Fig. 14.—Seismic element of simple vibrograph (recording-gear not shown). 

Let the displacements of the base and the mass be x and y 
respectively, x being a sinusoidal time-function of frequency 
p/2n cycles per second. The alternating forces acting on the 
mass are k(x — y) and c(x — y), and the equation of motion is 

k{x - y)+c(x — y) = my 

i.e. my+cy+ky = cx+kx . . . (i) 

Let the amplitude of x be X, so that 

x = X . sin pt. 

It is evident that the right-hand side of the equation (i) represents 
a , sinusoidal force of frequency p/2ji, and the displacement y 
can be obtained by a method similar to that adopted in Section 9. 
It is convenient here, however, to introduce an alternative 
method of solution involving complex-numbers. 

It has been shown in Section 3 that a sine-function of fire- 
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quency F can be regarded as the projection of a vector rotating 
with angular velocity 2nF (see Fig. 2). If the displacement 
x — X.ain.pt, the velocity % — Xp. cos pt, or 

x = Xp. sin (pt+Tj^y 

and the corresponding vector is phased n/2 radians ahead 
of the displacement vector. Elementary vector-analysis (see 
Appendix I, section 35) utilises the symbol j to represent a 
versor-operator which turns the associated vector through an 
angle n/2 radians in the positive sense. If the vector resulting 
from this rotation is operated on again by j it is turned through 
a further angle n/2 radians, so that it is finally phased n radians 
ahead of the original vector. As 

sin ( 0+n) = — sin 0, 

the operation j 2 is equivalent to changing the sign, i.e. 

j 2 = - 1 

and in general the operator j has the algebraic properties of 

i = V — 1. 

The form of equation (i) shows that the displacement y is 
a sinusoidal time-function of frequency p/2n. Thus the velocity 
y can be written jpy and the acceleration y — — p 2 y. (As a 
check on the method, it will be observed that in Section 2 it was 
found that the velocity x of a displacement x = X.sin cot is 
x = Dx = icox (see 2.4) 

and in the present section the operator j is being used in place 
of the number i.) 

Rewriting the equation of motion in terms of j, 

(k - mp 2 +jcp)y = (k+jcp)x. 

Hence - — -—- 

x k — mp 2 -\-jcp 

and V — x = _ m P 2 __ mp 2 {k — mp 2 — jcp) 

x k — mp 2 +jcp ~ (k — mp 2 ) 2 +c 2 p 2 
Substituting the sinusoidal form for x , 

V — x __ mp 2 [{k — mp 2 ) sin pt — cp . cos pt] 

X (k — mp 2 ) 2 +c 2 p 2 

cp 

k — mp *’ 


where tan <j> ~ 
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The amplitude of the relative motion is therefore 


rel. amp. = 


mp 2 X 

[(k — mp 2 ) 2 +c 2 p 2 ] 


or, making use of the substitutions co 2 = k/m, 
A = (p/co) 2 and p = (c/c 0 ) 2 , 


rel. amp. = 


XX 


t(l - A) 2 +4A /M ]* 


Cq 


2Vmk, 


Fig. 15 shows curves of the relative amplitude plotted against 
frequency (it is given that the natural frequency is 500 c.p.m.) 
for values of the damping coefficient from 0-1 to 0*7 of the critical 
value. For frequencies above twice the natural frequency, i.e. 



i a p/to 

Fio. 15.—Response-curve for Fig. 14. 


above 1,000 c.p.m., the response of the instrument is within 
2^% of the motion of the base if the damping has 70% of the 
critical value. At such frequencies the mass m is practically 
stationary in the gravitational field of the earth, and the recorded 
relative motion is therefore almost exactly the same as the 
absolute motion of the base in the gravitational field. 

Although it is no part of this present work to enter into 
engineering details, it is relevant to observe here that large 
damping fordes can be introduced into the system deliberately 
by means of an electro-magnetic device, the damping effect of 
which can be varied over a wide range by adjustment of the 
electric current flowing through the coil of an electromagnet. 
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Summary. 

Provided the damping force is less than a certain critical 
value, the effect of such a force on the free vibration of a simple 
system is to cause a decay in the amplitude. Mathematically 
the decay is expressed by multiplying the sine function by a 
decreasing exponential function. If the damping force exceeds 
the critical value there is no true vibration, but only a slow 
return to the position of rest after perhaps one swing. 

The response to an applied sinusoidal force is itself sinusoidal, 
and has the same frequency as the force. The amplitude depends 
upon the ratio of the frequency of force to the natural frequency, 
being greatest when this ratio is unity. The frequency cor¬ 
responding to a maximum response to a given amplitude of 
excitation force is termed the “ resonant frequency ” ; for small 
damping forces such as are commonly encountered, the resonant 
and natural frequencies are very nearly equal (for example, 
Table II on page 33 shows that if the damping force is less than 
10% of its critical value, the difference between the two frequencies 
is less than 1%). 

At resonance the amplitude of response to an applied vibratory 
force is very much greater than the response to a static force 
of the same magnitude. This “ dynamic magnification ” at 
resonance varies with the intensity of the damping force, as shown 
by equation (10.3) on page 30. In practice the damping force 
may be as small as 1% of its critical value, corresponding to a 
dynamic magnification of 50 : 1 at resonance. 

EXERCISES II 

1. Determine the condition distinguishing between the two forms of 
displacement curve illustrated in Fig. 7 and derived from equation (8.4), 
and show that the curve cannot cross the time-axis more than once. 

(Method: the type of motion depends upon the initial displacement 
and velocity; solve to find the time at which x = 0 in terms of these 
initial values.) 

2. Extend the method of Section 9 to prove the general theorem: 

If bp{D)]x = [a 0 +a 1 D-{-a 2 D 2 -f. . . +a n D n ]x, where the a’s are 

constants, and x = v is a particular integral of the equation 

[y,(D)> = /(*) 

then x = u+v is the general solution of this equation, where x = u is 
the general solution of the equation 

[y>mx - 0 . 
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Show also that this result is not true for non-linear equations such as 

z* = ax+bt. 

3. Show that if the force input in Fig. 9 is of the form P = Qp* sin pt, 
where Q is a constant, the resonant frequency formula analogous to 
(10.1) is 

o* 

p ~ [«»* - 2 

Draw a curve showing the relation between the displacement amplitude 
and the frequency of the applied force. Show that for both types of 
motion under the action of forces of the form P *= P 0 sin pt and P = Qp* 
sin pt this curve has no maximum if the damping coefficient has a value 
greater than 1/V2 times the critical value. 

(This type of motion occurs when the exciting force is the result of 
centrifugal action on an unbalanced rotor.) 

4. Verify that for damped free motion of a single-mass system the 
mass is effectively brought to rest in a shorter time if the damping is 
critical than if the damping is either greater or less than critical. 

(Consider a particular system in which the damping has the values 
1*1, 1-0, and 0*9 times the critical value, and the initial displacement 
and velocity are x Q and zero respectively. Find the times taken for the 
displacement (or in the case of 90% critical damping, the exponential 
“ envelope ” AB in Fig. 8) to fall to 10% of the initial value.) 

5. Show that when the damping is very much less than critical, the 
amplitude of the exponential envelope AB in Fig. 8 is only 0*19% of its 
initial value after a time 2 n/y, i.e. after co/y “ cycles ” ; and show further 
that after a time l/y the amplitude is 1/e or 36*8% of its initial value. 

6. The mass in Fig. 9 is held rigidly while a harmonic force P = P 0 
sin pt is applied ; it is released at an instant when sin pt — 1 ; if p = m /10, 
and y — co/5, draw a curve of displacement against time. 

(Let the amplitude of the forced motion be b as in (9.4); at the instant 
of release the displacement due to the forcing must have this maximum 
value, so that the displacement due to the free motion, at a frequency 
o/2 n, must be-6; similarly the velocity of the free motion is zero. 
The complementary function (9.5) therefore has an initial amplitude 6 
and a phase angle 3rc/2. Sufficient is now known to enable the required 
curve to be drawn. It will be seen that shortly after the instant of release 
the displacement is instantaneously nearly double the amplitude of the 
steady-state motion.) 

7. A load P is suddenly applied to the mass in Fig. 9, the system 
being in a state of rest prior to the imposition of the load. If the damping 
force is 5% of its critical value, calculate the response, and show that the 
mass finally reaches a static position at a deflection P/k below its normal 
position. (See reference 13.) 

8. The load P (see previous question) is suddenly removed after an 
interval T. Find how the free vibration after the load is removed depends 
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upon the relation between T and the period of free vibration. (See 
reference 13.) 

9. A constant vertical acceleration a is given suddenly to the base 
A in Fig 14, and ceases suddenly after an interval T. Investigate the 
response of the mass m. (See reference 14.) 



CHAPTER IH 

UNDAMPED MOTION WITH TWO DEGREES OF 

FREEDOM 


Introductory. 

In the previous chapters attention has been confined to the 
systems possessing one degree of vibrational freedom only : that 
is, to systems whose configurations can each be completely 
specified by means of one variable quantity, which is conveniently 
taken to be the displacement of the single mass from the position 
of rest. 

A more complicated system may require more than one 
variable for a specification of its configuration, and it is found 
that it can have more than one natural frequency. Similarly it 
is found that the response can have a maximum amplitude at 
more than one frequency, i.e. there are two or more resonances. 

It is shown in Chapter IV that provided the damping forces 
are reasonably small, the resonant frequencies of complicated 
systems are very nearly the same as the natural frequencies 
calculated as if damping forces were absent. So nearly identical 
are these two sets of frequencies that for all practical purposes 
it is sufficient to calculate the natural frequencies, assuming no 
damping. As outlined in Appendix IV, page 134, the practical 
treatment of vibration problems depends largely upon the avoid¬ 
ance of resonance, and so it is important to be able to determine 
the resonant frequencies ; it is much easier to calculate natural 
frequencies. If the damping forces are large, however, they must 
be taken into account, in the manner outlined in Chapter IV, 
page 76. 

For the systems discussed in this chapter there are two natural 
frequencies. It is found that the vibrating patterns are not the 
same at the two natural frequencies; the values of the two 
displacement amplitudes involved have different ratios at the 
two frequencies. It is well to keep in mind the fact that there 
must be some fundamental physical distinction between modes 
of free vibration at different frequencies; in some cases this 
distinction is to be found in the number of nodes, or points where 
there is no vibration. 
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13. Two-mass systems. 

In this chapter two different types of system with two degrees 
of freedom are dismissed; the first is typified by the system of 
Fig. 16, in which each of the two masses has one degree of free¬ 
dom—namely, motion in the vertical direction—and Fig. 21 
illustrates the second type, in which one mass has two degrees 
of freedom as denoted by the arrow diagram in this figure. 

The two types of system are fundamentally 
different and are treated in order. 

It has been shown in Section 11 that the 
resonant frequency of a lightly-damped system 
with a single degree of freedom is practically 
equal to the natural frequency calculated on the 
assumption of zero damping. This result will be 
extended to cover multi-mass systems in Chapter 
IV, section 18, and present attention is confined 
Fig 16 —Two- determination of the natural frequencies 

mass system with of the two-mass system in Fig. 16. Let the 
quencTes 111 ^ ^ di s pl acemen f s of the two masses below the posi¬ 
tion of static equilibrium be x x and x 2 . The 
extension of the lower spring is then x 2 — x 1} and the equations 
of motion for the two masses are : 



i.e. 

and 


m x x i = — k 1 x x +k 2 {x 2 — x x ) 
m x x x = — (k 1 +k 2 )x 1 +k^cS 
mjc 2 = k 2 pc 1 — k^x 2 j 


. (13.1) 


The spring stiffnesses are constant arid the motion of the system 
can therefore be assumed to be sinusoidal; the validity of the 
assumption is, as usual, to be checked by investigating whether 
the solution obtained does in fact satisfy the equations of motion. 
Without loss of essential generality the phase-angle for the dis¬ 
placement of the mass m l can be taken as zero, and the displace¬ 
ments can be written: 


X\ — (t\ sin COjt 
x t — a s sin ((o 2 t+cf>) 

Substitution of these values in (13.1) gives the equations : 

(— m 1 m 1 i +k 1 +k^)a lL sin co x t = k^n t cos <f> sin a> t t 

+&8«s sin <j> cos co t t . (a)' 

(— m a ft> s 2 +& a )aj cos <f> sin ojt 

+( — m,ft> g 2 +& 4 )a i! sin <f> cos co t t = lc& x sin a> x t (b) 



(13.2) 


(13.3) 
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These equations are true for all instants of time; put t = 0, 
then (13.3a) becomes 

k& t sin <f> = 0. 

Now, k a is not zero by hypothesis, and hence either a* or sin <£ is 
zero. If a 2 is zero, so is a x from (13.36), but if both a x and a % 
are zero there is no motion. Thus sin <f> = 0, and equations 
(13.3) become 

(— mxtOi 2 +&i+&s)«i sin co x t = k^a a cos <f> sin mj, (a) 

(— m^o a 2 +k a )a a cos <f> sin mjt = k ^ sin co x t (6) 

(where cos <f> = ± 1). 

The ratio a x /a 2 obtained from (13.4a) is 

a 1 _ k a cos <f> sin co a t 
a, —w x ft) x 2 -]“fc x -1—JS? 2 sin cojt 

This fraction being merely the ratio of two lengths, it must be 
independent of the time ; thus the right-hand factor must reduce 
to a constant, i.e. 

sin co a t = c.sin <o x <. 



where c is as yet undetermined. If possible, let c < 1, and put 
t — ji/2co a , then 


1 


c.sin 


jiojx 
2 co 2 


i.e. sin —- > 1, which is impossible. 

2o>, 


Thus c cannot be less than unity; similarly, it can be 
shown that c cannot be greater than unity; hence c = 1, and 
sin a) a t = sin o) x t. For this last result to be true at all times, 
<w 2 = ct) x . Let co = coi = ft) 2 , then the ratio a x /a 2 becomes : 


a x _ k a cos <f> 
a a n/t x o)^ ^~k x ~\~k a 

This same ratio can be evaluated from (13.46) as 
a x _ ( — m,ft) a -|-& 2 ) cos <f> 
d a k a 


. (13.5a) 


. (13.56) 


By equating the two expressions for the ratio, the so-called 
“ frequency equation ” is obtained : 


ft)" 


Lw 2 m 1 J 


k x k a ' 
m x m a 


= 0 


(13.6) 
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As this equation is a quadratic in co 2 there are two natural fre¬ 
quencies of the system, oj L /'2n and a> H / 2n o.p.s. where 


2(0^— *» + *l±*!._ 
m 2 m 1 

2co b * = + 

m i m l 


A _ f / fc, | t 1 +ib 1 \» _ 4^,-1* 
! L\^2 / m^aj 


&2 _|_^1 ^ 


Rearranging the expression in square brackets, 


2cor 2 = 


^2 


A/ 2 | JCi~\~Jc2 


m! / w^ratj 

i±**+cl 

J 


= - + 

m 2 


where (7 > 0 


^i+&a 


= 2*1+*! - (7. 


Hence aj £ 2 < *1+*? . . . (13.7a) 

Similarly it can easily be shown that 


<o H 2 >— 2 . . . (13.76) 



Fig. 17.—Displacement: time curves for the two masses of Fig. 16 : (a) lower 
natural frequency; ( b) higher natural frequency. 


It is convenient here to regard the amplitudes as being essentially 
positive (see note in Section 4); the results (13.7) when sub¬ 
stituted in (13.5) then give : 

cos <f>= 1, i.e. <£ = 0 if co = <o L 1 

cos <f> — — 1, i.e. <f> = n if co = a> H J • \ • ) 

(It has already been shown that cos <f> = ± 1.) 

Two very important results have now been obtained: 


(i) both masses vibrate with the same frequency, and 

(ii) there are two natural frequencies, at the lower of which 
the two masses move in phase (<ft = 0) as in Fig. 17a, whereas 
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at the higher frequency the two masses always move in opposite 
directions = n) as in Fig. 176. 

14. Deformation diagrams and nodes. 

The configuration of the system when vibrating is known 
as the “ deformation shape,” and is illustrated by a “ deforma¬ 
tion diagram ” similar to Fig. 18, in which the amplitudes of the 
motion of the two masses are plotted at right-angles to the axis 
of the system. As the springs are supposed to be composed of 
uniform material uniformly distributed, the deformation shapes 
for the springs are straight lines AB' and B'C', where BB' and 
CC' are the plotted displacements of the two masses and A is the 


A A 



Fig. 18.—Deformation diagrams for the two masses of Fig. 16 : (a) lower natural 
frequency; (6) higher natural frequency. 

suspension point. The deformation diagrams are quite different 
for the two natural frequencies. 

(i) Lower natural frequency {at = a> L ). 

From (13.56), when (o — co L 

«i _ i _ ra a ct) £ 8 j 
kf 

Moreover, a t and a, are positive; hence a x < a, and the deforma¬ 
tion diagram has the form indicated in Fig. 18a. 

(ii) Higher natural frequency {co — a» H ). 

The masses move in opposite directions, and there is a point 
N on BC such that the displacement at N is zero at all times, 
Fig. 186. Such a point is termed a node. 

At both natural frequencies the displacement of the suspension 
point A is zero, since the, support is taken to be rigid and so 
does not move. The point A is therefore a node in both cases. 
The shape of the deformation diagram, specified by the 



48 


FUNDAMENTALS OF VIBRATION STUDY 


— 0 



ji 


c 2 

n Ci 








u 

_ 


Fig. 19.—Three-mass system 
with two natural frequencies. 


position of the nodes, is termed the mode . The two modes of 
vibration of the two-mass system of Fig. 16 are therefore : 

(i) a one-node mode, the node being at A, and 

(ii) a two-node mode, one node being at A and the other 

between the two masses. 

In Chapter I, section 6, it was 
shown that the torsional two-mass sys¬ 
tem of Fig. 4 has only one natural 
frequency ; similarly it can be shown 
that the three-mass system of Fig. 19 
has only two natural frequencies. It 
is left as an exercise for the reader to 
show that the corresponding modes of 
vibration are : 

(i) a one-node mode, and 

(ii) a two-node mode. 

(See Exercise 1 at the end of the chapter.) 

It must be emphasised that the same result would hold for 
a linear system consisting of three masses, adjacent pairs being 
connected by springs, the motion being in a horizontal line on 
the surface of a smooth table ; the torsional case is more prac¬ 
ticable, however, and is therefore to be preferred as an illustra¬ 
tion. It will be seen from the results obtained that the number 
of natural frequencies is, in all cases of 
systems so far examined, the same as the 
number of springs, whether linear or tor¬ 
sional. This is true of any multi-mass 
system in which the masses and springs 
occur alternately (see Chapter IV, section 17). 


22 YMr/rsfr 




P=»P 0 sV\pt ’ 


m, 


15. Undamped vibration absorber. 




jnt2 


The system illustrated in Fig. 20 is fig. 20 .—Undamped 
a simple type of dynamic absorber lay- vibration absorber 
out. The mass m x is forced by the input 1 ^ 

P = P 0 sin pt } and it is desired to reduce the amplitude of 
vibration of m x to zero at a particular value of the forcing fre¬ 
quency p /2tz. The auxiliary spring-mass system k 2 , m 2 effects this 
by absorbing the input energy at its own natural frequency. By 
precisely the same method as used in Section 13 it can be shown 
that the two masses move with the same frequency p/2n as the 
exciting force, and further that they move either in phase or 
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n radians out of phase with each other. Allowing the amplitudes 
to take negative values, the displacements of the two masses 
can therefore be expressed as : 


x x = a x sin pi 
x 2 = a 2 sin pt 


(15.1) 


The equations of motion are : 


m l x 1 +(k 1 +k 2 )x 1 — k^c 2 = P 0 sin pt 
m^c 2 — k 2 x x -\-k^c 2 = 0 


} 


(15.2) 


Substitution of the sinusoidal expressions (15.1) in the equations 
of motion, and elimination of a 2 between the two resulting equa¬ 
tions, lead to an expression for the amplitude a x of motion of 
the mass m 1 : 


= _ Ppfe - rn 2 p 2 ) _ 

1 (k 2 — m 2 p 2 )(k 1 +k 2 — m x p 2 ) — k 2 2 


(15.3) 


If the forcing frequency p/2n is such that m 2 p 2 = k 2} the numera¬ 
tor of the fraction (15.3) is zero but the denominator is not; 
thus at the natural frequency of the auxiliary system, considered 
as a separate system similar to Fig. 1, the amplitude of motion 
of the mass m 1 is reduced to zero. The auxiliary system is thus 
absorbing all the vibratory energy in the system, and is termed 
a dynamic absorber . The name “ dynamic damper ” is fre¬ 
quently used but in error, as the property utilised is that of 
absorption, not of dissipation by way of damping. 

, The arrangement is of great utility in cases where it is required 
to reduce the amplitude of motion of a system without altering 
the suspension arrangements. In practice the damping forces 
must be taken into account; Den Hartog deals with the subject 
of the damped absorber in some detail (see reference 2 in the 
Bibliography). The dynamic absorber in its torsional form is 
of still greater utility. 

Two fundamentally different types of absorber are being 
applied to an ever-increasing number of troublesome torsional 
systems. One type is merely a torsional version of the linear 
absorber, the spring attachment being in some cases a patented 
rubber mounting; this variety of absorber is effective at a 
particular frequency in the same way as is the linear analogue. 
The other type is more interesting and at the same time more 
useful; it consists essentially of a pendulum hinged to a con¬ 
venient part of the system—in the case of internal combustion 
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engines, to a crankweb—and the dynamic conditions of the 
system are such that the pendulum absorbs energy, not at a 
certain constant frequency, but at a frequency which is a constant 
multiple of the rotational speed of the engine. Thus, for example, 
the absorber may be effective for frequencies x engine r.p.m., 
and it is then termed a “ 2\ x ” absorber. By suitable design of 
absorbers two or three harmonics (2J, 3£, and 6 x engine r.p.m., 
for example) may be absorbed simultaneously. A very full dis¬ 
cussion of this important subject is given by Dr. Ker Wilson 
(see reference 3 in the Bibliography). 


16. Coupled vibrations. 

Fig. 21 illustrates a system of the second type mentioned at 
the beginning of this chapter—one in which a single mass has 



lx 


Fig. 21.—One-mass system with two natural frequencies. 


two degrees of freedom. The centre of gravity of the mass is 
constrained to move vertically, and the springs react against 
this motion with restoring forces in the vertical direction ; and 
the mass rotates about its centre of gravity, also against the 
spring restraints. If 0 is the centre of gravity of the mass, the 
two degrees of freedom are therefore : 

(i) vertical motion of 0, and 

(ii) rotation of the mass about O. 

When both linear and angular displacements occur in the 
same problem it is convenient to denote them by capital and 
small letters respectively; in this system let the two displace¬ 
ments be X and y as shown in the diagram. The arrow-diagram 
on the right of the figure shows the positive directions of motion. 
In general both displacements occur at the same time, and for 
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small amplitudes the spring forces produced by the motion are 
proportional to the displacements. Suppose first that the mass 
moves vertically with no rotation, the displacement being X 
downwards. The spring k x exerts a restoring force k x X upwards, 
and this force has a moment k x d x X about G in the clockwise 
direction; similarly the spring k 2 exerts a restoring force k 2 X 
upwards, and this force has a moment — k 2 d t X clockwise about 
0 . Thus the total restoring force on 0 is (k x +k 2 )X upwards, 
and the total moment about 0 in the positive sense is 
(k 1 d 1 — k$ 2 )X. This moment tends to produce a rotation round 
0, and the two types of motion are therefore coupled . It is to 
be noted that the displacements are assumed small, so that the 
moment-arms d x and d 2 do not change appreciably during the 
motion. 

Now suppose the mass to rotate through a small angle y 
without linear motion of the centre of gravity G ; the rotation 
being about G , the displacements at the points of attachment 
of the springs are — d x y and d$ downwards, as y is positive 
when clockwise. The forces exerted by the springs are there¬ 
fore k x d x y and — k 2 d$ downwards, and these forces have moments 
— k^^ and — k 2 d 2 2 y clockwise about G. The total spring 
force on G is thus (k x d x — k^d^y downwards, and the total 
moment about G is — {k x d x 2 +k 2 d 2 2 )y clockwise. 

If now the mass is supposed to move generally—i.e. to have 
both the displacements X and y —the forces and moments induced 
are the sum of those induced by the two motions occurring 
separately* Let the moment of inertia of the mass about G 
be J, with respect to the rotation y, then the equations of 
motion are : 

wJ5l= — (k x +k 2 )X+(k x d x k 2 d 2 )y 1 na 

Jy = (W - kA)X - (k x d*+W)y ] ‘ { } 

In the usual manner (see Section 13) it can be shown that the 
displacements can be expressed as : 

X = sin cot 
y = a.sin (o)t+f>) 

where <f> = 0 or n. 

It is convenient here to allow the constant a to take negative 
values, and to write accordingly : 

X = A. sin cot 
y = a.sinarf 



( 16 . 2 ) 
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where A is positive but a may be either positive or negative. 
Substitution of these values in (16.1) leads to the equations : 


[meo 2 — (k 1 +k 2 )]A+[k 1 d 1 — k^d 2 ]a = 0 
[Jco 2 — (k 1 d 1 2 +k 2 d 2 2 )]a + [k x d x — k 2 d^\A = 0 


which can be written as : 

(mco 2 — o xx )A-\-c Xy a — Ol 

C yxA -{-(Jco 2 “ Cyy)® = 6J 

where : 

c xx = k 1 -\-k 2 1 

6 yy = k x d x 2 -\-k 2 d 2 2 > 

C X y = Cy X = ^ 1^1 k 2 d 2 J 


(16.3) 

(16.4) 


The substitutions indicated in (16.4) are very generally useful, 
for — c xx is the force parallel to the X-axis produced by a unit 
displacement along that axis, — c Xy is the force parallel to the 
X-axis produced by a unit angular displacement about the 
y-axis, etc. The equality of the two constants c Xy and c yX is 
not merely coincidental, but is an example of a well-known 
reciprocation theorem. 

By obtaining the ratio A /a in each of the equations (16.3) 
the “ frequency equation ” is obtained : 


co 4 


'CxX_i_ c yy\ 2 _^ C XxCyy ~ c Xy c yX = q 

m J) mJ 


or, as c Xy — c yXi 

co* _ (^+ e jg>\(oi+ c x&vv ~ C *L = o . (16.5) 

\m J J mJ 

(16.5) is very significant when expressed in terms of the “ un¬ 
coupled ” frequencies of the system. The constants c Xy = c yX 
are a measure of the degree of coupling between the two modes 
of vibration, i.e. between X and y. If these equal constants are 
zero, which occurs when 

kede = kjd t .... (16.6) 

the two degrees of freedom are uncoupled and each can take 
place without involving the other; for a displacement X then 
induces zero moment about 0, and an angular displacement y 
induces zero force parallel to the X-axis. (16.3) then become: 

(mco 2 — c X x)A = 0 

(Jco* — c w )a = 0 
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and if the corresponding values of co are denoted by co x , cd v , 


(U X Z 

COj, 2 



(16.7) 


From (16.7) the “ uncoupled ” frequencies (o x /2n and co y /2n can 
be calculated. (16.5) can be written as: 

co* — ((o x 2 +a) y 2 )(o 2 +(a) x 2 co 2 — B) = 0 . (16.8) 

where B = > 0. 

mJ 


Denoting the roots of (16-8) by oj l 2 and co H 2 , where m L 2 is the 
smaller root, the elementary theory of equations shows that 


“>l 2 +<°h 2 — co x 2 +co y 2 


o>l<»h 2 = co x 2 m y 2 


B <a>, 


CO-, 


(a) 

(b) 


(16.9) 


By using the identity (m+n) 2 = (m — n) 2 -\-4mn it is easy to 
deduce from (16-9) that 


(o H 2 — u> L 2 > to x 2 — o) y 2 numerically 


(16.10) 


(16.9a) shows that the mean of the two coupled values of on 2 is 
equal to the mean of the two uncoupled values, while (16.10) 


| UNCOUPLED 

1 

1 

_ 

COUPLED 1 

1 

■* ' "1 

1 

1 

1 
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I I I 


or 


as 


Fia. 22.—Frequency spectrum for Fig. 21. 


shows that the difference between the coupled values is greater 
than the difference between the uncoupled values. The lower 
coupled value is therefore lower than the lower uncoupled value, 
and the higher coupled value is higher than the higher uncoupled 
value, as shown in Fig. 22. Such a diagram is termed loosely 
a “ frequency spectrum.” 

If the mass in Fig. 21 is free to move in a horizontal direction 
against spring restraints, the system then has three degrees of 
freedom and the equation corresponding to (16.8) is a cubic in 
co 2 . The general case of a mass free to move in space under 
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spring restraints, as shown in Fig. 23, has six degrees of free¬ 
dom, namely, three linear and three angular displacements. 
The resulting “ frequency equation ” is of the sixth degree 
in co 2 , but symmetry in the system can reduce the equations of 
motion to separate groups by rendering some of the coupling 
coefficients zero. This occurs in the orthodox type of in-line 
aero-engine mounting ; in such a system the vertical motion, 
for example, is uncoupled from rolling motion (i.e. rotation 
about a fore-and-aft axis) while usually remaining coupled to 
pitching motion (i.e. rotation about an axis parallel to the wings 
of the aircraft).* The more symmetrical the system is, the 
simpler the solution becomes ; in the case of a solid of revolution 
such as an airscrew spinner, with spring mountings symmetrically 
disposed about the axis of revolution, the general solution reduces 



Fig. 23.—One-mass 



to two quadratic equations similar to (16.8) and two “ uncoupled ” 
equations similar to (16.7). 


Summary. 

All parts of a system with two degrees of freedom vibrate at 
the same frequency. There are two natural frequencies. In the 
case of the “ chain ” type of system illustrated in Figs. 16 and 19, 
there are respectively one and two nodes in the configuration at 
the two natural frequencies. In such a system, a final mass 
attached by a spring can eliminate forced vibration at its point 
of attachment by “ forcing a node ” there ; this fact is the basis 
of the dynamic absorber. 

In the type of system illustrated in Fig. 21, there js coupling 
between the two degrees of freedom possessed by the mass 
(horizontal motion being supposed to be prevented). It is found 

* See Appendix III, section 41. 
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that the mean of the two values of co 2 corresponding to the two 
natural frequencies is the same as the mean of the two values 
corresponding to the so-called “ uncoupled ” modes, in each of 
which one type of motion is supposed to exist in the absence of 
the other. 


EXERCISES III 

1. By the methods of Sections 6 and 13 show that there are only 
two natural frequencies of the three-mass torsional system of Fig. 19, 
and that the corresponding modes are : 

(i) a one-node mode, and 

(ii) a two-node mode. 

2. Derive an expression for the amplitude of motion of the lower 
mass (m 2 ) in Fig. 20, and draw curves showing how the amplitudes of 
motion of both masses vary with the forcing frequency p/2n. 

(It will be found that these amplitudes have infinite values when the 
forcing frequency equals the natural frequencies of the system. If small 
damping forces are included in the system the infinite amplitudes are 
reduced to the practical resonance “ peaks.”) 

3. —Find the condition that one of the coupled frequencies of the 
system shown in Fig. 24 should be zero, the centre of gravity of the mass 
being constrained to move in the direction X. 




Fig. 24.—Ex. Ill, 3. 


(Obtain a frequency equation analogous to (16.5) and find the con¬ 
dition that a) a is a root, i.e. that the term independent of co 2 is zero.) 

4. Derive the frequency equation for the system of Fig. 23 in terms 
of appropriate coupling-coefficients similar to (16.4). 

(See Appendix III, section 41.) 

5. Determine the force which mustlbe applied to the top of the spring 
k 2 in Fig. 16, page 44, in order to maintain unit amplitude of vibration 
thereat, ignoring the upper part of the system ; similarly, determine the 
force which must be applied to the upper mass (tn x ) to maintain unit 
amplitude thereat, ignoring the lower part of the system. Hence determine 
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the natural frequencies, as those frequencies at which zero total force is 
required to maintain a vibration (see page 74). 

6 . Vibratory torques of amplitudes T 1 and T 2 , and of the same 
frequency, act upon the inertias in the system of Fig. 4, page 14. Show 
that if the torques are in-phase and T 1 /T 2 = J 1 /J 2 , there is no twisting 
of the shaft; and that if the torques are anti-phased and T 2 = T v the 
vibration is in the mode corresponding to the natural frequency. Show 
also how the general case can be analysed by regarding the applied torques 
as resulting from the superposition of two sets of torques fulfilling these 
special conditions. 



CHAPTER IV 


MANY DEGREES OF FREEDOM 
Effective Inertia and Dynamic Stiffness 

Introductory. 

By forming a chain of alternate springs and masses, vibrating 
systems with any number of degrees of freedom can be constructed. 
The linear form of such multi-mass systems is of no particular 
interest, but the torsional form, in which a number of inertias 
are connected by light shafts, is of great importance in the study 
of the stresses in internal combustion engine crankshafts, and 
more generally in power-plant installations of all kinds. 

Of recent years there has been a great development of certain 
special methods of torsional vibration analysis. The problem is 
to determine the natural frequencies of a torsional system, given 
the values of the inertias and shaft stiffnesses. Once these 
frequencies are known, and provided there is sufficient information 
on the magnitude and nature of the forces causing vibrations, 
it is possible to analyse the system with a view to the avoidance 
of the dangerous dynamic magnifications due to resonance. For 
practical purposes it is very convenient to evolve a method of 
analysis in which various parts of the total system are first 
considered separately, with means for the immediate determina¬ 
tion of the natural frequencies of the complete system. Many 
different types of installation can thereby be treated on a unit 
basis. An aircraft power-plant can be split conveniently into 
two parts, the engine and the propeller ; an industrial compressor 
plant into two parts, the engine and the compressor; and so 
on. The greatest utility of such a method is seen where there is 
likely to be a considerable amount of interchanging of the units, 
different propellers (for example) being used with the same 
engine. 

Such a method can be based on the now familiar concept of 
effective inertia , and a very similar method proceeds by considering 
dynamic stiffness. The two methods are very closely related, and 
a choice between them depends largely upon the precise nature 
of the particular problem under consideration. They have been 
described in detail elsewhere by the author (reference 15 in the 
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Bibliography, page 151), and the application of the dynamic 
stiffness concept to the analysis of damped vibrations is described 
in reference 4. 

The system being regarded as composed of two sub-systems 
joined at a definite point in the shafting, for each part a character¬ 
istic quantity is determined either by calculation or by experi¬ 
ment. This quantity, which is a function of the frequency of 
vibration, is the amplitude of alternating torque which must be 
applied at the junction to maintain a vibration either of unit 
acceleration amplitude or of unit displacement amplitude. In 
the former case the characteristic quantity is then the effective 
inertia of the sub-system at the junction point: i.e. the inertia 
of a rigid flywheel which would vibrate with the same amplitude 
as does the actual sub-system under equal forcing torques. In 
the latter case the characteristic quantity is the dynamic stiffness : 
i.e. the stiffness of a shaft, one end of which is fixed, which is 
imagined to be acted upon by an alternating torque at the free 
end, and which vibrates with the same amplitude as does the 
actual sub-system under equal forcing torques. Whether effective 
inertia or dynamic stiffness is chosen, the procedure is to plot 
the variation of this quantity against frequency for each sub¬ 
system, reverse the sign of the graph for one part and superimpose 
it on the graph for the other part; the frequencies determined 
by the points of intersection of the two curves are the natural 
frequencies of the complete system, as shown on page 69. 

It is found that for practical purposes, in lightly-damped 
systems the damping forces can be neglected without causing 
any serious error except in the accuracy of amplitude determina¬ 
tions in the immediate neighbourhood of resonances, and the 
working is greatly simplified in this way. 

Several methods are given in the text for the determination 
of the deformation shape of the system at any frequency, from 
a knowledge of which the stresses in various parts of the shafting 
can be calculated. 

17. Multi-mass systems. 

The purely analytical methods employed in Chapters I-III 
prove to be cumbersome when applied to systems having more than 
two degrees of freedom. Particularly is this so in the calculation 
of resonant frequencies by differentiation of amplitude-against- 
frequency curves, but fortunately it is found that for lightly- 
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damped systems the damping forces can be disregarded and the 
natural frequencies of the resulting system calculated instead of 
the resonances. This useful fact has been proved for systems 
with one degree of freedom in Chapter II, 


section 11, and will be proved for multi¬ 
mass systems in Section 18 ; and when 
dealing with lightly-damped systems at¬ 
tention will be confined to determining 
the natural frequencies and deformation 
shapes. The analytical method can be 
used to show that a torsional system as 
that of Fig. 25a has one less natural fre¬ 
quency than the number of inertias ; if 


vht Jn-l ^3 Ja Ji 




(a) (b) 

Fig. 25.—Multi-mass systems: (a) n—1 natural frequencies; (b) n natural 

frequencies. 

the last inertia J n is removed and the system clamped at this 
point, thus replacing the finite inertia by the effectively infinite 
inertia of the earth, to which the clamp is rigidly attached, there 
will still be n — 1 natural frequencies. These results are easily 
derived, as follows : ■ 

Let the angular displacements of the inertias from the 
static position be 6 U 0 8 , 0 3 . . . 6 n . There are n equations of 
motion : 

J = — Ci(0i — 0 8 ) 

JA =- 0,(6, - 0a)+<^(0, - 0.) 

etc. etc. 

= — C n _ 1 (d n _ 1 — 0„)+<7 n _ 2 (0»_2 — 1) 

•AAi — +^n-l(0n—1 — ®n) 

By extension of the method of Section 6 it can easily be shown 
that all the inertias have sinusoidal motion at the same frequency^ 
and that all the phase-angles differ from each other by zero or 




FUNDAMENTALS OF VIBRATION STUDY 


a multiple of n radians, and that the displacements can therefore 
be put in the form 


Q x = A 1 .sin cot' 
0 a = -4 2 • sin cot 
etc. etc. 

0 n = A n . sin orf 


(17.2) 


where the constants A are allowed to take negative values (see 
note on the sign conventions, Chapter I, section 4). The 
equations of motion can then be written as : 

(C t - J 1 co 2 )A 1 = C X A 2 

(Ci~\-C 2 J 2 co z )A 2 = CA.+CA, ^ gv 

etc. etc. I • v • / 

Elimination of the amplitude ratios A 2 /A l9 A 3 /A l9 etc., leads to 
an equation of degree n in co 2 , of which co 2 — 0 is a solution 
corresponding to A x = A 2 = . . . = A n . As in Section 6 this 
solution refers to steady rotation, and there are thus n — 1 
natural frequencies of vibration. It has not been shown that 
these natural frequencies are all distinct, but that they are so 
will be evident in the course of the present chapter ; a proof of 
the validity of the assumption in the case of “ straight chain ” 
systems comprising alternate masses and springs, or alternate 
inertias and shafts, is given in the Appendix of reference 15. An 
important point to be noted is that the number of natural 
frequencies is equal to the number of springs, in systems similar 
to those of Fig. 25, wherein the inertias and springs occur 
alternately; whether the system is a linear or a torsional one 
is, of course, immaterial. 


18. Effective inertia. 

The multi-mass systems that occur most frequently in practice 
are torsional systems, and the discussion which follows is therefore 
concerned primarily with torsional vibration; the equivalent 
results for linear systems can be obtained by means of the linear- 
torsional analogues listed in Table I (Chapter I, section 5). 

When a torsional system such as that of Fig. 25a is vibrating 
at a natural frequency in the absence of damping forces, the 
motion is maintained by the energy stored in the system initially. 
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If the displacement at a point P in the system is 

0 = A . sin cot 

the acceleration is & = — Aco 2 sin cot, 

and there is no applied torque required to maintain this accelera- 



Fig. 26.—Partial inertias. 


tion. The effective inertia of the system towards vibrational 

motion at a natural frequency is thus zero, for 

. ,. torque 

inertia =-—-. 

angular acceleration 

If it is possible to determine the effective inertia of the system 
as a function of frequency, the natural frequencies can be found 
by investigating for what frequencies the effective inertia is zero. 

Consider the two-mass system shown in Fig. 26. The natural 
frequency is 

F = (see Seoti ™ 6) ’ 

Suppose that the inertia J 2 is divided into two parts, A, and p t , 
the partial inertia A* being such that the natural frequency of 
the sub-system J x — {C) — A* is xo/2n. Thus 

-JW 

GJi . A 

i e< Xt ~ - G (a) i . . ( 18 . 1 ) 

and — J t — A* (6), 
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jn i is then the effective inertia J e of the system at the frequency 
co/2n, for the remainder of the system is vibrating at its natural 
frequency and thus has zero effective inertia. The effective 
inertia J e of the whole system is thus given by : 


Jo — e/« — A 2 — J 2 


CJ X 


Ct) 2 J ! — C 


. (18.2) 


Putting J e = 0 naturally leads back to the formula (6.4). The 
method is, however, capable of immediate extension to multi-mass 
systems. In the system of Fig. 25a let 


J a = ^a+/^2 

J 3 = ^ 3+^3 > 

etc. etc. [ 
^n— K +/V 


(18.3) 


As before, let the sub-system J x — (C^) — A a have a natural 
frequency (o/2n ; then as in (18.1), 


— J 2 


G 1 J 1 

co 2 J 1 - G * 


(18.4a) 


Now determine A 3 so that the sub-system jn 2 — (C 2 ) — A 3 also has 
a natural frequency (jo/2n. Thus 


l* 3 — J 3 ^3 — «^3 — 


C 2 ^2 


(18.46) 


a) 2 // a — C 2 

In this manner the whole system is split up into sub-systems, for 
each of which co/2n is the natural frequency. In the next stage 
after (18.46), 

. . . (18.4c) 

and finally, for the last inertia J n , 


f^n — J n 


(18.4d) 


«>X-1 - Cn-\ ' 

The whole system can now be considered as consisting of two 
parts: 


(i) the system J 1 — J, — J a . . . J n _ t — X n 

(ii) the remaining inertia fi n . 

Part (i) is a system vibrating at a natural frequency, and has 
therefore zero effective inertia ; thus the effective inertia of the 
whole system at the point J n is fi n . 
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Apart from providing a method of determining the natural 
frequencies of multi-mass systems, this conception of partial- 
inertias and sub-systems enables two very important general 
conclusions to be drawn : 

(i) All the masses move so that their displacements can be 

expressed as in (17.2), for it has already been established 
(Section 6 in Chapter I) that the two partial inertias 
comprising each sub-system move in this manner, and 
the two partial inertias comprising each inertia J l9 J t , 
J 3 , etc., are rigidly connected. Thus any pair of 
adjacent inertias move either in phase or n radians out 
of phase with each other. 

(ii) The resonant frequencies of a lightly-damped multi-mass 

system are very nearly equal to the natural frequencies 




^3 




Cl 


r 

1-5 


LBS. INS. /MICRO- RADI AN 


2-5 3 2 LBS. INS. SECS. 

Fig. 27.—Four-mass torsional system. 


calculated on the assumption of zero damping; for 
this result has also been proved for a single-mass, 
single-degree-of-freedom system, and therefore by im¬ 
plication for a two-mass, single-degree-of-freedom 
system such as that of Fig. 4 with light damping 
forces included ; but the sub-systems are of the same 
type as that of Fig. 4, hence the general result. 

(A more rigorous demonstration of this second result 
is to be found in reference 4 in the Bibliography at the 
back of the book.) 

As an example of the application of the method, consider 
the four-mass system of Fig. 27. The stiffness values are given 
in the practical unit: lbs.ins. per microradian (i.e. per millionth 
of a radian) and must be multiplied by a million before being 
used in the formulae (18.4). 

At the frequency 13,500 C.P.M., i.e. 225 C.P.S., 

co = 2tc X 225 = 1,414 and co 2 = 2 x 10 6 . 
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At this frequency, 


Ctftt _ 2-220 


and 


Hi — J 3 
J e = Hi — J i 


(0 2 fl t — C g 

C 3H 3 


C0 3 fl a 


= 1 040. 


At the frequency 6,750 C.P.M., i.e. 112*5 C.P.S., 

(co = 2n X 112*5 = 707 and eo s = 0*5 X 10«) 
Hi = i °°> Hs = J e ~ H* ~ i c0 * 



At the frequency 9,550 C.P.M., i.e. 159 C.P.S., 

(<w = 2 n X 159 = 1,000 and <y 2 = 1*0 x 10 6 ) 

Hi = !, Ha = 1 - 6 , J e = H* ~ 0 . 

Thus 9,550 C.P.M. is a natural frequency of the system. The 
effective inertias at J 4 of the system are calculated in a similar 
manner for various frequencies so as to determine the other 
natural frequencies. Values are given in Table III, and the 
curve of effective inertia plotted against frequency in Fig. 28. 
The three natural frequencies are seen to be: 

4,100, 8,620, and 9,550 C.p.m. 

The curve in Fig. 28 is typical of effective-inertia curves. 
The zeros correspond to natural frequencies of the system of 
Fig. 27, and the discontinuities correspond to the natural fre¬ 
quencies when the inertia J 4 is clamped so that its displacement 
and hence its acceleration are necessarily zero, for the ratio 
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Table III 


Frequency 

(O.P.M.) 

o>* 

-r 10 « 

Je 

0 

0 

9 

955 

001 

10-6 

1,910 

004 

30-1 

3,020 

01 

- 6-83 

4,100 

0*1838 

0 

5,230 

0*3 

1*97 

6,040 

0*4 

4-26 

6,750 

0*5 

oo 

7,400 

0*6 

- 4-89 

8,000 

0*7 

- 0-84 

8,620 

0-8162 

0 

9,160 

0-9216 

0-98 

9,360 

0-9604 

- 1-03 

9,550 

1-0 

0 

10,500 

. 1*21 

0-58 

etc. 

etc. 



torque/acceleration is then infinite, and the effective inertia is 
equal to this ratio. Between adjacent zeros there must be one 
discontinuity, and between each pair of adjacent discontinuities 
there must be one zero, so that zeros and discontinuities occur 
alternately as the frequency is increased ; for the discontinuities 
correspond to the natural frequencies when an additional con¬ 
straint is applied (the clamping of J t ) and it is a general dynamic 
theorem that the natural frequencies of such a system must lie 
between the natural frequencies of the original system. (See 
E. T. Whittaker, Analytical Dynamics, Cambridge, 3rd Edition, 
pp. 191, 192 ; appendix of reference 15 ; and Biot, reference 5, 
from whom the reference to Whittaker is taken.) 

19. Swinging forms. 

The deformation curve of a torsional system undergoing 
vibration is called the swinging form, and it is frequently very 
important to know the shape of this curve, i.e. what are the 
relative amplitudes of motion of the various inertias. If the 
swinging form of a system at a particular frequency is known, 
then from experimental observations of the amplitude at one 
point in the system the amplitudes at all other points can be 
calculated. 


v 
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In the two-mass system of Fig. 26, let the displacements of 
the two inertias J x and J t be 

• = A x sin out 

and 0, — A t sin out. 

The velocities are 6 X — A x ou cos out 
and 6 t = A cos out. 


At a natural frequency of the system there is no applied force 
required to maintain the motion in the absence of damping 
forces, so that the rate of change of the total angular momentum 
must be zero, by Newton’s second Law of Motion 
(see Appendix II, section 37). 

Thus J x S x -\-J t St = 0, i.e. 

J x A x -\-J t A t = 0, 


or 


A, _ J, 
A x J t 


(19.1) 


Now, the partial inertia A 2 has been chosen so that the sub¬ 
system J x — (C) — A* has (o/2n as its natural frequency (18.1a), 
and therefore (19.1) holds true at any frequency ou/2n if A, is 
substituted for J t . Similarly, in the system of Fig. 26a, the 
following amplitude ratios are obtained : 


At _ Jx 

A x A, 

At = _ 

A t A 3 


(19.2) 


etc. etc. 


•^fi _ A*n— 1 

A n — i A n 

An arbitrary value of unity is usually assigned to the ampli¬ 
tude A x , and the equations (19.2) then give the corresponding 
values of the amplitudes of motion of the other inertias. As 
an example, Table IV gives the amplitudes of the four inertias 
of the system of Fig. 27 at the three natural frequencies. 


Table IV 


Frequency 

(C.P.M.) 

A x 

At 

At 

At 

4,100 

1 

0-6325 

- 0-7999 

- 1-2648 

S,620 

1 

- 0-6325 

- 0*7999 

1-2648 

0,550 

1 

- 1 

1 

- 1 
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The corresponding swinging forms are illustrated in Fig. 29 ; the 
three modes are: 

(i) one node between J, and <7,; 

(ii) two nodes, between J t and J t and between J, and J«; 
and (iii) three nodes, one in each spring. 




Fig. 29.—Swinging forms- for Fig. 27 (at natural frequencies). 

The method apparently breaks down when any of the partial 
inertias are infinite, as for example at the frequency 6,750 c.p.m. 
At this frequency 

A, = oo fi t = — oo 
A* = 1 = 1-5 

Thus ^* = —^1 = 0, and A t = 0 

Aj. A, * 

A u 

and _i = — c? = oo, whence A , = oo x 0 

A 2 A a 

and is indeterminate. The amplitude can however be determined 
in such cases by a consideration of the torques in the shafts to 
either side; in this example let T l and T t be the torques in the 
shafts CV-and C t . These torques are proportional to the twists 
in their respective shafts, and A t = 0. Hence 

— Ci(A t — A i) = — CiAi 

T% = C t (A» — A t ) = C t A t . 


and 
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Furthermore, as the amplitude A* is zero, the acceleration of J, 
is zero, and hence the resultant torque on is likewise zero ; i.e. 

Tt — T\ = C S A j-j-UjAi = 0, 


and hence 



Thus if Ax = 1, A, = — 2. In a similar manner the swinging 
form can be obtained at any frequency whereat there are nodes 
at any of the inertias. 

If it is desired to know merely how many nodes occur in the 
swinging form at a particular frequency, it is only necessary to 
count the number of discontinuities in the effective-inertia curve 
between that frequency and zero frequency. At these discon¬ 
tinuities the effective inertia is infinite and thus A t = 0. As the 
frequency is increased through such a value, an additional node 
is introduced into the swinging form. Thus for the system of 
Fig. 27, the effective inertia curve for which is shown in Fig. 28, 
if F is the frequency, 


for 0 < F < 2,100 c.p.m., no nodes ; 

2,100 < F < 6,750 c.p.m. , one node ; 
6,750 < F < 9,320 c.p.m., two nodes ; 
9,320 c.p.m. < F, three nodes. 


20. Rigidly-coupled systems. 

The system illustrated in Fig. 30 consists of two definite 
parts A and B, rigidly coupled at P ; part A is a system in 



Fig. 30.—Rigidly-coupled systems. 
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which all the inertia and torsional stiffness values are known, 
and part B is a system in which some or all of these values are 
for some reason unknown. The effective inertia of system B at 
the coupling-point P can be determined experimentally by apply¬ 
ing at this point a sinusoidal torque 7.sin cot and measuring the 
amplitude A of the resulting motion at P ; for if the displace¬ 
ment is 6 = A .sin cot , the acceleration is S = — Aco 2 sin cot, and 
the torque required to produce this acceleration is 

Jj) = — J^Aco 2 sin cot = y.sin cot. 


The effective inertia J e is thus given by : 

Jt Aco ■' 


( 20 . 1 ) 


and in this formula the quantities T and A are determined 
experimentally for various values of the frequency co/2n. It 
is therefore possible to obtain experimentally a curve of J e 
plotted against frequency. 

Let the effective inertias of the parts A and B be J A and 
J B at P. Part A consists of a system vibrating at a natural 
frequency together with a residual inertia J A at P; part B 
consists of a system vibrating at a natural frequency together 
with a residual inertia J B at P. The total effective inertia of 
the whole system at P is therefore J A +J B , and the natural 
frequencies of the complete system are determined by the 
condition 

J A -\-J B = Q • ... ( 20 . 2 ) 

Suppose that part A is the system of Pig. 27, P being at J t in 
this diagram, so that Pig. 28 is the curve of effective inertia 
J e at P; and suppose further that the effective inertia curve 
of part B obtained experimentally is that shown in Pig. 31a. 
The function J A +J B can be plotted by adding together the 
corresponding ordinates in Pigs. 28 and 31a, but the frequencies 
at which J A +J B = 0 are more conveniently determined by 
plotting J A and — J B on the same axes, as in Fig. 316. The 
two curves intersect where J A — — i*©* Jj+Jb — 0. The 

natural frequencies of the complete system are thus: 2,960, 
6,000, 7,300, 9,200 and 10,060 c»p.m. As the two parts of the 
system are rigidly connected at P there are five spring con¬ 
nections and hence five natural frequencies. The node distribu¬ 
tions at these frequencies are listed in Table V. 
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Je 



Fig. 31.—Determination of natural frequencies of Fig. 30 (J A from Fig. 28). 


This method of calculating the natural frequencies of coupled 
systems is of great use even when all the dynamic constants 
(inertias and stiffnesses) are known, particularly when it is 


Table V 


Frequency 

(O.P.M.) 

Nodes in 
part A 

Nodes in 
part B 

Total 

2,960 

1 

0 

1 

6,000 

1 

1 

2 

7,300 

2 

1 

3 

9,200 

2 

2 

4 

10,600 

3 

2 

6 


desired to find the effect of making alterations to one part only 
of the system, as it is unnecessary by this method to recalculate 
over the whole system; in such cases it is usual to split the 
system as near as possible to the point where the alterations are 
to be made. 
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21. Frequency (Torque-Summation) Tables. 

The treatment, given in Sections 18 and 19, of the concept 
of effective inertia is a convenient one for the purpose of intro¬ 
ducing the method, and also for enabling the two important 
general conclusions (Section 18) to be deduced. By far the most 
convenient practical method of calculating the effective inertia 
of a system, however, is afforded by the Frequency Tabulation, 
or Torque-Summation Table. By this method the effective 
inertia, swinging form, and variation in vibration-torque through 
the system are obtained by means of a systematic tabulation. 

Referring to Fig. 25a, the torque T x sin cot in the shaft G x is 
proportional to the difference between the displacements of the 
inertias J t and J x . Thus 

T x sin cot = C x (A t — A x ) sin cot 
or T x = G x (A t - A x ). 

This torque produces the acceleration S x = — A x co 2 sin cot, hence 
T x = — J x co 2 A x = G x (A t - A x ). 

Let the difference between A x and A % be d x A, so that 

A % = A x — d x A 

then d x A = J -g^A x .... (21.1) 

G x 

Similarly the torque T, sin cot in the shaft G t is given by: 

T t = C t (A t — .A,) 

and the acceleration of the inertia J t is effected by the resultant 
torque acting on </„ i.e. 

— J t co 2 A t — T x — T x 

hence T t = — J x co 2 A x — J t co 2 A 2 = C t {A 3 — A t ). 

Let Ax — Ax — &%A, 

then ),A . , 21 . 2 ) 

k a 

Proceeding in this manner, the amphtude A n of the last inertia 
J n is calculated in terms of A t ; similarly the total torque which 
must be applied to J n to maintain the motion is obtained as the 
sum 

— {JiCo^Ai~{-J^co^Af-\-J 9 (d^A f-f* . . . -f 
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Denoting this torque, by — ST n , the effective inertia J e of the 
system at J n is given by putting — ST n for T in (20.1), i.e. 



(21.3) 


As the amplitude is a constant multiplier throughout the 
calculations it is usually given unit value for convenience. 

As an example the frequency tables for the system of Fig. 27 
at the frequencies 4,100 and 7,990 c.p.m. (co 2 = 0-1838 x 10® and 
0*7 x 10® respectively) are given in Tables VI and VII. 


Table VI 

Frequency table for Fig. 27 at 4,100 C.p.m. 


J 

J CO* 

4- 10® 

A 

J<o*A 
-4- 10® 

(o l A 

-r 10® 

C 

+ 10* 

6A 

2 

0-3675 

1 

0-3675 

0-3675 

1 

0-3675 

3 

0-5513 

0-6325 

0-3487 

0-7162 

0-5 

1-4324 

2 5 

0-4594 

- 0-7999 

- 0-3675 

0-3487 

0-75 

0-4649 

1-5 

0-2757 

- 1-2648 

- 0-3487 

0 

— 

— 


ST 

J e (o 2 = -j-2 = 0, J t — 0 (natural frequency). 
A m 


Table VII 

Frequency table for Fig. 27 at 7,990 c.p.m. 


J 

Jo»* 

* 10 * 

A 

Jm*A 
* 10 * 

XJa>*A 

- 7 - 10 * 

C 

-r 10 ® 

dA 

2 

1-4 

1 

1-4 

1-4 

1 

1-4 

3 

2-1 

- 0-4 

- 0-84 

0-56 

0-5 

M 2 

2-5 

1-75 

- 1-52 

- 2-66 

- 2-10 

0-75 

- 2-8 

1*5 

1-05 

1*28 

1-344 

- 0-756 





- 0-756 
1-28 


= - 0-591, J e = 


0-845. 


22. Dynamic stiffness. 

A quantity having a close relation to the effective inertia of 
a system is the dynamic stiffness. The torsional stiffness of a 
shaft is defined as the ratio of the torque transmitted by the 
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shaft to the twist induced by that torque. If one end of the 
shaft is clamped as in Fig. 32 and a torque T is applied at the 
free end, the angular displacement at the free end is T/C, where 
C is the torsional stiffness. If this displacement is A, then 

0=2 ■ ■ ■ ■ ( 22 . 1 ) 

Now consider a torque T . sin cot to be applied at a point P in 
a vibrating system, and let the resulting displacement at P be 
6 = A.sin cot. The ratio torque/displacement is defined as the 
“ dynamic stiffness ” of the system at the point P, and is denoted 
by Z. Thus 

Z = 2 ' ' • • (22,2) 

The name “ mechanical impedance ” is sometimes given to this 

$ 

T 

ANGULAR DISPLACEMENT A* T /c 

Fig. 32.—Static stiflhess of shaft. 

quantity, by a loose analogy with electrical quantities. Electrical 
impedance is the ratio of e.m.f. to current, and the analogous 
mechanical ratio is that of force (or torque) to velocity, not to 
displacement; while it is quite justifiable to choose any name 
for the quantity defined, it is perhaps unfortunate that the word 
“ impedance ” should be used to denote two ratios of different 
types. Biot (reference 5 in the Bibliography) uses the term 
“ dynamic modulus,” but this choice again is unfortunate as the 
word “ modulus ” has certain other specific meanings in mathe¬ 
matics, and it is indeed necessary to refer to the “ modulus of 
the dynamic modulus ” in the case of damped systems—see 
(22.8) and reference 4 in the Bibliography. The quantity Z 
defined by (22.2) is, however, the precise dynamic analogue to 
static stiffness, and the name “ dynamic stiffness ” has been 
adopted in this book for this reason. 

In engineering practice, torques are normally measured in 
lb.in. units; the basic unit for dynamic stiffness Z is therefore 
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lb.in./radian. As the radian is a large angle, however, it is 
usually more convenient to express the stiffnesses in lb.in./mioro- 
radian. 

Comparison of the equations (20.1) and (22.2) shows that the 
dynamic stiffness Z and the effective inertia J e are related by 
the formula: 

Z = - J e o> 2 . . . (22.3) 

and a curve of dynamic stiffness can easily be derived from a 
curve of effective inertia, or vice versa, by multiplying or dividing 
all the ordinates by the corresponding values of w 2 and changing 
the sign. The condition for the natural frequencies of two systems 
rigidly coupled together is therefore 

Zi+Z t = 0 .... (22.4) 

where Z x and Z t are the dynamic stiffnesses of the two systems 
1 and 2 at the coupling point. (Note : the distinguishing sub¬ 
scripts have deliberately been changed from letters to numbers, 
as both notations are useful on different occasions.) 

Both the effective inertia and the dynamic stiffness are 
measures of the response of a system to an applied sinusoidal 
torque ; the dynamic stiffness is frequently to be preferred when 
it is desired to know the response at one point in a system in 
terms of the response at another point, and more particularly 
when two systems are coupled through flexible shafts or through 
gearing. For these purposes three formulae relating to transfer¬ 
ence of the reference point are very useful. 


(i) Transference through an inertia. 

In Fig. 33a let the dynamic stiffnesses of the systems to the 
right of the points A and B immediately each side of the inertia 
J be Z A and Z B ; the corresponding torques and angular displace¬ 
ments are: 

^) 8 in^ and 

As the inertia J is rigid, A a = A B = A (say) and the difference 
between the torques T B and T A maintains the acceleration 
6 = — Aco 2 sin cot of the inertia. Thus 


j sin cot. 




t b = t a 

and 

a b A a 

whence by (22.2) 

%B = %A 


Jco 2 A 
Jco 2 

— Jco 2 . . . (22.6a) 
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(ii) Transference through a shaft. 

Using the same notation for the conditions of the systems to 
the right of each end of the shaft in Fig. 336, T A = T B = T (say), 
as the shaft is supposed to be light. Furthermore, the amplitudes 
are connected by the relation: 


C{A b — A a ) = T 
T 

Hence A B = A a -\-— 



Fig. 33.—Transference of reference point for dynamic stiffness: (a) through 
a mass ; (6) through a shaft; (c) through gearing ; (d) and (e) application 
to flexibly-coupled systems. 


(iii) Transference through gearing. 

Let Z A in Fig. 33c be the dynamic stiffness of the system to 
the right of, and including, the gear-pinion, and let Z B refer to 
the system to the right of, and excluding, the gear wheel. If 
the number of teeth on the pinion is n times the number of 
teeth on the gear-wheel, the gear-ratio is »:1. If the displace¬ 
ment of the pinion is 0 A , the displacement of the wheel is n0 A ; 
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whereas if the torque at the pinion is T A , the torque at the 
wheel is T A /n. From the equation (22.2) the relation between 
Z B and Z A follows : 

Zb = ~Z A . . : (22.5 e) 

71 


The same method (of multiplying displacements and dividing 
torques by the gear-ratio) can be used in frequency tables. 

Successive use of the two equations (22.5 a, b) alternately 
enables the dynamic stiffness at the inertia J 4 in Fig. 27 to be 
determined as a function of the frequency, for the dynamic 
stiffness of the (zero) system to the right of J x is zero, and the 
point of reference can be shifted from right to left until the last 
inertia J 4 is included. This affords a third method of calculation 
of the effective inertia of the system. 

The application of the method of dynamic stiffnesses to the 
problem of flexibly-coupled systems is illustrated in Figs. 33d, e. 
Two systems, 1 and 2, are coupled through the shaft C c . The 
dynamic stiffnesses are as shown in the exploded view, Fig. 33c. 
The quantities Z 1A and Z 2B are known as functions of the fre¬ 
quency and it is required to find the natural frequencies of the 
coupled system. The stiffness Z lB is given in terms of Z XA by 
(22.56), i.e. 



and the condition for natural frequencies is 


i.e. 


%1B + %2B — A 

--L.-L+JL 

ZiB Z lA C c 


( 22 . 6 ) 


A striking example of the power of the dynamic stiffness 
method, and the convenience of its use, is found in the problem 
of systems coupled by single helical gearing (see reference 15 in 
the Bibliography, page 151). 

The method of dynamic stiffnesses is capable of extension to 
damped systems, and provides the easiest way of calculating the 
resonant frequencies of systems with large damping forces. The 
full details of the method are given in reference 4 in the Biblio¬ 
graphy, but a brief outline is not out of place here. It has 
already been shown (Chapter II, section 9) that in the presence 
of damping forces the displacement in forced motion has a non- 
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zero phase-angle with respect to the applied force. Thus if a 
torque T .sin cot is applied at a point P in a damped system, 
the resultant displacement is of the form 

6 = A.ain (cot+<f>) 

= .4.cos <^> sin a>£+4.sin <f> cos cot. 

In Appendix I, section 35, it is. shown that the cosine function 
can be expressed as : 

cos cot = j. am cot, 

j being a versor operator turning the associated vector through 
an angle n/2 in the positive sense. (This conception has already 
been utilised in the discussion of the damped seismic vibrograph, 
Chapter II, sectiop 12, example II.) Thus the displacement 0 can 
be written : 

6 = (A. cos ^ +Aj. sin <f>) sin cot. 

The dynamic stiffness at P is defined, as before, as the ratio of 
the torque to the displacement at P, but to emphasise the complex 
nature of the quantity the small letter z is used. Thus 

T 

z = -. 

A(coa<j>+j.ain <f>) 

Multiplying top and bottom of this fraction by (cos <f> — j .sin <f>), 
T 

z = -j(cos <f> —j.am<f>), as cos 2 <f> —j 2 avar<f> = 1. 

A. 

T T 

If a = -j cos <f> and /? = — sin <f>, 

JL jA. 

then z = a — jft . . . . (22.7) 

Formulae similar to (22.5) can be evolved for transferring the 
reference point through a system including dashpots which 
produce forces proportional to velocities. The “ modulus ” of 
the dynamic stiffness is defined as 

\z\ = Va 2 +£ 2 . . . (22.8) 

which, by the definitions of a and /?, can be written as 

1*1=2 • • • • (22.9) 

and clearly the amplitude of motion A is a maximum when the 
modulus of the dynamic stiffness is a minimum. In order to 
determine the resonant frequencies of a damped system it is 
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necessary therefore to obtain a curve of the modulus of the 
dynamic stiffness plotted against frequency; the resonances 
correspond to the minima in this curve. 

Note .—The present chapter has been concerned primarily with torsional 
vibration, as the multi-mass systems occurring most frequently in engin¬ 
eering are torsional ones, and all the results can be converted for use in 
linear systems by means of the table of analogous quantities given in 
Chapter I. It is no part of this present work to consider how the inertias 
and stiffnesses of, say, an internal combustion engine are calculated; 
in this work it is assumed that such data are available and to hand. Full 
details of the methods of calculation of these dynamic constants are 
given by Dr. Ker Wilson—see reference 6 in the Bibliography at the 
back of the book. 

Summary. 

The use of effective inertia and dynamic stiffness facilitates 
the determination of the natural frequencies of complicated tor¬ 
sional systems. As demonstrated on page 63, the resonant 
frequencies are very nearly equal to the natural frequencies 
calculated on the assumption of zero damping. Thus the first 
part of the practical problem of avoiding resonance is solved if 
the natural frequencies are determined. 

The effective inertia of a system at any point is the amplitude 
of the alternating torque which must be applied at that point 
to maintain unit amplitude of vibratory acceleration there. The 
dynamic stiffness is similarly defined, “ displacement ” being 
substituted for “ acceleration.” Natural frequencies are the 
frequencies for which the effective inertia or dynamic stiffness 
is zero. 

If two systems are joined together, the natural frequencies 
of the combination are those frequencies for which the total 
effective inertia (or dynamic stiffness) is zero ; this total is the 
sum of the effective inertias (or dynamic stiffnesses) of the two 
systems at the junction. The method of dynamic stiffnesses can 
also be extended to deal with heavily damped systems, and 
systems joined by flexible shafts or couplings or by gearing, both 
spur and helical. 

EXERCISES IV 

(Numerical examples on the calculation of effective inertia and dynamic 
stiffness of systems are rather tedious unless a calculating machine is 
available. Such examples have, therefore, been omitted from this set 
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of exercises, except for No. 1, the solution of which does not require more 
than slide-rule accuracy and for which the majority of the work has 
already been done in the text.) 

1 • Calculate- the dynamic stiffness of the system of Fig. 27 from the 
figures of effective inertia given in Table III, section 18. 

Suppose that two systems similar to that of. Fig. 27 are coupled by 
means of a shaft joining the two inertias J v If the stiffness of this shaft 
is 1,000,000 lbs.ins. per radian, determine the natural frequencies of the 
complete system by means of (22.6). 

2. Obtain the natural frequencies of the system of Fig. 16 by calcu¬ 
lating the effective inertia at m % by the method of Section 18, and also 
by calculating the dynamic stiffness at m a by the method of Section 22. 

3. Obtain the resonant frequency of the system of Fig. 9 by deriving 
an expression for the modulus of the dynamic stiffness at the mass (see 
Section 22). 



CHAPTER V 

CONTINUOUS SYSTEMS 

(Heavy Shafts and Beams) 

Introductory. 

All the vibrating systems considered in the previous chapters 
have only a finite number of degrees of freedom, and hence of 
natural frequencies. The present chapter is concerned with 
continuous distributions of mass and flexibility such as" occur in 
heavy shafts and beams. 

All material structures have a continuous distribution of mass 
and flexibility; a cantilever spring, for example, should be 
regarded not only as a flexible unit but also as a quantity of 
matter. In many cases, where the mass suspended on such a 
spring is very heavy in,comparison with the spring, it is convenient 
to neglect the inertia effects of the mass of the spring and to 
consider it to be “ light,” so that the system is essentially the 
same as that in Fig. 1 or Fig. 9. Simplification in this manner, 
however, is not always possible; when the mass of the spring 
cannot be considered negligible in comparison with that of the 
suspended body, it is necessary to study the effects of this con¬ 
tinuous distribution of mass in the spring. 

It is frequently important to be able to calculate the natural 
frequency of systems including units such as diaphragms, rings, 
plates, etc. It is not possible here to describe in detail the 
analysis of such systems; the reader is referred to Timoshenko 
(reference C in the Bibliography, page 151) for a full treatment 
of these matters. 

In this chapter the torsional vibrations of a heavy shaft of 
uniform circular cross-section are studied by the aid of an 
extension of the dynamic stiffness method. Heavy beams are 
treated in normal fashion, and also by Rayleigh’s approximate 
method, which is based on a consideration of energy. Rayleigh’s 
method enables the lowest natural frequency of loaded heavy 
beams, and similar systems, to be determined in an easy 
manner. 


80 
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23. Torsional vibration of a heavy circular shaft. 

The torsional systems discussed in Chapter IV are composed of 
rigid inertias connected by light shafts, i.e. the flexibility of the 
inertias and the polar inertias of the shafts have been neglected. 
Many torsional systems encountered in practice can be con¬ 
sidered as built-up in this manner; in internal combustion 
engines, for example, torsional vibration of the crankshaft pro¬ 
duces vibration of the reciprocating parts (pistons, connecting- 
rods, etc.) and the polar inertia of the crankshaft itself is very 
small compared with the inertia at each cylinder due to the 
motion of these reciprocating parts. In the present section the 
torsional vibration of a heavy uniform circular shaft will be 
considered. 

Fig. 34 shows longitudinal and transverse sections of such 
a shaft, the length being l and the radius r. Torsional vibrations 


*0 


Fig. 34.—Uniform heavy circular shaft. 



take place about the longitudinal axis AB, and it is assumed that 
during the motion plane cross-sections of the shaft remain plane 
and undistorted. Consider a thin slice normal to the axis AB 
at a distance x from one end B, and of thickness dx. Let the 
angular displacement of the cross-section at x be 0 from the 
static position, and let the angular displacement of the cross- 
section at x+8x be 6+80. The twist over the length 8x of the 
slice is then 80, and dx and 80 are connected by the formula 


80 = T_ 
8x 01 ’ 


. (23.1a) 


where T is the torque transmitted by the slice, O is the shear 
modulus of the material, and I is the polar moment of area of the 
cross-section, defined thus : 

I = Er 2 8a 


where 8a is an elemental area of cross-section at a distance r 
from the axis AB, and the summation is extended over the 
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whole cross-section. (See Appendix II, section 39.) The polar 
inertia of the slice about AB is 


dJ=Ipdx . . . (23.16) 

where p is the density of the material in appropriate units; in 
the engineers’ system of units, dx is in inches, I in inches 4 , p in 
slugs per cubic inch and SJ in lb.in.sec. 2 (see Appendix II, section 
37). The difference 8T between the torques T and T+6T at 
either end of the slice maintains the acceleration B. Thus 


From (23.16, c) 


8T = BdJ . 


(23.1c) 


the partial notation being used as there are two variables x and t. 
Differentiation of (23.1a) gives : 


and hence 


?Z = GI— 
dx ~ U1 dx 2 


B = - — 

p dx 2 


(23.2) 


Equation (23.2) is a partial differential equation, and 0 is 
evidently a function of both x and t. Assume that 0 is such 
that it can be written 

d = X.f(t). . . . (23.3) 

* 

where X is a function of x only and f(t) is a function of t only. 
Then 


and 


*- r 5“ wl 

te 2 JK) dx* 


Substitution in (23.2) gives 




S’X 


Now, as X is independent of t, so must be 

Hence g(JWI - [J g j]/(0. . . ^.4) 
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and the expression in square brackets is independent of t. 
this expression by — co 2 , so that 


6WZ_1_ 
p dx 2 X 


— to 2 


Denote 

(23.5) 


Then 




which yields the normal simple-harmonic-motion result 

f(t) = a.sin (wt+<l>) . . . (23.6a) 

Equation (23.5) can be written 


where 


d*X 
dx 2 



. (23.66) 


the solution of which is evidently of the form 
X = 6.sin ( cokx+y>). 

It is more convenient in this case to write X in the form 

X = 6j sin mkx +6, cos cokx . . (23.6c). 

From equations (23.3) and (23.6), 

0 = (6j sin a)kx+b t cos a>fcc)a.sin (cot+<f>) 

or, writing A for b x a, and B for 6*a, 

6 — (A .sin cokx+B.cos cokx) sin (a)t+<f>) . (23.7a) 

The torque T is given by (23.1a) and (23.7a) as: 

T = GIook( — B. sin mkx + A . cos cokx) sin (cot+<f>) (23.76) 

From the equations (23.7) the frequency of vibration co/2n (i.e. 
the natural frequency) can be determined in accordance with 
specified end-conditions. 


24. Dynamic stiffness and natural frequencies of heavy 
shaft. 

Without loss of essential generality the phase-angle ^ can 
be made zero by a suitable choice of the instant t = 0. In the 
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absence of damping forces, the torques and displacements at the 
ends A and B of the shaft can be written: 

(^j — si® tot at B (x = 0) 

s™ tot at A (x — l) 

Putting x — 0 in (23.7) yields 

P = B . . . (24.1o) 

T b = GIcokA . . . (24.16) 

and putting x = l, and writing y for cokl, 

a = A.sin y-\-B.cos y . . . (24.1c) 

T a = GI(y/l){ — B.ainy+A.cosy) . (24.1d) 

Solving these last two equations for A and B, 

A = a sin y+KT A cos y . . (24.2a) 

B = — KT a sin y+a cos y . . (24.26) 

where K = -J— 


The dynamic stiffnesses Z A and Z B at A and B are 


Z A = 


Zr — 


i.e. Z B = GIcokr= 

X) 

_ qjV « sin y+KT A cos y 
l ’ a cos y—KT A sin y 

This last formula reduces to 


Z D = 


cot y+ 


cosec y 

cos y — KZ a sin y 


(24.3) 


where K = —— 

Gly 

y = cokl 
and k 

Equation (24.3) gives the dynamic stiffness at B in terms of 
that at A ; it can be applied immediately to systems similar to 
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those discussed in Chapter IV except for the inclusion of heavy 
shafts, the formula (24.3) being used in place of (22.5 b). 

Four special cases are worthy of particular consideration: 

(i) Z at one end of the shaft, the other end being free (Fig. 35a). - 



Fro. 35.—Various end conditions and special cases for Fig. 34. 


In this case Z A — 0, 

j „ 1 ( , , cosec wl 

and Zg = -=< — cot y-\ - 

K\ cos y J 

= Itany 

Z* = Vi- 


Hence 


(24.4a) 
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(ii) Natural frequencies of shaft, both ends free (Fig. 356). 
At the natural frequencies, Z B = 0 


Hence 


0)1 J I = °’ *’ 


2n, etc. 


Apart from the zero solution, the natural frequencies are 

a)/2n 


where 



(24.46) 


(iii) Z at one end of the shaft, the other end being clamped 
(Fig. 35c). 

In this case 


and 

Hence 


Z A = co 

Z B — ~( — cot y) unless sin y = 0. 

A 



(24.4c) 


(iv) Natural frequencies of shaft, one end clamped (Fig. 35 d). 
At the natural frequencies Z B = 0 


Hence 


- 1 


371 5 71 

T’ T’ etc - 


or 


co = 0. 


Apart from the zero solution the natural frequencies are 

co/2n 

i 7t jO 3tt fO 5tt Jo , ic%A aj\ 

where ” = aVj' • (24 '* ,) 


Three other special cases, illustrated in Fig. 35e, /, g, are of 
interest. The method of calculating the natural frequencies of 
these systems is indicated in Exercise 4 at the end of the chapter. 


25. Flexural vibration of uniform beam. 

Fig. 36 represents an elemental portion of a uniform heavy 
beam vibrating transversely. The section is at a distance x 
from one end of the beam and is of length dx ; the total length 
of the beam is l. It is assumed that each point in the beam 
moves at right-angles to the static direction of the beam-axis 
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Ox. Let the displacement of the section from the static position 
be y. Elementary theory of elasticity (see Appendix II, 
section 38) gives the formulae: 


M — Bending moment = 

OX* 


F = Shear force 


BM 
Bx J 


(25.1) 



Fig. 36.—Elemental portion of vibrating beam. 

where E is Young’s modulus of elasticity for the material and 
I is the second moment of area of the cross-section defined as 
Zr 2 da, da being an elemental area of cross-section at a distance r 
from the neutral axis of bending (see Appendix II). 

Let the mass per unit length of the beam be m. The resultant 
shear force on the section is — 8F in the positive direction of y ; 
hence 

dF = — mdxy 
3F 

or Tx = ~ my 

-o . f /oki\ VF 3 2 M WT 3*y 

But from (25.1) ■ _ = 

as the quantities E and I are constant over the whole length 
of the beam. 

Hence EI^ - my . . . (25.2) 

OX 

(25.2) is a partial differential equation similar to (23.2) except 
that the partial differential with respect to a; is the fourth instead 
of the second. The solution y can therefore be written 

y = X. sin cot . . . (26,3) 
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, where X is a function of x only, i.e. X is independent of t. The 
phase-angle required for the most general solution is omitted as 
it is of no, importance here. Substituting in (25.2), 

EI^--=mco*X 


or 

where 


^ = k*X 
dx* 


k* = 


moo* 


EIJ 


(25.4) 


Writing equation 25.4 in operator form, 

D*X = k*X, where D = A. 

OX 


Hence D = ± k or ± ik (i = V— 1). 

The general solution in x must have four arbitrary constants, 
as (25.4) must be integrated four times to eliminate differential 
coefficients. By extension of the method of Section 2, Chapter I, 
the general solution is found to be of the form: 

X - C 1 e kx +C i e- kx +C 3 e ikx +C t e- ikx . (25.5) 

It is shown in Appendix I, sections 33 and 34, that 

e e = cosh 0 +sinh 6 
e~ 9 = cosh 6 — sinh 0 
e ie = cos 0+i.sin 0 
e~ ie = cos 0— i.sin 6 

Putting kx for 6, (25.5) can therefore be written in the form 
X — A ± sin kx-\-A 3 cos kx+A 3 sinh kx-\-A x cosh kx 
and hence 

y — (A t sin kx-\-A % cos kx+A 3 sinh kx-\-A x cosh kx) sin cot (25.6a) 

The bending moment M is given by (25.1) and (25.6 a) as: 

% 

M = EIk 2 (— Ai sin kx — A % cos kx+A 3 sinh kx 

+A t cosh kx) sin cot (25.6 b) 

and the shear force F = i.e. 

ox 

F = EIk 3 (— Ai cos kx+A t sin kx+A t cosh kx 

+■ A t sinh kx) sin cot (25.6 c) 

From equations (25.6) the natural frequencies can be determined 
aocording to the end-conditions. 
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26. Natural frequencies of beam. 

(i) Cantilever. 

If the beam is built-in at one end, as in Fig. 37a, the end- 
conditions are that at x = 0, the displacement and slope are zero, 



x’o (a) 



Fig. 37.—(a) Cantilever ; (6) free-free beam. 

while at x = l the bending moment and shear force are zero, as 
there is no loading beyond the end of the beam. 

I.e., at x = 0, y — 0 and ^ = 0 

ox 

at x = l, = 0 and ~ = 0. 

dx* dx 3 

Substituting these values in (25.6), 

A 3 +A t =0 . . (26.1a) 

A 3 +A 3 = 0 . . . (26.16) 

— A x sin kl — A 3 cos M-\-A 3 sinh kl+A t cosh kl = 0 (26.1c) 

and , 

— Ax cos sin kl+A 3 cosh kl+A 3 sinh kl — 0 (26.1d) 

Elimination of the constants A lt etc., from these four equations 
leads to the result: 

cos kl .cosh kl +1 = 0 •. . . (26.2) 

The first four roots of this equation are 

kl = 1-875, 4-694, 7-855, 10-996 
and the first four natural frequencies are therefore: 

a)/2n, 

where a> = and A = 3-52, 22-4, 61-7, 121 . (26.3) 
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(ii) Free-free beam. 

If the beam is free at both ends, as in Fig. 376, the end- 
conditions are : 

M = F = 0 at both ends, 

i.e. ^ = 0 at both ends lx = 0 and x — l) 

ox 2 dx s 

Hence — A a +A t = 0 . . (26.4a) 

-^+^3 = 0 . . . (26.46) 

— A x sin kl — A a cos kl+A 3 sinh kl+A a cosh kl = 0 (26.4c) 

and 

— A t cos kl-\-A a sin kl-\-A 3 cosh kl-\-A t sinh kl — 0 (26.4d) 

Elimination of the constants A u etc., from these equations leads 
to the result: 

cos H.cosh kl — 1 = 0 . . . (26.6) 

The first four roots of this Equation are 

kl = 4-730, 7-853, 10-996, 14-137 
and the first four natural frequencies are therefore 

c o/2n , 

where o> = A and A = 22-4, 61-7, 121, 200 . (26.6) 

The natural frequencies of beams under other end-conditions 
can be found in a similar manner; the subject is treated fully 
by Timoshenko (reference 7 in the Bibliography at the back of 
the book), who also discusses the validity of the various assump¬ 
tions that have been made. 

The deformation shapes of the cantilever, Fig. 37a, at the 
various natural frequencies are illustrated in Fig. 38. It will be 
seen that the continuous beam has the same property as multi¬ 
mass systems as regards the different modes of vibration, namely 
that each successive natural frequency involves an additional 
node. 

The analytical solution of the problem of the uniform beam 
is relatively simple, owing to the fact that the quantities E and 

dF 

I in equations (25.1) and m in the equation — = — my are 

dx 

constant over the whole length of the beam. In a beam of 
variable cross-section or of non-uniform material the quantities 
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E, I and m are not all constant, and the differential equation 
corresponding to (25.2) is 



The difficulty of solution is not due to any abstruseness of the 
theory but to the complexity of. the calculation. The equation 





Fig. 38.—Deformation shapes of cantilever at first three natural frequencies. 


(26.7) cannot be solved in terms of analytic functions if E, I and 
m are unspecified functions of x , and an approximate numerical 
method has to be developed to enable a solution to be obtained 
at all. Numerous authors have published methods of solution; 
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an adequate treatment of the subject is given by Timoshenko 
(reference 8 in the Bibliography). 

27. Rayleigh’s method for heavy beams. 

Rayleigh’s approximate method for finding the lowest natural 
frequency of heavy uniform beams depends upon the energy 
properties of vibrating systems. When the system is at its 
maximum displacement the kinetic energy is instantaneously 
zero, as the system is instantaneously at rest, and the potential 
energy is at its maximum value ; for the potential energy is the 
work done against the restoring forces, and this is clearly a 
maximum at the maximum displacement. Similarly, when the 
system passes through its mean position (zero harmonic displace¬ 
ment) the kinetic energy is a maximum and the potential energy 
is zero. 

Consider for example a single-mass system, the displacement 
of the mass m being x = a .sin cot against a spring force kx, as 
in Chapter I. The kinetic and potential energies are given by 

K.E. = \mx 2 = \ma 2 co 2 cos? cot (a)l 

P.E. = J (kx)dx = \kx 2 = \ka 2 sin 2 to t (6)J 

At the maximum displacement sin cot = ± 1, cos cot — 0 
Thus (P.E.) maa; , = \ka 2 , K.E. = 0 . . (27.2 a) 

At the mean position sin cot — 0, cos cot = ±1, 
and (K.E.) max _ = \ma 2 co 2 , P.E. = 0 . (27.26) 

The total energy in the system is the sum of the kinetic and 
potential energies, and this total remains constant throughout 
the motion in the absence of damping forces, as there is no 
dissipation. The maximum kinetic energy is thus equal to the 
maximum potential energy, and thus 

\ka 2 = \ma 2 co 2 
k 

whence to 2 = — . . . . (27.3) 

m 

and the frequency co/2n is thus determined. 

In the case of heavy beams, the expressions for the kinetic 
and potential energies are more complicated. Consider the 
elemental section of beam in Fig. 36. For sinusoidal motion 
y = y 0 sin cot the maximum kinetic energy 6(K.E. L,„ r of the 
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section is i(K.E.) maa . = \mu> 2 y n 2 6x, as the mass of the section 
is mix. The maximum kinetic energy of the whole beam is 
therefore 

(K.E.) maa . = £mei> 2 f y 0 2 dx . . (27.4a) 

Jo 

Referring to Fig. 39a, if the left-hand end of the section is fixed 
the bending moment M turns the right-hand end through the 
angle <30, and M is proportional to dd. The work done against 
M during this rotation is given by the shaded area in Fig. 396, 



Fig. 39.—Derivation of potential energy formula for beam. 


and is equal to \M6Q. The potential energy of the section is 
thus 

<5(P.E.) = pf<50. 

The slope of the displacement curve at the left-hand end of the 
dy 

section is and the slope at the right-hand end is 

OX 


and 


i + P)- fa ' hen “ * - 

«(P.E.) _ pfgfa. 


dx 2 


Substituting M = 

OX 

a<p.E.) = 

The maximum potential energy of the whole beam is therefore 

(P.E.W = \EI^(^jdx . . (27.46) 

The natural frequency a>/2n can be found if the deformation 
shape is known, by equating the two expressions (27.4). The 
procedure is to guess a deformation curve which satisfies the 
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end-conditions and corresponds to the mode of vibration con¬ 
cerned. As an example, consider the system of Fig. 40, consisting 
of a heavy uniform beam simply supported at each end (so that 
the displacements at the ends are zero but the slopes are not 


Fig. 40.—Simply-supported beam with central loading. 



necessarily zero) and carrying a concentrated load at the centre 
of the span. 

Let l — length of beam 

m = mass per unit length of beam 

M b = mass of beam 

M l = mass of central load 

y c = maximum deflection at centre of span. 

The end-conditions to be satisfied are that the displacements and 
bending moments at each end are zero, i.e. 

dhj 

y = 0 and —^ = 0 at x = 0 and x = 1. 
a dx 2 

For the fundamental, or lowest, natural frequency the curve 

TtX 

Vo = y c sin -j- satisfies the end-conditions and the mode. From 
(27.4a) the maximum kinetic energy of the beam is 

Cl TTCTj 

(K.E.) roax . = £mw 2 | y c 2 sin 2 —dz = $m<o*y c H 

= i Mtfotyc* 

and for the central load (K.E.) mM = \M I p % y c ‘ t . 

Hence for the whole system, 

total (K.E.) maa . = . (27.5a) 

From (27.46); 

(P.E.W = 

The value of the definite integral is l, and hence 

(P.E.W = ®|2^ . . (27.56) 
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Equating the expressions in (27.5 a, b), 


whence 


\«> 3 y c \\M B +M L ) = 

2 _ 48-7.0/ 

" l\\M B +M L ) ’ 


(27.6) 


From this formula the frequency (o/2n is obtained. 

Two special cases deserve attention: 

(i) If the central load is zero, the system reduces to a heavy 
beam simply supported at both ends, and (27.6) becomes 

, 97-4.0/ 

“ M B l* ' 


This result agrees with the exact expression derived by the 
method of Section 26. 

(ii) If the beam is light compared with the central load, 
(27.6) becomes 

. 48-7 El 

01 M l V » ' 

Comparing this formula with a> 2 = Te/m for the system of Fig. 1 
the effective spring-constant of the light beam is here given as 
48-7 EI/l 3 . If this spring-constant is calculated according to 
normal methods (see Appendix II, section 38) the numerical 
factor is obtained, as 48 ; thus Rayleigh’s method gives a value 
for a) 2 which is 1-4% in excess of the true value ; the frequency 
error is thus only 0-7%. 

Rayleigh’s method is very useful when it is desired to know 
the effect of neglecting the mass of a beam carrying concentrated 
loads; in the case of a simply-supported beam carrying a central 
load, (27.6) shows that half the mass of the beam should be 
added to the central load if it is required to find the lowest 
natural frequency assuming the system to be simply a one-mass 
system similar to Fig. 1, k being the effective stiffness of the 
beam at the centre of the span. 

Table VIII gives various formulae and results for three 
different types of loaded heavy beams : 


(i) cantilever with end-loading, 

(ii) simply-supported beam with central loading. 
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LENGTH OF BEAM l, DEFLECTION y , MASS OF BEAM Mg 
MASS OF LOAD Mu, MAX: DEFLECTION A 


SYSTEM: 






3-°*k- 

& 



■°*fe 



ASSUMED 

FUNCTION 


( KE W. 

< PE W 


c o z 


FREQUENCY 

FOR LIGHT F=— |^- 
BEAM ^ 


CORRECTION 

FOR MASS 0 226 M B O S M B 
OF BEAM 


FREQUENCY 

ERROR 


y * A 

ya^(|-COS^ 

«W(iM^iMa) 


487 El 

195 El 

^(M t +4M e ) 





075 % 
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(iii) clamped-clamped beam with central loading. The table 
gives: 

(а) end-conditions, 

(б) Rayleigh function, 

(c) maximum kinetic and potential energies, 

(d) value of co 2 , 

(e) frequency for light beam (M B = 0), 

(/) correction for mass of beam, 

and (g) approximate frequency error. 

Summary. 

Vibrations of heavy beams and shafts are analysed by obtain¬ 
ing the general solution to the partial differential equations of 
motion, and determining the constants from a consideration of 
the end-conditions relevant to the problem. 

Rayleigh’s approximate method for beams consists of assuming 
a deformation shape for vibration at the lowest natural frequency, 
and then calculating the maximum values of kinetic and potential 
energies for sinusoidal vibrations in this mode. Since when the 
kinetic (or potential) energy is a maximum, the potential (or 
kinetic) energy is zero, and the sum of the two, being the total 
energy content of the system, is constant in the absence of 
damping forces, the maximum values of the kinetic and potential 
energies may be equated. This equation yields the value of the 
natural frequency. 


EXERCISES V 

1. Check by Rayleigh’s method the expression given for the lowest 
natural frequency of a free-free beam in (26.6), and so obtain data for 
another column in Table VIII. 

(The appropriate function is y 0 =* a — Y sin in which a has to 

be found by determining the condition for zero total inertia force (my) 
for the whole beam.) 

2. Check the data given in columns (i) and (iii) of Table VIII. 

3. Draw the swinging forms of a uniform heavy circular shaft at its 
first four natural frequencies, the shaft being free at both ends. 

(Use the formula 23.6c in conjunction with the special end-conditions.) 

4 . Determine the natural frequencies of the heavy-shaft systems 
shown in Fig. 35c, /, g. 

(Use the formula 24.3, putting in the particular end-conditions required 

H 
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(e). Z A = 0, Z B « Jo)* giving y = ^ tan $/, 

where y = f and «7, is the polar inertia of the shaft = pZJ. 

(/) ^4 = oo, = Jco 2 giving y = -—coty. 

(9) Z A = — J^ca 2 , Z 5 = J 2 o 2 , giving 

iC 2 JiJ r a a) 4 +iC(Ji+J r 2 )co 2 c o t2/ — 1=0 
where iC = l/QIy.) 

5. Some variable-pitch aircraft propellers can be “ feathered,” i.e. their 
blades can be turned edge-on to the direction of flight. Idealising the 
blades in the form of uniform beams of rectangular cross-section, and 
supposing these to be built-in to a shaft of negligible diameter, derive an 
expression for the torsional dynamic stiffness of the propeller as a whole, 
the blades performing flexural vibrations in a direction perpendicular to 
the shaft. 

(In equations 25*6, page 88, substitute the end conditions : at x = l 
(outer end of blade), M = F = 0; at & = 0 (root end of blade) y = 0, 
(dy/dx) = A .sin cot, M — T .sin cot. The required dynamic stiffness is by 
definition nT/A, where n is the number of blades.) 



CHAPTER VI 

COMPLEX VIBRATIONS 
(Fourier Analysis) 


Introductory. 

Attention has been confined, in the previous chapters, to 
vibrations of sine-wave form. In a great number of practical cases 
vibrations of this type are of considerable importance, but it very 
frequently happens that the form of serious or even dangerous 
vibrations is not that of a sine-wave. Moreover, it is found that 
many vibratory forces are not sinusoidal (for example, the gas- 
pressure variations in an internal combustion-engine cylinder, or 
the regular blows of a power hammer) and it becomes necessary 
to study their effect. 

The practical result of Fourier’s Theorem (see Section 28 
p. 100) is that any periodic variation whatever can be expressed 
in the form of a series of sinusoidal components. Since a sinusoidal 
force causes a similar vibration at the same frequency as itself, 
a periodic force which is not sinusoidal can cause vibrations at 
more than one frequency. Its components, which are generally 
termed “ harmonics,” can have frequencies which are any 
multiple of its own frequency, the presence and relative severity 
of each harmonic depending upon the form of the periodic 
variation. Thus, for example, the periodic variations in Figs. 41 
and 42 are such that the corresponding Fourier series contain 
only odd harmonics, whereas that in Fig. 43 consists of the sum 
of two sine-waves whose frequency-ratio is 5:6. 

The practical importance of Fourier’s Theorem is great. If 
a structure is acted on by a sinusoidal force with a frequency 
of 2,800 C.P.M., it will only undergo sustained vibrations at this 
frequency, and these will be small if this value does not nearly 
or quite coincide with any natural frequency of the structure; 
but if the pattern of the force is similar to that in Fig. 416, the 
load being suddenly applied and suddenly removed, there are in 
effect sinusoidal forces at 2,800, 8,400 (= 3 X 2,800), 14,000 
(= 5 X 2,800), etc., c.p.m., and if any of these frequencies nearly 
or quite coincides with any natural frequency there may be severe 
vibration. 
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The practical details of Fourier analysis cannot be described 
here at length, and in the present chapter attention is confined 
to the simpler parts of the subject. The author has described 
elsewhere the various methods of Fourier analysis, and particu¬ 
larly of waveform analysis (the analysis of periodic variations 
into their sinusoidal components, the data being a vibration 
record or similar waveform); see reference D in the Bibliography, 
page 151. 


28. Fourier series. 


If a function /(<) of any variable t is periodic with a ^-period r, 
so that if the variable t is increased by any multiple of x the 
value of the function is unaltered (i.e. f(t) = f{t+kx), where k is 
any integer), then the function can be expressed in the form : 

. 2nt . n 2nt , „2 jvt . 

f(t) = a„-|-ai cos- \-a t cos 2-|-as cos 3—, etc. 

XXX 

, , . 2 nt , , . „2 nt . , . „2 nt , 

+ &i sin —-fo* sm 2—+o 3 sin 3—, etc. 
xxx 

The series can be expressed concisely as : 


00 00 

f(t) = a t+ Ya n cos Yb n sin rfoj? 

x r 


where n is integral . . . (28.1) 

The result is known as Fourier’s Theorem, and the series is 
commonly termed a Fourier series. 

If the variable t represents time, so that the function has 
a time-period r, it is convenient to write co for 2jt/x, and then 

00 00 

f(t) = «»+ 2> n cos not -h y ^ b n sin ncot . (28.2) 

n=l n=l 

•The terms a n cos ruot-\-b n sin ncot can be written as r n sin n(a>t-\-<f> n ), 
where r n 2 = a n 2 -h6 n 2 and tan n<j> n = a n /b n so that the series can 
be put in the alternative form: 

00 

f(t) =o 0 + 2J r n sin n(a>t+tf> n ) . . (28.3) 

n=l 

The quantity <o/2n, being the reciprocal of the time-period r, 
is turned the fundamental frequency of the function ; the term 
in the series for which n = 1 is called the fundamental component, 
and the other terms (n = 2, 3, etc.) are known as the harmonics. 
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As an example of the synthesis of a periodic function by the 
addition of sine-waves, consider the function 
/(0) = sin 0+J sin 30+^ sin 50+ . . . 

. . . 2^—- sin (2 n — 1)0 + . . . (28.4) 

In Table IX are listed values of the various terms in (28.4) as 
far as the term involving 70, i.e. the 7th harmonic. 


Table IX 


e/n 

sin 0 

1 sin 3 0 

1 sin 56 

\ sin 76 

sum 

0 

00000 

0-0000 

0-0000 

0-0000 

0-0000 

01 

0-3090 

0-2697 

0-2000 

0-1156 

0-8943 

0*2 

0-5878 

0-3170 

0-0000 

- 0-1359 

0-7689 

0-3 

0-8090 

0-1030 

- 0-2000 

0-0441 

0-7561 

0-4 

0-9511 

- 0-1966 

0-0000 

0-0840 

0-8385 

0-5 

1-0000 

- 0-3333 

0-2000 

- 0-1429 

0-7238 

0-6 

0-9511 

- 0-1966 

0-0000 

0-0840 

0-8385 

0-7 

0-8090 

0-1030 

- 0-2000 

0-0441 

0-7561 

0-8 

0-5878 

0-3170 

0-0000 

- 0-1359 

0-7689 

0-9 

0-3090 

0-2697 

0-2000 

0-1156 

0-8943 

10 * 

0-0000 

0-0000 

0-0000 

0-0000 

0-0000 

11 

- 0-3090 

- 0-2697 

- 0-2000 

- 0-1156 

- 0-8943 

1-2 

- 0-5878 

- 0-3170 

0-0000 

0-1359 

- 0-7689 

1*3 

- 0-8090 

- 0-1030 

0-2000 

- 0-0441 

- 0-7561 

1-4 

- 0-9511 

0-1966 

0-0000 

- 0-0840 

- 0-8385 

1-5 

- 1-0000 

0-3333 

- 0-2000 

0-1429 

- 0-7238 

1-6 

- 0-9511 

0-1966 

0-0000 

- 0-0840 

- 0-8385 

1-7 

- 0-8090 

- 0-1030 

0-2000 

- 0-0441 

- 0-7561 

1-8 

- 0-5878 

- 0-3170 

0-0000 

0-1359 

- 0-7689 

1-9 

- 0-3090 

- 0-2697 

- 0-2000 

- 0-1156 

- 0-8943 

20 

0-0000 

0-0000 

0-0000 

0-0000 

0-0000 

21 

0-3090 

0-2697 

etc. 

0-2000 

etc. 

0-1156 

0-8943 


The values in the last column of Table IX are plotted in Fig. 41a. 
If more terms of the series are taken, the curve approximates 
more nearly to the shape of Fig. 416,* i.e. 

/(0) = jt/ 4 for 0 < 0 < 2ji \ • „ . 

/(0) = - *r/4 for n < 0 < 2nj ' ' K ) 

It will be seen that the function (28.5) is such that 

m = -/(- 9 ) « -/( 2 * - 0) . . (28.6) 

* See note on Gibbs’ phenomenon, p. 106. 
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so that the curve, Fig. 416, is skew-symmetric about the values 
0 = 0 and 6=3 r. The result is due to the fact that all the 
terms in the series (28.4) are sine-terms and obey the law (28.6). 



-*/ 4 L I-1 L_ 

(C) 

Fig. 41.—(a) Graph of the function (28*4), first four components : (6) sine-term 
function ; (c) cosine-term function. 

Let the variable 0 be changed to Q lt where 0 = 0 4 -f-jr/2 ; the 
resulting graph is that of Fig. 41c, which is seen to be symmetrical 
about the lines 0 r = 0 and 6 l = n, so that 

/(01) =/(- 0i) =/(2 tt - 6,) . . (28.7) 

Making the substitution in the series (28.4), 

/(0,) = cos 0 X - i cos 30j+| cos 50, etc. . (28.8) 

as sin (0 X + n/2) = cos 0 

sin 3(0!+jr/2) = — cos 30! 
sin 5(0i-f 3t/2) = cos 50 x , 
etc. etc. 

The forms of (28.4) and (28.8) point to a general conclusion that 
if the function /(<) in (28.2) is symmetrical about the lines t = C, 
C-\-3i/<o, etc., the series can be put into a form containing cosine 
terms only by making the substitution t x = t — C\ and similarly 
if the function is skew-symmetric about these lines the series 
can be .put into a form containing sine terms only. 

29. Fourier analysis. 

The process of determining the coefficients a and 6 in the 
series (28.2) for any particular function f(t) is known as Fourier 
Analysis. The procedure depends upon the results: 
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sin meat. sin neat dt = 0 


(6) sin www*.cos rmt dt = 0 * mand M are unec l ual 

J o integers 

» r2n/(o 

(c) I COS WWW*. COS ?WW* dt = 0 

Jo 

r2n/<o 

(d) I sin moot, cos mot dt = 0 if m is an integer 
J o 

r2n/(o 

(e) I sin www* dt = 0 
Jo 


cos www* dt = 0 


p2jr/co 

(g) I sin 2 mcotdt = n/w 
Jo 

r2n/ai 

(h) I cos 2 mot dt = n/o 

Jo 


if m is an integer. 



These results are proved in Appendix I, section 36. 

Integrate (28.2) with respect to t over a cycle of the function, 
i.e. between the limits 0 and 2 n/o for t. 


r2n/a> r2n/oa 

I f(t) dt = I ( a o+ /. a n cos / b n sin neat) dt. 

® ® n= 1 n=l 


The equations (29.1) show that every term on the right-hand 
side will be zero except that due to the integration of a Q . Hence 



and a 0 = —l f(t) dt . . (29.2a) 

InJ o 


Now multiply (28.2) by cos meat and integrate over a cycle ; thus 

r2n/<u r2n/(o n 

l f(t) cos mot dt = I a m cos 2 mot dt = —a m 

Jo Jo ® 

as all the other terms integrate to zero by reason of the equations 
(29.1). 

Hence a m = — I f(t) cos mot dt . (29.26) 

71J o 


Similarly, by multiplying (28.2) by sin mot and integrating over 
a cycle, it is found that 
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<u f **/• 

b m = -I /(0 sin moot dt . . (29.2c) 

rcJo 

The procedure for determining the Fourier coefficients a and 
b is thus as follows: 

(i) The constant o 0 is found by evaluating the average value 

of the function over a cycle, for that is what the value 
(29.2a) is. 

(ii) The cosine coefficients “ a ” are found by multiplying 

the function by the appropriate cosine function, 
integrating over a cycle, and multiplying the result 
by ao/ji. 

(iii) The sine coefficients “ b ” are found in a similar manner, 

except that the appropriate sine function is used 
instead of the cosine function. 

The equations (29.2) are collected together below for easy 
reference. 



Fig. 42.—Periodic functions analysed in section 29. 


o° = ^j o (a) I 

0 ) r2n/a) I 

a m = —\ f(t) cos mat dt ‘ (&)> . (29.2) 

n,Jo I 

/» f 2n/<o \ 

b m = — I /(f) sin moot dt ( c) 

n J o J 

As an example, the function graphed in Fig. 42a will be 
analysed in this manner. 
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cos mtdt = 0. 


sin mt dt 


f(t)= 1 for 0 < t < n \ 

/(<) = — 1 for n < t < 2aJ ' ' ' ' 

Here a> = 1. By (29.2a), 

a o = ^[£(1) dt +£“(- 1) dt ] = 0. 

By (29.26), 

1 rn l r2n 

a m = — I cos mtdt — -1 cos mt dt = 0. 

71J o 71 J „ 

By (29.2c), 

If" 1 f 2 " 

b m = - 1 sin mt dt - 1 sin mt dt 

^Jo 

IT — cos mt 1" 1|“ — cos mfV 71 

~~ 7t[_ m Jo m J* 

= —(1 — cos mn). 
mn 

If m is even, cos mn = 1 and b m = 0 

4 

If m is odd, cos mn — — 1 and b m — — . 

m mn 

The Fourier series representing the function of Fig. 42a is 
therefore : 

f(t) = -(sin t+\ sin sin 5t . . .) . (29.4) 

n 3 5 


This result agrees with the previous synthesis (Section 28, 
Kg. 41). 

The Fourier analysis of a discontinuous function gives results 
which are accurate except at the discontinuities, where the 
Fourier function “ jumps ” too far; this is known as Gibbs’ 
phenomenon, an account of which is given elsewhere by the 
author (reference 9 in the Bibliography). 

A slightly more complicated function is shown in Fig. 426. 
By a suitable choice of the point t = 0 the function can be made 
either symmetrical or skew-symmetrical about t = 0, so that the 
series can be made to consist of cosine terms only or of sine terms 
only. Let the time-origin be as shown, so that the function is : 


Of 

/(<) = 1- for 0 < t < n 

31 > 

O Of 

f(t) = _ l+_(f - n) = - 3+ - a<t<2n 

n Jt 


{29.5) 
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Dividing the integration into two stages, as before, and using 
the results: 


and 




cos mt dt — t 


= t 


sin mt 
m 

sin mt 


-f 


sin mt 


m 
cos mt 


dt 


m 


m 2 




sin mt dt = — t 


cos mt 


\ 


— cos mt 


dt 


= - t 


m J m 
cos mt , sin mt 


m 


m l 


the Fourier series is obtained as : 

/(<) = -~( c °s cos cos 5< . . .) . (29.6) 

When it is required to carry out a Fourier analysis of a 
function which cannot be readily expressed in an analytic form 
such as (29.3) or (29.5), the integration is performed numerically 
as a summation. The ordinates, or values of the function, at 
a number of time-values equally spaced over the cycle are mul¬ 
tiplied by the corresponding values of the appropriate cosine or 
sine function, and the products summed over the cycle. The 
results of this process are sufficiently accurate so long as the 
number of ordinates is greater than twice the index-number of 
the highest harmonic present. Thus a 48-ordinate analysis will 
give true results so long as there are no harmonics higher than 
the 23rd present. The numerical work is easy but tedious ; 
in order to reduce it to a minimum, various schemes of com¬ 
putation have been evolved, the most convenient one being that 
due to Runge, a full description of which is given elsewhere 
(reference 10 in the Bibliography at the back of the book). 

A particular example of a complex wave with only two com¬ 
ponents is afforded by the phenomenon of beating . The effect 
occurs when two sine-waves of slightly different frequencies are 
added. Let the two functions be a. sin cot and b . sin (co +A (o)t, 
where Aco is small compared with co. The sum of the two 
functions can be expressed in the form 

(a-f&.cos Aa>t) &ma>t+b.&mAa)t.cd&(ot 9 

or r sin (art +<f>) 9 where r 2 = (a+b.cos Awt) 2 +(b.sin Awt)* 

= a, 2 +b 2 +2ab. gob Acot. 
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The amplitude r varies between the limits (a+6) and (o — 6), 
the frequency of variation being Aa>/2n, which is the difference 
between the frequencies of the two components. As an example 
of this phenomenon, the function sin6$+2.sin 5t is plotted in 
Fig. 43. 

ifM .f(fc) = sin 6t +2sia5t 



Fia. 43.—Boating. 

30. Vibration under non-harmonic forces. 

A sinusoidal force is termed a harmonic force. Periodic forces 
which are not sinusoidal, but which can however be expressed 
as a Fourier series of harmonic forces, are termed non-harmonic , 
and they are of very frequent occurrence in engineering. In the 
internal-combustion engine, for example, the gas-force acting on 
the piston in each cylinder is a periodic function with many 
harmonics present; the frequency of the fundamental component 
is the frequency of firing in each cylinder, i.e. half the number 
of revolutions per minute of the crankshaft in the case of four- 
stroke engines. These gas-forces are transmitted to the crank¬ 
shaft through connecting-rods of finite length, and the non-linear 
relation between the force on the piston and the resulting torque 
exerted at the crank-throw introduces further harmonics. The 
resultant torque in the crankshaft is therefore a Fourier series, 
the constant term (a 0 ) in which represents the steady torque 
available for the power-output of the engine, and the sine and 
cosine terms in which induce torsional vibration of the crankshaft 
and its associated inertias. As the frequency of the fundamental 
component in c.p.m. is | x r.p.m. of the crankshaft, there will 
in general be harmonic torque components at 

(J, 1, 1J, 2, 2£, 3, . . . etc.) x engine r.p.m. 

Many of these harmonics are practically balanced-out by. suitable 
design of a multi-cylinder engine, but all the harmonics are 
present in a single-cylinder engine. Further information on 
this important subject may be found in textbooks on torsional 
vibration; the reader is referred to that by Dr. Ker Wilson 
(reference 11 in the Bibliography). 
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Suppose the mass in Fig. 9 to be acted upon by a complex 
force P = P x sin pt +P S sin qt, where p and q are not equal. 
Using the method of section 9, 

P P 

[D 2 +2yD+co 2 ]a: = — sin pt -\—- sin qt 
m m 


and the displacement x is found to be of the form: 

x = 6 X sin (pt — f 1 )+b t sin (qt — ip t ) (30.1) 

where the quantities b and ip are expressions similar to those in 
(9.3), subscript 1 indicating the frequency p/2n and subscript 2 
indicating the frequency q/2ai. The motion of the mass is 
therefore the sum of two sine functions whose frequencies are 
the same as those of the two components of the exciting force. 

In general, if the exciting force is a Fourier series in the 
form (28.3), i.e. 

00 

si nn(pt+(f> n ) . . (30.2) 

n=1 

then the displacement of the mass is given by 

00 

X = a 0+^ a n sin n(pt+ip n ) . . (30.3a) 

?l=l 

where the amplitude coefficients a n are given in terms of the 
modulus of the dynamic stiffness | z n \ at the frequency np/2n as 



(see Chapter IV, end of Section 22). 


(30.36) 


The phase-angles ip n can easily be calculated, but are generally 
of far less importance than the amplitudes. 

An important consideration is the work done by the force 
(30.2) on the displacement (30.3). Consider any particular 
harmonic of the displacement-function (30.3a), say the mth 

harmonic x m — a m sm.m(pt+f m ). The work done is j F-dx m ; 

-but dx m = a m mp . cos m(pt -\-ip m )dl, and the work done during a 
cycle of the fundamental component is 


W m =l a m rn^F 0 +^\ n sin n(pt +^ n )Jp.cos m(pt±ip m )dt (3Q.4) 


and equations (29.1) show that all the terms on the right-hand 
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side of (30.4) integrate to zero except that due to the mth harmonic 
in the force input. For any particular harmonic in the displace¬ 
ment function, the input energy that is required to make good 
the loss due to dissipation is therefore derived exclusively from 
the corresponding harmonic in the force. This result is of great 
practical importance, for it means that if one harmonic com¬ 
ponent of a periodic force input is varied, only the corresponding 
harmonic in the displacement-function will be altered. 

Summary. 

Equations (29.2) on page 104 enable the sinusoidal components 
to be found for any periodic variation which can be expressed 
mathematically, either as a single function or as a number of 
functions each extending over part of the cycle. These com¬ 
ponents have frequencies which are multiples of the frequency 
of the original variation, the presence and relative intensity of 
each component being dependent upon the form of the variation. 

A form of numerical approximation to the integrals of 
equations (29.2) enables the components of periodic variations 
to be found when these are given as sets of values or as graphs 
(e.g. recorded waveforms). 

Whereas the response to a sinusoidal force is also sinusoidal, 
at the same frequency as the force, the response to a more complex 
periodic force is itself complex ; the response has a component 
at the same frequency as each of the harmonics of the force. 
Consequently, a system with a single degree of freedom, and 
hence one natural frequency, will be in resonance at various 
different fundamental frequencies of a complex force. Thus if 
the force varies after the manner of Fig. 42a, and the natural 
frequency is 1,500 c.p.m,, there will be resonances at frequencies 
of 1,500, 500, 300, 214, etc., force c.p.m. ; in each case, however, 
the large vibrations due to resonance will be at 1,500 o.p.m. 
(See Appendix IV, page 134.) 

EXERCISES VI 

1. Determine the Fourier series for a function /(£) which is such that 
f(t) — a sin kt for 0< « ^ 

f(t) = 0 for t< 2 ji 

k being integral. 
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2. Determine the Fourier series for a function /(<) which is such that 

/(<)=■ 1 . for — 

f (t) = 0 for t l <t<2n — 

being less than n. 

3. Develop a method of Fourier analysis by numerical integration 
with six ordinates, to evaluate the constants r and \p in the function: 

/(<) = o 0 + r i sin (co<+Vi) +»»sin 2 (cot+y,). 
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31. Exponential function. 

/ i\» 

The number e is defined as the limit of the expression (1+- 1 
as n tends to an infinite value. 


e = Lt 

n —>-oo 


Thus 

K)' 


Expanding by the Binomial Theorem : 

n(n — 1) t n{n — l)(n—2) 


e = Lt 

71 - 


. ji+-- 

> l n 


n 


11 


etc. 


i.e. 


e ~ 1+1+ ]4 + || + ]i + etc - 


(31.1) 



The function e x is similarly defined: 

-o(. + i)7=nr 

and the result is obtained: 

e* = l +a; +j| a + j| + j| 4 + etc. . . (31.2o) 


Putting — * for x, 


/y2 /*»3 <»»4 

e-. = 1 - a;+ _-_ + __ etc * 

* A graph of e - * is given in Fig. 44. 


(31.26) 


It is shown in textbooks on analysis that these series are con- 
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vergent for all values of x. The function e x is termed the ex¬ 
ponential function and obeys the normal laws of indices. 

Thus - e x .e v = e x+v , (e x f = «*», etc. 

The exponential function is such that the series form can be 
differentiated term by term (see Whittaker and Watson, refer¬ 
ence 12 in the Bibliography). Performing the differentiation, it 
is seen that the resulting function is identical with the original, i.e. 

iL(e*) = e x ... (31.3) 


32. Sine and cosine series. 

Maclaurin’s formula for the expansion of a function as a 
power series is: 

/(*) =f(0)+xm+£f"(0)+ etc.. . (32.1) 

II 

where dashes denote differentiation with respect to x and /(0) is 
the value of fix) when x — 0, etc. 

If fix) = sin x, /(0) = 0 

f'\x) = cos®, /'(0) = 1, etc. 


In this manner the function sin x can be expanded by Maclaurin’s 
formula as : 


4*3 /y>5 ' /y*7 

• M/ W •!/ . J 

sin x = x — + etc. . 


and similarly, 


11 H 11 


ry* 2 /y> 4 /y»6 

co 8a: = l-(J_ + [4 ~[6 + etc. 


(32.2a) 

(32.26) 


33. Exponential form of sine and cosine functions. 

Let i stand for V— 1; then 
** = — 1, i 3 = — i, i* = 1, i 5 = i, etc. 

Thus «■* _ [l - g+j| - etc.] + ,[* - ||+|£ - «tc.] 

i.e. e te = cos x+i.sin x . . . (33.1a) 

Putting — x for x, 

e -ix = cos x — i.sin x . . (33.16) 
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and hence 


COS X — 


sin x 


e«+e' 


—ix 


-ix 


2 i 


(a) 

(b) 


(33.2) 


From (33.2) it can be shown that cos 2 x + sin 2 x = 1. If the 
co-ordinates of a point in plane Cartesians are x — r. cOs 0, 
y — r.sinfl, the point lies on the circle x 2 +y 2 —r 2 . The sine 
and cosine functions are therefore termed “ circular functions.” 


34. Hyperbolic functions. 

The functions cosh x and sinb x are defined : 


cosh x = 

e x +e~ x 

2 

(a) 

sinh x = 

e x _ e -X 

2 

('b) 


J 


(34.1) 


and it can easily be shown that cosh 2 x — sinh 2 x = 1. If the 
co-ordinates of a point in plane Cartesians are x = r cosh 0, 
y = r sinh 0, the point lies on the hyperbola x 2 — y 2 = r 2 . The 
functions cosh and sinh are therefore termed “ hyperbolic 
functions.” 


35. Vectors. 

A scalar quantity is one having magnitude but no direction; 
examples are : numbers, temperatures, lengths, etc. A vector 
quantity is one having direction as well as magnitude ; displace¬ 
ments, velocities, accelerations and forces are all vectors. A 
scalar can be represented by the distance between two points ; 




Fig. 45.—(a) Representation of scalar quantities on a straight line; 
(6) Representation of vector quantities on a plane. 


I 
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if one point 0 is fixed, a set of scalars can be represented by the 
distances OA , OB, 00, etc., from 0 to points A, B, C, etc., on 
a straight line OP (Fig. 45a). A vector can be represented in 
magnitude and direction by the line joining two points. If one 
point 0 is fixed, a set of vectors in a plane XOY can be repre¬ 
sented by the lines OP, OQ, etc., joining 0 to points P, Q , etc., 
in the plane XOY. Let the co-ordinates of the point P be (x, y) 
(Fig. 456) and let N be the foot of the perpendicular from P to 

OX. A displacement from 0 to P is equivalent to a displacement 
from 0 to N together with a displacement from N to P ; if x 
represents a displacement of magnitude x in the direction OX, 
and y represents a displacement of magnitude y in the direction 

OY, then 

OP == x+y .... (35.1) 

As the point P is determined by the co-ordinates (#, y) the vector 
OP can be denoted by the set (x, y). 



It can be seen from Fig. 456 that the vector — (x, y) is in 
the same straight line POP' as the vector (x, y) but is in the 
opposite direction. Thus multiplication by — 1 is equivalent to 
a rotation through two right-angles or n radians. By (33.1a), 
e* 71 = cos = — 1 . . (35.2a) 

Hence rotation of a vector through an angle n radians is effected 
by multiplying it by e in . By convention, rotation of a vector 
through one right-angle or n/2 radians is said to be effected by 
multiplying it by e*" /2 . Now, from (33.1a), 

e i7t/2 = cog ^4-i.sm ? = i . . (35.26) 

and h&ice rotation through one right-angle is by convention 
supposed to be effected by multiplication by i. That this con- 
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vention is consistent can be easily shown. Beal scalar quantities 
being measured along OX, Fig. 46, the quantities i, 2 i, 3 i, etc., 
are to be measured off along OY. Thus 

OP = x+iy. 

Let x' +iy' — i(OP) = i(x+iy) — (— y+ix) = OP'. 

It is clear that OP' in Fig. 46 is at right-angles to OP. 

In order to distinguish between the number i ( = V— 1) and 
the operator i (turning the associated vector through an angle 
n/2 in the positive sense), the symbol j is frequently used for 
the versor-operator. (A rersw-operator is a turning-operator.) 
If the polar co-ordinates of P are (r, 6), then 

x+iy = r cos 0+i.r sin 6 = re® . . (35.3) 


This exponential form is sometimes useful in calculations. The 
following alternative forms are of frequent occurrence : 

1 — cos 0 +i sin 0 = e° 

— 1 = cos n +i sin n = e“ 


i — cos ^ +i sin ^ = e“ /2 
2 2 

i = cos ~sin ^ = e <3ll/2 

2 2 


. (35.4) 


Thus, for example, i(x+iy) = ire if> = re ^ 6+n / i 2 \ a form of state¬ 
ment which shows clearly the turning effect of the operator i. 


36. Integral formulae. 

The equations (29.1) are easily derived. 


r2n/m 1 r “l 2 n/<* 

I sin mcotdt = — — cos meat = 0 . (29.1c) 

Jo ml Jo 

By ordinary trigonometry, 

sin most . sin na>t = J cos (m — n)wt — £ cos (m+n)cot. 

■ If m and n are unequal, 


r2n/to 

* 

Jo 


sin meat . sin ncot dt 




2 n/a 
0 


n)(ot dt 

r2n/m 


= —— -Tsin (m — 

2 (to — »)L 


cos (to - 

— 11* cos (to +n)(ot dt 
Jo 

“l2n/°> 1 r ~\2n/tn 

i -^or (M+n H.v- 
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As both m and n are integers, both these terms are zero and the 
whole integral is zero. ...... (29.1a) 

•If, however, m — n, then as 

sin 2 meat = £ — \ cos 2mcot, 

C2n/<o f2 n/(o r2n/<o 

(1) dt — \ I cos 2 moot dt 

Jo 

.(29.1?) 

The other equations (29.1) are obtained in a similar manner. 


Jo 


sin 2 mcot dt = 




n 

ID 
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37. Newton’s second Law of Motion, units. 

Newton’s second Law of Motion states that the rate of change 
of momentum of a body is proportional to the impressed force. 
If the force in a certain direction is F, and the displacement in 
this direction is x, the resulting equation of motion is 

F = ^-(mx) . . . (37.1a) 

at 

where m is the mass of the body. If the mass remains constant, 

F = mx . . . . (37.16) 

Use of the equation (37.16) leads to two fundamentally different 
systems of units. In the scientific system the basic quantities 
are mass, length, and time (lbs., ft., secs.); (37.16) then serves to 
define unit force as being that force which produces in unit mass 
a unit acceleration. This unit force is called the poundal, and 
since the weight of unit mass produces in the mass an accelera¬ 
tion of 32-2 ft./sec. 2 , 1 lb.wt. = 32-2 poundals. In the engineers’ 
system, however, the basic quantities are force, length, and time 
(lbs. wt., ins., secs.) (37.16) then serves to define unit mass as 
being that mass in which unit force produces unit acceleration. 
This unit mass is called the slug, and 1 slug = 386 lbs., since the 
acceleration due to gravity is 386 ins./sec. 2 . 

(In aerodynamics a “ slug ” equal to 32-2 lbs. is used, the 
basic units being lbs., ft., secs., as it is more convenient to 
make measurements in feet when such quantities as Young’s 
modulus—usually given in lbs./in. 2 —are not involved.) 

Either system, the scientific or the engineers’, may be used, 
but care must be taken to express all the quantities consistently 
in one system. 

When a body rotates about an axis under the action of a 
torque T about that axis, the angular acceleration 6 is propor¬ 
tional to the torque in the same way as a linear acceleration is 
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proportional to the applied force. If the constant of propor¬ 
tionality is J, as yet undetermined, then 

T = IB . . . . (37.2) 

Consider a particle A of mass m at the end of a light arm OA 
of length r, rotating about 0 under the action of a force F at 
.A normal to OA (Fig. 47). The velocity of A is r& normal to OA, 



Fig. 47.—Derivation of I = mr 1 for a particle. 

and the acceleration in the direction of the force F is therefore 
r8 as r is constant. By Newton’s Second Law, 

F = mrS. 

The moment of the force F about O is Fr = T, and hence 

T = mr 2 0, 

and by comparison with (37.2) it is seen that 

I = mr 2 .(37.3a) 

This quantity I is termed the “ moment of inertia ” of the particle 
about O. 

A continuous rigid body can be regarded as composed of 
elemental masses dm { . The moment of inertia of the body about 
any axis is then defined as 

I = . . . (37.36) 

where ryis the distance from the mass dm,- to the axis of rotation, 
and the summation is extended over the whole body. The 
symbol J is frequently used for moments of inertia of solid 
bodies about an axis of rotatipn ; if the axis of rotation is also 
an axis of symmetry J is also termed the “ polar inertia.” 

The dimensions of moments of inertia are those of mass X 
length*; the scientific unit is therefore lbs.ft. 1 , while the engineers’ 
unit is slugs.ins. 1 ; in terms of the basic units a slug is lbs.ins. -1 
secs. 1 , being a force divided by an acceleration, and the engineers’ 
unit of inertia is therefore lbs.ins.sec8; 1 

A somewhat analogous quantity is, the “ second moment of 
. area ” ; this is for a cross-sectional al?ea what |he moment of 
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inertia is for a body; i.e. the second moment of area of any 
area about any axis is defined as 

I =,-Tr^da, .. .. (37.4) 

where da i is an elemental area at a distance r t from the axis, 
and the summation is extended over the whole area. The units 


of I are clearly ins. 4 or ft. 4 

The polar second moment of area of a circular cross-section 
is found to be \nR*, where R is the radius; this polar moment 
is with respect to an axis normal to the section and passing 
through the centre (Fig. 48a). 


(a) 



61 = 2Jtr.6ry 
I-Zemr’dr 
- iirff 


Fig. 48. — (o) Determination of polar 
second moment of area for circular 
cross-section; ( b) Rolling oscillation 

of a solid half-cylinder. 



Example. 


(b) 


o w 
6a 

6 



To illustrate the method of application of Newton’s Law in 
a complex problem, the natural frequency of small rolling oscilla¬ 
tions of a solid half-cylinder (Fig. 486) will be determined by 
ordinary dynamical methods, although the solution is most 
readily obtained by the energy method. Let the radius be r, 
the length l and the density p. The centre of area of the cross- 
section is on the line of symmetry OA, at a distance 00 = 4r/3?i 
from the geometric centre 0. (This result follows easily from 
one of the theorems of Pappus, which states that the volume 
of a solid of revolution is equal to the product of the generating 
area and the length of the path of the centre of this area; in 
this case the half-circle is rotated about the diameter PQ and 
generates a sphere. The generating area is and the length 

of the path of O is 2 n. (00 ); thus 2 n.(00) — %tr*, whence 

* 



120' FUNDAMENTALS OF VIBRATION STUDY 

The moment of inertia J a of the body about the axis 00' is 
pllo, where J 0 is the second moment of area of the cross-section 
about 0 and is thetefore Jjw 4 . Now J 0 = E(OW) 2 da, where W 
is the position of a typical elemental area da. The second 
moment of area about Q is easily shown to be I 0 — ( OQ)*.A , 
A being the area of the cross-section = Eda. For 

7 0 = E{OG*+GW* - 2.0G.GW.coa \QGW) da 

= 0Q*.A+I o - 2.0GEGW.coa \QGW .da 

where Ig = EGW 2 .da — second moment of area about G. The 
third term vanishes as G is the centre of area of the cross-section, 
and thus I a — I„ — (0G) 2 .A. (This is a perfectly general result, 
which holds true for any area.) 

Let OG = x ; the moment of inertia J Q of the body about 
an axis parallel to 00' and passing through G is given by: 

J Q = I G pl = (|7tr 4 — x 2 \nr 2 )pl 

= nr 2 pl{^ - = 7ir*pl($ - 8/9n 2 ). 

The solid is rolled to one side of the static position, without 
slipping on the supporting horizontal plane. Let B be the point 
of contact when the angular displacement is 8, and let the resolved 
reactions at B be R and S as shown. The vertical distance 
moved by G is x(l — cos 0) = \xB 2 , as 0 is small; the horizontal 
displacement of G is r.tan0 — x.sin0 = (r — x)6. Applying 
Newton’s Law for three displacements, i.e. for vertical and 
horizontal displacements of G and for the rotation about G, 

Ai 

S = m-r—[(r — x)6] = m(r — x)b (a) 

at “ 

R — mg = m d —^xB 2 ] = m^-pcdd] 

■— mxb 2 +mx0.6 ( b) 

Rx.aind+8(r — x.coa 0) = — J 0 B 
Rxd +S(r — x) = — Jgb (c) 

as 0 is small. 

Elimination of R and S leads to the equation: 

. m\gx0+x*6d*+x*6 t 9+(r - x)*8] - - Jjb . (37.6) 



and 

i.e. 
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Putting 0 = a sin cot, and considering the instant when cos cot = 0, 
sin ail = 1, the equation is obtained: 

m(gxa — x 2 a*co 2 ) — m(r — x)*aco* = +J c flco i , 

and omitting terms containing powers of a higher than the first, 

mgx — m(r — x) s co a = JqCo 8 

whence a> 8 = ■ _ . . . (37.7) 

9ji — 16 v ’ 

by substituting for m and J Q ; the natural frequency co/2n is 
thus obtained. 

38. Bending of Beams. 

Fig. 49 represents a short section of a beam which is deformed 
by bending; circular deformation is assumed so that all plane 
cross-sections remain plane and undistorted and all pass through 
a line normal to the plane of the longitudinal section Fig. 49a, 



Fig. 49.—Bending of beams : (a) longitudinal section ; ( b ) perspective view; 

(c) transverse section. 

the intersection of this line with the plane of Fig. 49a being O. 
Considering the beam to be composed of fibres originally parallel 
to AB and DC in the undeformed condition, these fibres are 
deformed into circular arcs, and those lying in the plane of the 
distgram are concentric at 0. The fibres nearer 0 are compressed 
and those further from 0 are extended. There is a plane EFRS 
such that all the fibres in it are merely bent, without extension 
or compression. This plane is termed the “ neutral plane,” and 
the line of intersection of the neutral plane with a cross-section 
is termed the “ neutral axis.” 
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Let E = Young’s modulus for the material 

p = radius of curvature of deformed neutral plane = OE 
h = distance of slice PQ from neutral plane 
dh — thickness of slice 
e = extension per unit length of fibre PQ 
l = undeformed length of PQ = EF 
M = bending moment at BCKL 
8 = stress on area QQ'T'T 
t = width of area QQ'T'T = QT 


By the properties of similar figures, 

arc PQ _ OP _ p+h _ A 
arc EF OE p — p 


but 

arc PQ 

= l(1 + €) = l+e . 


* arc EF 

1 

Hence 


h 

e = 



P 

Now 

, , tensile stress 

e = tensile strain = „- 

E 



Eh 

hence 


s = — 


P 


8 

E 


(38.1) 

(38.2) 


(38.3) 


The tensile force on the area QQ'T'T is s(area) = stdh and this 
force has a moment 6M about FR, where 


6M — sthdh. 


The total moment for the whole area BCKL is the bending 
moment M ; thus 

M = Eslhdh 

or, by (38.3), 

M = -Ethm . . . . (38.4) 

P 

JPT 

or M — — where / = Eth*dh . . (38.5) 

P 

(see Section 37) 

The resultant tensile force on the area BCKL is zero, henoe 

Estdh — 0 

i.e. -EhtSh - 0. 

P 
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This last equation expresses the fact that the neutral axis FR 
passes through the centre of area of the cross-section BCKL. 

If the distances along the undeformed beam are denoted by 
x, and y is the deflection of the beam from its undeformed position, 
the radius of curvature is given by 

d*y . 

1 _ dx* 

; ' 

For small displacements for which is small, this formula 
, 1 dhi 

reduces to - = and (38.5) becomes 

’ M = El d p- . . . . (38.6) 

CUE 

The deflection curve of a loaded beam is easily obtained. 
Consider, for example, a light loaded cantilever, the loading 
being concentrated at the free end (Fig s 50). At a distance x 



from the built-in end the bending moment is given by 
M = P(l - x) 

i.e. Big = P(l - x) 

hence eA = Plx - A'+ A 

dx '2 



Px 8 
6 


-\-Ax -f- B. 


At the built-in end (x = 0) the displacement and slope are zero; 
hence A — B 0, 



Pa* 

T• 
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If A is the end-deflection at x — l, it is found that 



The linear stiffness of the beam to an end-load P is thus 



(38.7) 


39. Torsion of a circular shaft. 

The shear modulus 0 of a material is defined thus : let ABCD 
be a small cube of the material (Fig. 51a), and while the face 




Fig. 51.—(a) Shearing deformation of small cube, (6) Torsion of circular shaft. 

AB is fixed, let the face CD be subjected to a shearing stress q, 
so that the cube deforms into a rhomboid prism ABC'D'. Let 
the angle CAC' be ip, then 

G = - ... . (39.1) 

V 

Now consider a short section 6x of a uniform circular shaft 
(Fig. 516) of radius r. Under torsion the generator PQ becomes 
PQ' ; let the angle QOQ' be 66, Q being the centre of the cross- 
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section through Q and Q'. Imagining the small cube of Fig. 51a 
to be at P ,.the shear strain ip is given by 

* M 

ton v =rg- 


i.e. ip = r—, as ip is small. 

dx 

Let the area of that side of the cube which is included in the 
cross-section of the shaft at P be da, on which the shear stress 
is q — Qip (by 39.1). The shear force on da is then qda, and by 
the above formula this is equal to 

add. * 

Gr-~da ; 
ox 

the moment of this force about the axis AB is 

6T = OrAa, 
dx 

and the total torque transmitted by the section is 

T = ZGr*daf = Q ^-. . . (39.2) 

dx dx v ' 


where I is the second moment of area of cross-section (see 
Section 37), as the summation is extended over the whole 
cross-section. 
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FREQUENCY EQUATIONS AND NUMERICAL 
SOLUTIONS 


40. Determinants. * 


The equations of motion for a complicated vibrating system 
contain a number of variables representing linear and angular 
displacements; elimination of these variables in order to obtain 
a frequency equation similar to (16.5) is often best performed 
by means of a determinant. The mathematical properties of 
determinants are fully treated in standard intermediate text¬ 
books on Algebra; sufficient will be restated here to enable the 
method to be explained. 

Hie set of simultaneous linear equations in three variables 


aiX+b{y-\-c t z = O' 

a&+b$+c# = 0> . . . (40.1a) 

a&+b$+c& = OJ 

has a solution x = y = z — 0, and may have another solution 
in which at least one of the variables x,y,z is not zero. Suppose 
z is not zero; the set can then be written: 


= 0 

a 1 X+6,F -|-c, — 0 
a,X-|*6,F-l-c, — 0 
where X = x/z, Y = y/z 


(40.16) 


Solving the second and third equations of this transformed set 
for X, Y, 

X - Y . 1 


6|C| — cjb t a,c, ~~ CgUg a,6, 6,a, 

X - Y 1 


6, c, 


a, c, 


a, 

b t 

b,. c. 


a* c * 


a, 

b t 


where 

(40.2) defines a 


6j c, 
6, c, 


— 6,c, — c,6j, etc. 


(40.2) 


second-order determinant, represented symbolic* 
126 
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ally by the left-hand side of the equation. Substitution of the 
values of X and Y in the first equation of (40.16) leads to 


o, 


6 t c, 
6 * c* 


- 6 , 


a * c * 


+®l 


a, 6 
®* 


= 0 


and this last equation is expressed.in the form: 


a x 6 X c x 

Qj | 63 Cf 

a3 63 C3 


= 0 


(40.3a) 


(40.36) 


This arrangement of letters in three rows and three columns is 
called a third-order determinant, and is defined to be identical 
with the left-hand side of (40.3a). (40.36) expresses the condition 
that the set of equations (40.1a) should have a solution inx,y,z 
other than x = y — z — 0. 

The rule for expanding a three-row determinant, as indicated 
in (40.3a), can be stated as follows : “ Associate with each 
element a x , 6 X , c x of the top row the signs +, —, + alternately, 
always giving the left-hand element the positive sign ; multiply 
each such element, with its associated sign, by the determinant 
formed by omitting the row and column containing the element, 
and sum the products so obtained.” Thus the determinant 


formed by omitting the row and column containing 6 X is 


®* c, 
a, c. 


and the sign to be associated with 6 X is negative ; similarly for 


a x and c x . 


It is convenient to express a determinant concisely by writing 
down the elements of the diagonal which runs from the top left 
comer to the bottom right comer; thus the determinant in 
(40.36) is represented by |a x 6, c 3 | and the three-row determinant 
is defined: 


where 


|a x 6, 


Cf|= ®i|6* Cs|- 6 x |a. c,|+c x |a, 6,| 
| 6 t Cjj = 63 C 3 C 363 , etc. 


(40.4) 


The following rules of manipulation are proved in Algebra 
textbooks, and may easily be demonstrated in particular cases 
by performing them on numerical determinants: 

(i) The value of a determinant is unaltered if columns and 
rows are interchanged. Thus 
7 5 3 
0 1 0 
5 8 2 


= 7(2 - 0) - 5(0)+3(0 - 5) = - 1 
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and 705 

5 18 = 7(2 - 0) - 0+5(0 - 3) = - 1. 

_ 3 0 2 

. (ii) If two rows or columns are interchanged the value of the 
determinant is multiplied by — 1. Thus 

0 10 

7 5 3 =- 1(14 - 15) = 1. 

5 8 2 

. (iii) If two rows or two columns are identical, the determinant 
has zero value ; this follows from (ii). Thus 

0 1.0 

7 5 3 = 1(21 - 21) = 0. 

7 5 3 

(iv) If all the elements of a row or column are multiplied by 
the same constant, the value of the determinant is multiplied 
by that constant; thus 


14 10 6 


7 5 3 


7 5 3 

0 10 

= 2 x 

0 1 0 

= 2 x 3 x 

0 1 0 

15 24 6 


15 24 6 


5 8 2 


(v) If the elements of a row or column are each expressed 
as the sum of two numbers, the determinant can be expressed 
as the sum of two determinants, the other rows or columns 
remaining unaltered. Thus 


7 

5 

3 


5+2 

8 

— 

3 

2 + 1 




0 

1 

0 

= 

0 


i 


0 





5 

8 

2 


5 


8 


2 








5 

8 2 


2 

- 

3 

1 

0 

1 

0 



= 

= 0 

1 0 

!+ 

0 


1 

0 

= - 2 - 

- 3 

1 




5 

8 2 


5 


8 

2 

5 

8 

2 


■ (Vi) If to the elements of any row, or column, are added the 
same multiples of the corresponding elements in any other row, 
or column, the value of the determinant is unaltered. Thus 

117 21 71 17+10 5+16 3+41 


0 10 = 0 

1 

0 


5 8 2 

8 

2 


7 5 3 


5 8 2 


7 6 3 

= 010 

+2 

0 1 0 

= 

0 1 0 

5 8 2 


5 8 2 


5 8 2 
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These rules are very useful for determining the value of a numerical 
determinant. For example, let 

1 2 3 
A = 4 6 5 
7 8 9 

Subtracting the top row from the bottom row (rule vi) 


A = 

1 2 3 
4 6 5 

= 6 

1 2 3 
4 6 5 

= - 6 

1 1 1 
4 6 5 


6 6 6 


1 1 1 

(rule iv 

i 

1 2 3 

(rule ii) 


Subtracting the left-hand column from the other two, (rule vi), 
1 0 0 

A = - 6 4 2 1 = - 6(4’- 1 ) = - 18. 

112 

All the results stated in this section are capable of extension 
to determinants of higher order. Thus, for example, the definition 
of a fourth-order determinant is : 

= Ui | & jC 3 ^ 4 1 — — dijUghgC 4 j . (40.5) 

and the equation la^^d^ = 0 is the condition that the set of 
four simultaneous linear equations 

'a 1 x+b 1 y+c 1 z+d 1 w = 0 
ja^c+b^y+c^+diW = 0 
\a&+b&+c i z+d a w = 0 

jUtX+bty+CtZ+diW = 0 

has a solution in x, y, z, w other than x=y = z — w = 0 . 

41. Frequency equations. 

Suppose that the equations of free undamped motion of a 
vibrating body having three degrees of freedom are: 

a l X+ 6 1 r+c 1 2 +mJ£ = O' 
a t X-\-b t Y+c&+mY = 0 > . . (41.1a) 

d a X-\-b a Y -\-ctfi-\-Jz = 0 ^ 

where m is the mass of the body, J its moment, of inertia with 
respect to the rotation z, and X and Y represent linear displace¬ 
ments. For sinusoidal motion X = A sin cot, etc., X = — eo a X, 
etc., and the equations of motion can be written as 

(oti — moj a )X-\-b\Y -f-CiZ = 0 ") 

atX+( 6 » — mco 2 )Y+c& = 0 > . . (41.16) 

a,X+ b t Y +(c, — J(o*)z = OJ 


K 
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and there is a solution other than X = F = z= 0 if 

(a x — nuo 2 ) bx c x 

" a , (6, — ma> s ) c» =0 

o 3 6* (c s — Jen 2 ) 

Expansion of this determinant by (40.4) gives an equation of 
the form 

Pa> 6 +<2a> 4 +-Rtt> a +$ = 0 . . (41.2) 

where P, Q, R, S are expressions involving the constants m, J, o„ 
etc. 

Similarly, the equations of motion (free and undamped) of 
a body having six degrees of freedom (Fig. 23, Chapter III) are 
of the form 

diX+bxY+c 1 Z+d 1 x+e 1 y+f l z+mX = 0, 
etc. etc. 

The resulting frequency equation analogous to (41.2) is of degree 6 
in co 2 . As noted in Chapter III, section 16, however, symmetry 
in the spring-mounting scheme may result in certain modes of 
vibration being uncoupled. In a standard type of in-line aero¬ 
engine mounting, for example, if X, Y, Z represent displacements 
vertically, transversely and fore-and-aft respectively, and x, y, z 
are rotations about corresponding axes (see Fig. 23), the following 
pairs of displacements are uncoupled : 

X and x X and Y 

Y „ y X „ Z 

Z „ z Y „ Z 

X „ z x „ y 

Z „ x y „ z 

The equations of motion then reduce to the set: 


a x X 


+e t y 

-(-jW'JSl = 0 

b»Y 

+dfX 


+/jZ +mY = 0 


c 3 Z 

+e s y 

+mZ = 0 

b t Y 



+/t* = 0 

a s X 

+c f Z 

+e t y 

’i’JyV = 0 

b t Y 

+d,x 


+/«Z +Jji = o 


where J x , J v , J z are the moments of inertia of the body with 
respect to rotations x, y, z respectively. The first, third and 
fifth equations contain only X, Z and y, and the second, fourth 
and sixth equations contain only Y, x and z ; the two resulting 
determinantal equations are therefore : 



(a t — mco 2 ) 
0 
(Is 


APPENDIX m 


131 


0 e x 

(fit — mco 2 ) c 8 =0 
Cs (e 5 — J y co*) 

(b t — mca 2 ) <£, /, 

6 4 (d 4 — J x (o 2 ) f t =0 

bt dt (ft - J z (o 2 ) 

and these lead to two cubic equations in co 2 ; the six natural 
frequencies are thus determined. 

42. Numerical solutions. 

Suppose that a cubic equation in eu 2 is obtained by the method 
of Section 41. It is required to find the three natural frequencies. 
Let the equation be, for example, 

w # - 6-51 x 10 6 o> 4 +13-19 x 10 12 cd 2 - 7-99 X 10 18 = 0 (42.1) 

Put x = co* x 10 _ * 

then X s - 6-61z 2 +13-19* - 7-99 = 0 . . (42.2) 

Let the left-hand side of (42.2) be represented by y. 

When x — 0, y = — 7-99'| 

1 — 0-31 > ; there is therefore a root between 

2 0-35J 

1 and 2 . Put x = 1 +A, and substitute in (42.2), neglecting 
powers of h higher than the first. Thus : 

1+3 h - 6-51(1 +2h) + 13-19(1 +h) - 7-99 = 0 
i.e. — 0-31+3-17A = 0, h = 0-1 approximately, 

and * = 1 +A = 1-1 approximately. Now put x = 1 - 1 +4 and 
solve for k in a similar manner ; it is found that k is 0-01 approxi¬ 
mately, so that a nearer approximation to a; is Ml. The process 
is one of applying a series of decreasing corrections to an approxi¬ 
mate value of the root found by inspection or graphical methods ; 
the next correction in the series 0 - 1 , 0 - 01 , ... is found to be 
— 0 - 002 , so that the value 1-11 is correct to three significant 
figures. The equation is now divided through by (* — 1-11), 
i.e. the left-hand side is factorised with (x — 1 - 11 ) as one factor : 

(* - 1-11 )(* 2 - 5-40*+7-20) = 0, 
and the other two roots are found by solving the quadratic 
equation 

* 2 - 6-40*+7-20 = (* - 3-0)(* - 2-4) = 0 . 
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Hence the three roots in x are : I'll, 2*4, 3*0, and the correspond¬ 
ing values of to 2 are : 

<o 2 = 1 11 x 10«, 2-4 x 10 6 , 3-0 X 10« 
and F = co/2n = 167, 247, 275 O.P.S. . . . (42.3) 

The method is often stated in another form. Let x = a x be 
an approximate solution to the equation f(x) = 0. Let the 
function be represented in the neighbourhood of x — a x by the 
curve of Fig. 52; Q X R = f(a x ) and if the approximate root o, 



Fio. 52.—Approximation to solution of a numerical equation. 


is sufficiently close to the true value x, a closer approximation 
is found by drawing the tangent RQ t at R and finding the point 
Q t in which the tangent cuts the x-axis. Now Q X R — Q x QJ'(a x ) 

where /'(a,) means the value of ^/(x) when x = a x . Thus 


a t — Oj — QiQi — ®i 


/(q i) 

/'(«i) 


(42.4) 


Putting a 2 for a u a still closer approximation x = a, is found 
in a similar manner. As an illustration, consider this method 
applied to the equation (42.2). 


Hence 


/(x) = x 3 - 6-51x 2 +13-19x - 7-99 
/'(x) = 3x 2 - 13 02x+13-19 
e a x = 1, f(a x ) = — 0-31 and /'(o x ) = 3-17 



- 0-31 
3-17 


= 1-1 approx., as before. 


It is evident that the two methods are equivalent; the second 
form of the procedure is more easily applied to higher order 
equations. It is necessary in any case to obtain as a first 
approximation a value which is not too widely different from 
the. true solution, particularly if the curvature of the graph 
y =/(x) is large in the neighbourhood of the root. 
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43. Normalising of numerical equations. 

It frequently happens in the analysis of an engineering 
vibration problem that a large number of fairly simple equations 
(quadratics and cubics) require to be solved. In such circum¬ 
stances much time can be saved by normalising the equations 
and making use of special tables. The process, which was 
known to the Babylonians circa 2000 b.c. for certain types of 
cubic equations, will be illustrated here by applying it to the 
general quadratic equation, 


ax 2 +bx+c = 0 


Let x = yd, so that ay 2 d 2 +byd+c = 0, 

ie - ),+ ^ + ^ = °’ 


(43.1) 


and choose d so that the coefficient of y is unity, i.e. d — b/a. 
Then 


y 2 +y+R = 01 
where R = ac/b 2 J 

The roots in y are 


(43.2) 


and 


y = - i(l ± Vl - 45)] 
6 


x = -y 
a 


(43.3) 


Special tables can be constructed giving the two values of y 
corresponding to any value of R ; for real roots it is only neces¬ 
sary to consider values of R < 0-25, and the range of negative 
values included will depend to some extent on the nature of the 
problem which gives rise to the equation. If the roots in * are 
essentially positive, as will be the case if a; is proportional to a> 4 , 
a)/2jt being a natural frequency, then R must be positive and 
is restricted to the range 0 <5 <0-25. Tables giving the 
values of y corresponding to this range of values of R are easily 
prepared. 

The method of normalising equations is of no great use unless 
a large number of equations of the same type have to be solved, 
but the time saved in such cases is considerable, apart from the 
fact that the results are easily tabulated together with the 
incidental working. 
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PRACTICAL TREATMENT OF VIBRATION 


Although practical details have in general been excluded from 
the body of the text, as falling outside the scope of the book, 
this appendix describes the treatment of simple vibrations as an 
illustration of the method of applying the fundamental principles. 

The basic theory of anti-vibration mounting and insulation 
is outlined, and its application discussed. Various forms of 
resilient mountings are described. 


44. Theory of anti-vibration mounting. 

For practical purposes it is often sufficient to confine attention 
to the response of a single-degree-of-freedom system without 
damping, under the following conditions: 

(а) Flexible suspension fixed to a rigid immovable structure, 

a sinusoidal force being applied to the suspended mass 
(Fig. 9, page 25). 

(б) Flexible suspension fixed to a movable base, to which is 

imparted a sinusoidal vibration (Fig. 14, page 37). 

Case (a). 

From formula (9.3) on page 26 it is seen that, disregarding 
the effects of damping forces, the amplitude of vibration of the 
mass is 


b = 


Pj _ 

m(co 2 — p*) 


(44.1) 


where the symbol | | denotes that the numerical value of the 
fraction is to be taken without sign, 

P 0 = amplitude of disturbing force applied to mass, 
m = suspended mass, 

a)/2n = natural frequency of mass on suspension, 
p/2n = frequency of disturbing force. 

The equation (44.1) can be interpreted in two different ways. 
First, if the variation of the response amplitude b with changes 
in the forcing frequency is considered, the relation is found: 


b_ 

C, 11 - X* 

X =p/(o 
Ci — PJnuo* — PJk 
134 


where 


> 


. (44.2a) 
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For very small values of the frequency-ratio X the ratio b/C 1 
approaches unit value, and in this condition the amplitude is 
equal in magnitude to the deflection that would be produced by 
a static force P 0 acting upon a spring of stiffness k. For very 
large values of the frequency-ratio, the ratio b/C 1 tends to the 
value 1 / X a , so that the amplitude b approaches the value P 0 /mp a , 
which is equal to the amplitude of response of an unsupported 
mass to the same vibratory force. At intermediate frequencies 
the response amplitude displays the resonance phenomenon; 
some values of the ratio b/C v for various value of the frequency- 
ratio X, are given in Table XII on page 146. (In using this 
table it should be remembered that damping has been neglected 
in the present treatment; the effect of damping forces will always 
be to decrease the amplitudes given by the figures in the table, 
and so long as the forces are small their effect is confined to 
frequency-ratios very nearly equal to unity.) 

Secondly, if the forcing frequency remains constant and the 
stiffness of the suspension is changed, it is often convenient to 
express equation (44.1) in the form: 

b _ 1 

C 2 Y - 1 

where Y = k/mp 2 = 1/X a 

and C t = P 0 /mp 2 

For very small values of the ratio Y (i.e. for a very flexible 
suspension or a high-frequency force) the ratio b/C t approaches 
unit value, and in this condition the amplitude is equal to that 
produced by the same force acting on an unsupported mass. On 
the other hand, if the suspension is very stiff the ratio Y is 
large, and in this condition the amplitude approaches the value 
P 0 /k, which is equal in magnitude to the deflection which would 
be produced by a static force P 0 acting on a spring of stiffness k. 

It must be noted that although the provision of a very flexible 
mounting (corresponding to a small value of Y, and a large value 
of X) is frequently desirable for the purpose of avoiding resonance 
with the disturbing force and with any possible harmonics (at 
higher frequencies) if this force is not a pure sine-wave, the 
advantage is not as great as at first sight it appears to be. The 
ratio b/C x is then very small, but C x = PJk is large. It is 
often stated that a good measure of vibration prevention is 
obtained if the frequency-ratio is made not less than 3:1; but 
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there is. no effective standard for comparison. Clearly, a rigid 
mounting would entirely prevent vibration, although it has 
disadvantages on account of its high degree of transmission of 
vibration through to the supporting structure. If a flexible 
mounting is fitted, then in order to avoid resonance with possible 
harmonics it must be more flexible than a mounting which would 
give a condition of resonance with the fundamental component 
of the force ; and equation (44.26) shows that the amplitude of 
response cannot then be less than C t . The fact that this value 
is decreased by increasing the mass gives the reason for the 
common practice of mounting machinery on massive concrete 
plinths which are themselves supported on flexible members. 

Example. 

A horizontal shaft carrying an unbalance of 1 oz. at 19J inches 
radius and running at 1,800 r.p.m. is mounted on a frame which 
is supported on four flexible mountings, between which the total 
weight of 760 lbs. is evenly distributed. If the rotor is at the 
mass centre of the supported body and the mountings each have 
a vertical stiffness of 12,000 lbs. per inch, calculate the amplitude 
of vertical vibration. Calculate also the amplitude if (i) the 
running speed is increased to 2,100 r.p.m., and if (ii) the mountings 
each have a stiffness of 6,000 lbs. per inch, the speed being 
1,800 r.p.m. (Assume that the body can only move vertically, 
with no “rocking.”) 

Here m = 750/386 = 1-94 lbs.ins. _1 secs. 2 For 1,800 R.P.M., 
p 2 — 3-56 x 10 4 secs. -8 (from Table XI, page 145). 

P„ = 19-25p 2 /(16 X 386) = 111 lbs., 
and with the original mountings, 

h = 4 x 12,000 = 4-8 x 10 4 Ibs.ins., 
o>* = k/m = 2-47 x 10 4 secs. -2 The natural frequency of ver¬ 
tical vibration is therefore 1,500 c.p.m. (from Table XI). Then 
X = 1,800/1,500 = 1-2, b/C t = 2-27 (from Table XII) and 
C t = P 0 /k = 2*31 x 10 -3 inches. Hence 6 = 0-0053 inch. 

For a speed of 2,100 r.p.m., P 0 = 151 lbs., X = 1-4, 
b/C x = 1-04 and 6 = 0-0033 inch. 

For 1,800 r.p.m. and 1 = 4 x 6,000 = 2-4 x 10 4 Ibs.ins. -1 , 
Y — 0-348, b/C % = 1-53 and C t = 1-6 x 10“ 3 inches. Hence 
6 = 0-0025 inch. 

Case (b). 

For the system of Fig. 14, with c = 0 (no damping) and a 
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vibratory movement of amplitude a and frequency p/2n c.p.s. 
imparted to the base A, the equation of motion (i) on page 37 
yields the result: 


where 

and 


b 


k 


1 

1 

a 


k — mp 2 


1 -X 2 



x = a sin pt . 

y = b sin pt 
X = p/(o 


(44.3) 


The numerical value of the ratio b/a indicates the ratio of the 
amplitude of response of the supported mass to the amplitude 
of movement of the base A. Comparison of equations (44.3) and 
(44.2a) shows that the values given in Table XII are applicable 
to this case also. For very small values of X (i.e. for very low 
frequencies of excitation or for a very stiff mounting) the response 
of the mass is equal to the applied vibration. For very large 
values of X (i.e. for very high frequencies or for a very flexible 
mounting) the response is negligible. It is thus seen that in order 
to provide a good measure of insulation the mounting should 
be made very flexible. If the frequency-ratio is greater than 
V2 : 1, the response will be less than that with a rigid mounting ; 
for example, if X is greater than 3:1, the response will be reduced 
to less than an eighth of that with a rigid mounting. 


Example. 

A galvanometer is to be mounted on a table which is subject 
to vertical vibrations of amplitude 0-010 inch at a frequency of 
3,000 c.p.m. The galvanometer, which weighs 8 lbs. complete 
with a suitable base-plate, must not vibrate with an amplitude 
greater than 0-0005 inch. Relying solely upon a flexible suspen¬ 
sion to achieve this measure of isolation, how flexible should it be ? 

From Table XII it is evident that the value of X must be 
made somewhat greater than 4-5 to obtain a value 0-05 for the 
ratio | b/a |. Since X is greater than unity, a and b are of opposite 
signs, and the critical value of X can be found from the equation 


Hence X = 4-6 approximately, and the natural frequency of the 
suspension must not be greater than 3,000/4-6 or about 650 c.p.m. 
From Table X it is found that the gravity defleotion corresponding 
to this frequency is 0-0833 inch. If the load is evenly distributed 
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over four springs, each must deflect at least 0*0833 inch under a 
load of 2 lbs., so that the stiffness of each must not be greater 
than 24 lbs. per inch. 


45. Theory of vibration isolation (continued). 

One type of isolation has been described in the preceding 
section ; this is the isolation of suspended bodies from vibrations 
of the supporting structure. Another type is of importance, 
viz. the prevention of the transmission of vibratory forces from 
a vibrating body through its suspension to the supporting structure 
or to the ground. 

Considering the system of Fig. 9, page 25, but with no 
damping, the force transmitted by the spring to the support is 
always equal to the product of the spring stiffness k and the 
displacement x of the mass. Consequently, if P t is the amplitude 
of the force transmitted, 




k 

k — mp 2 


1 _ 7 

1-X 2 7-1 

where X = p/m, 7 = k/mp 2 = l/X 2 


. (45.1) 


This ratio R is a useful quantity for consideration, as it expresses 
in a convenient form the effectiveness of an isolating installation. 
If R is unity the force transmitted to the support is the same as 
in the case of a rigid suspension ; if R is less than unity, a degree 
of isolation has been achieved ; while if R is greater than unity, 
a greater vibratory force is transmitted than in the rigid case. 
The variation of R with changes of 7 is plotted in Fig. 53c, and 
it can be seen that in order to obtain a degree of isolation 7 
must have a value less than \ ; i.e. the suspension must have a 
stiffness less than half that value which would correspond to a 
state of resonance with the given applied force. 


46. Application of theory. 

From the foregoing sections it will be apparent that three 
graphs are useful for the appreciation of the properties of flexibly- 
mounted single-degree-of-freedom systems. These are the graphs 
of the numerical values of the quantities : 




(46.1) 
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against X, Y and Y respectively, where 

g _ p _ frequency of disturbing force 
co natural frequency of system 

Y _ 1 _ & _ actual stiffness of suspension 
X 2 mp 2 /stiffness corresponding to 

^resonance with given force 

These quantities are graphed in Fig. 53, the scales of which are 
lettered to correspond to (46.1). 

(a) shows the variation of the amplitude of response of the 
mass to a vibratory force applied directly to the mass, assuming 
that the only quantity altered is the frequency of the applied 
force. (Note .—If the frequency-ratio X is altered by varying the 
spring stiffness, the appropriate graph is (6), as such a variation 
alters the value of C x in (44.2a).) It shows also the variation 



of the amplitude of force transmitted by the spring to the support, 
or of the amplitude of response of the mass when the support is 
subjected to a vibration, for changes in the frequency-ratio X 
brought about in any manner. The unit value, approached at 
very low values of X, is in the first case the magnitude of deflection 
which would be produced in the spring by a static force equal 
in magnitude to the applied dynamic force, and in the other cases 
it is the amplitude of the applied force or vibration itself, which 
in the case of a rigid suspension is transmitted through the 
suspension unaltered. 
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( b) shows how the amplitude of response to a directly applied 
vibratory force depends upon the stiffness of the suspension. 
The unit value, approached at very low values of Y, is the 
amplitude of response of the mass to the same force if it were 
unsupported. 

(c) shows the variation of the force transmitted through the 
suspension when a vibratory force is applied directly to the mass, 
and also of the response of the mass when the support is subjected 
to a vibration, against changes in the stiffness of the suspension. 

From a study of these graphs it will be apparent that for the 
purposes of isolation, i.e. preventing the transmission of vibratory 
forces from the mass (to which they are directly applied) through 
the suspension to the support, or preventing the transmission of 
vibrations in the support through the suspension to the suspended 
mass, it is sufficient to provide a suspension which is flexible 
enough to give a value not less than 3 to the frequency-ratio X ; 
i.e. so that the natural frequency of the mass on its suspension 
is not more than one-third of the frequency of the disturbing 
force. Naturally, the more flexible the suspension is, the greater 
will be the degree of isolation. A very flexible mounting scheme, 
however, has the disadvantage that the mounted mass is not very 
precisely located, and has a good deal of freedom in small but 
possibly undesirable displacements. 

When it is desired to prevent vibration in the mass which is 
directly acted upon by the disturbing force, without at the same 
time causing large or relatively large vibratory forces to be trans¬ 
mitted to the supporting structure, it is customary to add sub¬ 
stantially to the mass of the supported body, by means of a large 
concrete block, and then to mount the whole assembly on flexible 
units. The effect of the mass increase is to reduce the unit 
value C x in (44.2 a). 

The effect of the flexible units is to give a large value to X, 
so making both the ratios b/C 1 in (44.2a) and b/a in (44.3) small; 
at the same time the large mass ensures that the unit value C 
in (44.2a) is small. 

In the great majority of cases it is not possible to alter the 
forcing frequency. The most common causes of vibration in 
engineering practice are unbalance of rotors and gas forces in 
internal-combustion engines and compressors ; in these cases the 
frequency of the disturbing force is dependent on the running 
speed of the machinery, and it is only seldom that any radical 
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change can be made in this speed. Consequently, the desired 
large value of X can only be obtained by the use of fairly flexible 
mounting units, various types of which are described in the 
following section. 

In assessing the vibration characteristics of any actual system 
it is important to realise that the idealised types of system which 
are necessarily studied in mechanics are over-simplified in order 
to make a mathematical analysis possible ; a “ rigid suspension ” 
does not exist. Even apparently rigid structures may be appreci¬ 
ably flexible so far as vibrations are concerned. Thus if a 
“ rigidly mounted ” machine is noticeably vibrating as a whole, 
there must be a certain amount of flexibility in the supporting 
structure. The degree of flexibility can be assessed quite accur¬ 
ately enough in most cases by observing the natural frequency 
and substituting in the formula (3.5) on page 8. Any additional 
flexible unit is to be regarded as being in series with this flexibility 
in the structure. 

Fortunately, for a variety of reasons the majority of vibrating 
systems which cause trouble in general engineering are for 
practical purposes equivalent to single-degree-of-freedom systems, 
and there are usually not more than two possible frequencies of 
disturbing forces. Consequently, the theoretical considerations 
outlined above enable a successful solution to be found for most 
problems. In certain specialised fields, however, and particularly 
in aircraft engineering, a large number of natural frequencies may 
be significant and at the same time there may be many different 
disturbing frequencies, so that the treatment of vibration problems 
becomes very complicated. The author has given elsewhere 
(reference 16 in the Bibliography on page 151) a general outline 
of the vibrations in aircraft. In these more complicated cases 
a tremendous amount of actual experimental work is usually 
necessary, mathematical analysis proving to be of little use owing 
to the large number of unknown quantities, or to the complexity 
of the calculations. Nevertheless, an appreciation of the broad 
principles underlying the treatment of simple cases is very 
helpfrd as guidance in all problems. 

47. Flexible mounting units. 

The most commonly employed type of flexible unit for anti¬ 
vibration and insulation purposes embodies rubber as the flexing 
medium. The rubber is usually stressed in shear or in oompres- 
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sion; its strength in tension is too small to allow of tensile 
stressing in general, although in special cases this is quite possible 
and sometimes convenient. Rubber units can be divided into 
two classes: those in which the metal parts are bonded to the 
rubber, and those in which the rubber is pre-compressed between 
the metal parts, as in the Silentbloc flexible bearings and anti¬ 
vibration mountings. 

Bonded units are available in a wide variety of shapes and 
sizes. The design of such units calls for specialized knowledge 
based on experience; clean and simple outlines, and avoidance 
of stress concentration in the rubber, are among the essentials. 
Several firms specialise in the manufacture of such mountings. 

Pre-compressed units can take advantage of extremely rugged 
construction, to withstand heavy static and dynamic loads. The 
initial state of compression enables heavy loads to be carried by 
mountings of small volume; the standard sizes range from a 
mounting to take a static load of 20 lbs., with a length of inch 
and a diameter of § inch, to one capable of carrying 7,000 lbs., 
with a length of 9£ inches and a diameter of 9 inches. 

Metal springs are also very often used for anti-vibration and 
insulation purposes. A mounting scheme for a 680 b.h.p. 
generator set is a typical example. The load is taken, via 
girders projecting from the concrete foundation block (attached 
to the set) and vertical rods, by semi-elliptic springs, the outer 
ends of which bear on plinths forming extensions of the base raft. 
Various other forms of spring supports or “ anti-vibrators ” are 
supplied by the same firm ; loads from 28 lbs. to 5 tons can be 
taken on each unit, depending upon the size. 

In another method of foundation isolation, the isolating 
medium consists of specially selected and treated virgin cork, 
cut into strips and arranged in a packed formation inside iron 
frames. The resulting units, which are 2£ inches thick and up 
to 18 square feet in area, are then impregnated with a preservative 
to maintain resilience and prevent deterioration. The base raft 
is first laid, and the whole area covered with the units, over which 
thin felt is placed to seal the surface for the pouring of the main 
concrete block. An air space is maintained round the sides of 
the block, but in the case of belt-driven installations the space 
on the side of the belt-pull may be filled with the cork units. 

In ail types of installation, whether the insulation is achieved 
with rubber, metal springs or cork, it is important to prevent the 
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insulation being “ short-circuited ” by holding-down bolts, rigid 
pipe-connections or other solid paths for the transmission of 
vibration forces. When there is a massive concrete block, the 
mass is usually sufficient to make holding-down bolts unnecessary, 
and in the case of rubber mountings the units act themselves 
as holding-down devices; when it is essential to pass a bolt 
through an insulating unit, short-circuiting may be avoided by 
providing a flexible washer under the bolt-head or nut and a 
bush or grommet round the shank where it passes through a 
metal part. Rigid pipe-connections can usually be avoided by 
fitting flexible pipe-couplings; with recent advances in the 
manufacture of synthetic elastomers (rubber-like substances) it 
is possible to prevent deterioration of the coupling by most of 
the fluids commonly passed through pipes. 



APPENDIX V 

TABLES 


The tables included in this appendix are of quantities frequently 
used in vibration calculations. Brief notes are given below. 

Table X. Static deflection under gravity and natural frequency. 

For suspensions that can be regarded as essentially equivalent 
to the system of Fig. 1, page 2, the relation between the natural 
frequency and the static deflection under gravity is: 

where <S f = static deflection in the suspension under gravity, in 
inches, 

and F = natural frequency of vertical vibration, in c.p.m. 
The values given in the table are in thousandths of an inch, and 
at intervals of 100 c.p.m. over the range 100-9,900 c.p.m. 

Example. 

An engine-bed weighing 20 tons is to be supported on 10 
flexible mountings so that the natural frequency for vertical 
vibration is 280 c.p.m. Assuming all the mountings to be equally 
loaded and to have the same stiffness, what must this stiffness be ? 

Static deflection for 2,800 c.p.m. (from table) is 0-00449 inch ; 
hence deflection for 280 c.p.m. is 0-449 inch. Load on each mount¬ 
ing is 2 tons; hence stiffness must be 2 x 2,240/0-449 = 9,980 
or about 10,000 lbs. per inch. 

Table XI. Values of (o 2 for various frequencies. 

The values of to* are given in the table, at intervals of 100 
c.p.m. over the range 100-9,900 c.p.m. The quantity given in 
the table must be multiplied by the power of 10 given in the 
column headed (A). 

Examples. 

(i) For F = 5,200 c.p.m., to* = 29-7 x 10 4 (sec. - *). 

(ii) For F = 520 c.p.m., co 2 = - - 7 ^ 8 ■— = 29-7 x 10* (sec. -2 ). 

(iii) For F = 800 c.p.m., co* = 70-2 x 10* (sec. -2 ). 
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Table XII. Values of the dynamic magnifier (no damping). 

The table gives the numerical values of the dynamic magnifier 
(10.2), page 30, m the absence of damping forces. X is the 
ratio of the vibration frequency to the natural frequency. 


Table X 

Static deflection under gravity (thousandths of an inch) for various natural 

frequencies (c.p.m.). 

(Note .—If F is multiplied or divided by 10, deflection is divided or 

multiplied by 100.) 


F 

0 

100 

200 

300 

400 

500 

600 

700 

800 

900 

0 

_ 

3,519 

880 

391 

220 

141 

97*8 

71*8 

55*0 

43*5 

1,000 

35-2 

29*1 

24*4 

20*8 

18*0 

15*6 

13*7 

12 *2 

10*9 

9*75 

2,000 

8*80 

7*98 

7*27 

6-65 

6*11 

5*63 

5*20 

4*83 

4*49 

4*18 

3,000 

3*91 

3*66 

3*44 

3*23 

3*04 

2*87 

2*72 

2*57 

2*44 

2*31 

4,000 

2*20 

2*09 

1*99 

1*90 

1*82 

1*74 

1*66 

1*59 

1*53 

1*47 

5,000 

1*41 

1*35 

1*30 

1*25 

1-21 

1-16 

M2 

1-08 

1*05 

1*01 

6,000 

0*978 

0*946 

0*915 

0*887 

0*859 

0*833 

0*808 

0*784 

0*761 

0*739 

7,000 

0*718 

0*698 

0*679 

0*661 

0*643 

0*626 

0*609 

0*594 

0*579 

0*564 

8,000 

0*550 

0*536 

0*523 

0*511 

0*499 

0*487 

0*476 

0*465 

0*454 

0*444 

9,000 

0*435 

0*425 

0*416 

0*407 

0*398 

0*390 

0*382 

0*374 

0*366 

0*359 


Table XI 

Values of co 2 (radians 2 /sec. 2 , radians/sec. 2 , or sec. -2 ) for various 
frequencies (c.p.m.). 


F 

0 

100 

200 

300 

400 

500 

600 

700 

800 

900 

(A) 

0 

0 

M0 

4*39 

9*87 

17*5 

27*4 

39*5 

53*7 

70*2 

88*8 

X 10 2 

1,000 

M0 

1*33 

1*58 

1*85 

2*15 

2*47 

2*81 

3*17 

3*55 

3*96 

X 10* 

2,000 

4*39 

4*84 

5*31 

5*80 

6*32 

6*85 

7*41 

7*99 

8*60 

9*22 

» 

3,000 

9*87 

10*5 

11-2 

11*9 

12*7 

13*4 

14*2 

15*0 

15*8 

16*7 

» 

4,000 

17*5 

18*4 

19*3 

20*3 

21*2 

22*2 

23*2 

24*2 

25*3 

26*3 

99 

fr,000 

27*4 

28*5 

29*7 

30*8 

32*0 

33*2 

34*4 

35*6 

36*9 

38*2 

99 

6,000 

39*5 

40*8 

42*2 

43*5 

44*9 

46*3 

47*8 

49*2 

50*7 

52*2 j 

99 

7,000 

63*7 

55*3 

56*8 

58*4 

60*0 

61*7 

63*3 

65*0 

66*7 

68*4 

99 

8,000 

70*2 

71*9 

73*7 

75*5 

77*4 

79*2 

8H 

83*0 

84*9 

86*9 

99 

9,000 

88*8 

90*8 

92*8 

94*8 

96*9 

99*0 

101 

103 

105 

107 

99 


(Note .—The quantity given in the table must be multiplied by the factor 
in col umn (A). If F is multiplied or divided by 10, co 2 is multiplied 
or divided by 100.) 


L 
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Table XII 

Dynamic magnifier [numerical value of M = 1/(1 — X*)] for various 

, „ „ vibration frequency, 

values of X , „ 

natural frequency 


X 

M 

X 

M 

X 

M 

01 

101 

11 

4-76 

2-4 

0-210 

0-2 

104 

1-2 

2-27 

2-5 

0-190 

0-3 

110 

1-3 

1-46 

2-6 

0174 

0-4 

119 

1-4 

1-04 

2-8 

0-146 

0-5 

1-33 

1-6 

0-800 

1 3-0 

0-125 

0-6 

1-56 

1-6 

0-641 

3-5 

0-088 

0-7 

1-96 

1-8 

0-446 

4-0 

0-067 

0-8 

2-78 

2-0 

0-333 

4-5 

0-052 

0-9 

1 

5-26 

! 2-2 ; 

: 1 

0-260 

j 6-0 

0-042 






MISCELLANEOUS exercises 

(Numerical examples on the work of Chapters I-IV.) 

All frequencies should be given in c.p.m. 

Chapter I 

Find the natural frequency of a system similar to Fig. 1, in which 
the constants are: 

1. k « 4,000 lbs./in., m = 2 slugs. 

2. k = 3,000,000 poundals/ft., m = 150 lbs. 

3. k = 7,500 lbs./in., m = 200 lbs. 

4. k = 3,750 lbs./in., m » 100 lbs. 

Find the natural frequency of a system similar to Fig. 4, in which the 
constants are: 

5. J x «= V* - 2, C - 3,500,000. 

6. JTj = 3, J a - 4, C = 1,700,000. 

7. « 8, J t = 6, C = 3,400,000. 

(inertias in lbs.ins.secs. 1 , stifEnesses in lbs.ins./radian.) 

8. Determine the period of oscillation of a compound pendulum 
similar to Fig. 3c, in which the constants are: 

m = 10 lbs., I A = 3-73 lbs.ins.secs. 2 , h = 8 inches. 

Chapter II 

Find the resonant frequency, and the natural frequency if the damping 
is neglected, of a system similar to Fig. 9, 'in which the constants are: 

9. k = 200 lbs./in., m = 2 slugs, c = 40% critical. 

10. k = 3,500 lbs./in., m = 1*2 slugs, c = 12 lbs.ins.-i-secs. 

11. k =* 40,000 lbs./in., m = 3 slugs, c = 20% critical. 

12. Find the resonant frequency of a torsional system equivalent to 
Fig. 9, in which the inertia is 10 lbs.ins.secs. 2 , the shaft stiffness is 4,500,000 
lbs.ins./radian, and the damping force is 45% of the critical value. 

13. Find the force transmitted to the base in a system similar to 
that of Section 12, Example I, in which the mass is 5,000 lbs., the springs 
have a total stiffness of 12,500 lbs./in., the damping force is 25% of the 
critical value, and the applied force has a maximum value of 750 lbs. 
and a frequency of 1,500 c.p.m. 

Chapter III 

Find the natural frequencies of the system of Fig. 16, in which the 
constants are: 

14. k x = Jfc, =» 1,000 lbs./in., = m, = 200 lbs. 

15. k x * 300, k t = 450 lbs./in., * 20, m a = 30 lbs. 

147 
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16. Jbj .= 600, k t = 900 lbs./in., m l = 40, m, = 60 lbs. 

Determine the ratio a ,/Oj at the natural frequencies in nos. 14-16. 
Find the natural coupled and uncoupled frequencies of the system 

of Fig. 21; in which the constants are: 

17. m = 1-2 slugs, J = 40 lbs.ins.secs. 1 , Jkj = 2,000, k t = 3,000 lbs./in., 
di =» d, = 6 ins. 

18. m = 1*6 slugs, J = 200 lbs.ins.secs. 1 , Jk t * 3,000, k t = 4,600 
lbs./in., d x = 6, d 2 = 4 ins. 

19. m = 2-25 slugs, J = 300 lbs.ins.secs. 1 , fc x = jfc, = 4,000 lbs./in., 
d x *= 6, d 2 = 4 ins. 

20. m = 2-475 slugs, 7 = 330 lbs.ins.secs.>, k t = jfc, = 8,800 lbs./in., 
d x = 6, dj = 4 ins. 

Chapter IV 

21. Find the lowest natural frequency of a seven-mass torsional 
system in which: 

J 1 =J t =J 3 =J i -J 6 =J t = 0-39, J 1 = 11-1 Ibs.ins.secs. 1 , and 
Cj = G t = C 3 = C 4 = C { = 7-25, C, = 4-55 lbs.ins./micro-radian. 

Determine the swinging form at this frequency. 

22. Find the lowest natural frequency of a geared torsional system 
which is as follows: 

(inertias in lbs.ins.secs. 1 , stiffnesses in lbs.ins./micro-radian.) 

J l = 0-505 C x = 8-788 

J t = 0-521 C t = 8-793 

J 9 = 0-501 C, = 8-793 

J t = 0-501 C« = 8-793 

J t = 0-521 C 5 = 8-788 

J t = 0-556 C t = 2-751 

J 7 = 0-121 

(J 7 geared to J s with a 0-477:1 reduction gear) 

J t = 1-027 C 7 = 11-182 

J, = 3-060 

Determine the swinging form at this frequency. 



ANSWERS TO EXERCISES 

Chapter I 

1. * - C^+Cfl*, where (j) = *(- P±Vp* - 4Q). 

2. (1+Jj/J.Mx = Vfi*+<SVo> , » tan <j> = o>B/S, 

il = d - 5/(1 +Ji/J t ), B = ©! - R/(\+J x /J t ) 9 where i? = - © 2 , 

£ = — d 2 , and a)* = (/(Jj + For no vibration, <0*22* +5* = 0; 

for no steady rotation, J 1 & 1 +J i 0 t = 0 = J x 6 t +J t d r 

3. y = C^+Cje-^+Cae^+^e-^. 


Chapter II 


1. For type (6), x 0 /a? 0 < — y, where z 0 and x 0 are initial displacement 
and velocity. 

Chapter III 

2. a, = j k ------ 2 a x , % being given by (15.3). 

3* fcj 53 0. 

Chapter IV 

1. First six natural frequencies are approximately 1,900, 4,100,6,100, 
8,620, 9,200, and 9,550 c.p.m. ; seventh is outside range of Table TTT 
(Note : Z iA = Z 2B = 0 is a solution of (22.6) giving three of the frequencies.) 


d*y 


Chapter V 


nx 


I. (a) y = a, g^- 2 = 0 at each end; (b) y 0 = a — Y sin —. a is found 


l 


to be 2Y/n by determining condition for zero total inertia force (my) for the 
whole beam ; (c) (K.E.)^. = 0 048 o>'Y*M B , (P.E.W = EIn*Y*/ti *; 
(< d) cd 1 = 519 EI/M B l\ 

Chapter VI 


1 . 


a 0 = a/hn 9 a k = 0, b k = a/2 k ; and for m^k 9 


a m - 


a * (1 + cos rrm/k); b m = 


n (m % — &*)^ ”* n (m*-*rk % ) 

2. a 0 =* ^/n ; a n = (2/nn) sin ; 6 n = 0. 


sin m rc/t 


Miscellaneous Set 
(Frequencies in c.p.m.) 

1. 427. 2. 1,350. 3. 1,150. 4. 1,150. 

6. 9,500. 7. 9,500. 8. 1*36 secs. 

9. Resonant, 78*6 natural, 95*5 
10. „ 511 „ 516 


5. 17,900. 
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11. Resonant, 1,070 natural, 1,110. 

12. 4,940. 13. 560 lbs. 

14. 269, 679; ratios: 1-615, - 0-618. 

15. 407, 1,297 ; „ 1-98, - 0-498. 

16. As for 15. 

17. Coupled, 503, 642; uncoupled, 534, 615. 

18. „ 286, 675; „ same as coupled. 

19. „ 206, 585; „ 251, 570. 

20. „ 292, 828; „ 355, 807. 

21.9,940; 

A, 4, A 3 4, 4, 4, A, 

1 0-942 0-829 0-667 0-467 0-239 - 0145 

22. 11,990; 

A, A 2 A 3 Ag A 3 A 3 A 7 A 3 A 3 

1 0-909 0-734 0-492 0-207 - 0-099 - 1-042 - 0-497 - 0-875 
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Absorber, undamped vibration: 

linear, 48 ; torsional, 49 
Alternation of zeros and discon¬ 
tinuities, 65 

Amplitude, 7 ; convention regard¬ 
ing, 10 

Analogues, linear-torsional, 12 
Analysis, Fourier, see Fourier 
Anti-vibration mountings, details, 
141 ; theory, 134 

Approximations: natural frequen¬ 
cies of beams (Rayleigh), 92 ; 
numerical solutions of equa¬ 
tions, 131 


Complex number method, 38, 77 
Complex vibrations, 99 
Constant excitation, 30 
Constant of integration, 3 
Continuous systems, see Beams and 
Shafts 

Convention regarding amplitude 
and phase, 10 

Coupled systems : flexibly-coupled, 
76 ; rigidly-coupled, 68 
Coupled vibrations : two degrees of 
freedom, 50 ; three degrees of 
freedom, 53, 129 ; six degrees 
of freedom, 54, 130 


Beams, theory of bending of, 12 
Beams, heavy, vibrations of, 86; 
Rayleigh’s method, 92 ; special 
cases—cantilever, 89, 96 ; 

clamped-clamped, 96; free- 
free, 90 ; simply supported, 94, 
96 

Beams, light cantilever : vibration 
with mass on free end, 11 ; 
stiffness formula, 123 
Beating waveform, 106 
Bending of beams, theory, 121 
Bending vibrations in beams, see 
Beams 
Biot, 65, 73 

Cantilever, heavy, vibrations of, 89 ; 

deformation curves, 91 
Cantilever, light: vibration with 
mass on free end, 11; stiffness 
formula, 123 
Circular functions, 112 
Circular shafts, see Shafts 
Clamped-clamped heavy beam, 96 
Complementary function, 27 


Damped systems, 18; dynamic 
stiffness of, 76 
“ Damper,” 49 

Damping, 18; coefficient, 20; 
critical, 23, 41 ex. 4; effect 
on insulation, 36; electro¬ 
magnetic, 39; experimental 
determination of coefficient,25 
Dashpot, 20 

Decay of motion, 2, 19, 24 
Decrement, logarithmic, 25 
Deformation diagrams : heavy can¬ 
tilever, 91; light cantilever, 
123 ; spring-mass systems, 47. 
See also Swinging forms 
Den Hartog, 10, 4# 

Determinants, 126 
Differential operator D, 4, 88 
Dissipation of energy, 2, 19 
Dynamic absorber: linear, 48; 
torsional, 49 

Dynamic magnifier, 29; table of 
values, 146; value at reson¬ 
ance, 30; relation between 
d.m. and phase angle, 34 
Dynamic modulus, 73 
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Dynamic stiffness, 72; damped 
systems, 76; modulus of, 77 ; 
transference formulae, 74 

Effective inertia, 00 ; curves, 64, 
70; experimental determina¬ 
tion, 69 

End conditions : heavy beam, 89 ; 
heavy shaft, 84 

Energy : dissipation, 2,19 ; kinetic, 
of beam, 92; potential, of 
beam, 93 

Engines: dynamic absorbers in, 
49 ; gas-forces in, 107 ; mount¬ 
ings for, 54 

Equations: frequency, 14, 129; 
normalising of, 133 ; numerical 
solution of, 131 

Equations of motion : compound 
pendulum, 13 ; coupled vibra¬ 
tions, 51,129,130 ; heavy beam, 
87 ; heavy shaft, 82 ; multi¬ 
mass systems, 59 ; one degree of 
freedom, 3, 21, 26; seismic 
vibrograph, 37 ; simple pendu¬ 
lum, 12 ; torsional pendulum, 
13; torsional two-mass system, 
14; two degrees of freedom, 44 ; 
undamped linear absorber, 49 

Excitation : “ constant,” 30, 36 ; 

gas-force, 107 ; “ inertia,” 31, 
36 

Experimental determinations : 
damping coefficient, 25 ; effec¬ 
tive inertia, 69 

Exponential function, 111 

Flexibly coupled systems, 76 

Flexural vibrations of beams, see 
Beams 

Forced motion, 25 

Fourier: analysis, 102, 106 ; series, 
100; synthesis, 101 

Free-free heavy beam, 90 

Frequency, 8 ; equations, 14, 129 ; 
natural, 8; natural damped, 
25 ; resonant, 30, 41 ex. 3; 
tables, 71. See also Natural 
frequency 


Function : complementary, 27 ; 
exponential, 111 ; hyperbolic, 
113 ; sine and cosine, 112 
Fundamental frequency, 100 

Gas-force harmonics, 107 
Gearing in torsional systems, 75 
Gibbs’ phenomenon, 105 
Gravity deflection and natural fre¬ 
quency : formula, 9; table, 
145 

Harmonic, 100; analysis, see 
Fourier; gas-forces, 107; 
motion, 11 

Heavy beams and shafts, see Beams 
and Shafts 

Hype$>olic functions, 113 

Impedance, mechanical, 73. See 
also Dynamic stiffness 
Inertia: effective, see Effective 
inertia; excitation, 31, 36; 
moment of, 12, 118 ; partial, 61 
Initial conditions in single degree of 
freedom system: undamped, 
3, 6, 17 ex. 2 ; damped, 40 exs. 
1 and 6 

Insulation, 35, 134; effect of 

damping on, 36 
Integral, particular, 27 
Integration, constants of, 3 

Kinetic energy of beam, 92 

Light(-weight): beams, see Beams ; 

shafts, see Torsional systems 
Linear spring characteristics, 10 
Linear-torsional analogues, 12 
Logarithmic decrement, 25, 

Magnifier, dynamic, see Dynamic 
magnifier 

Manley, 32, 57, 58, 105 
Measurement systems and units, 9, 
12, 73, 117 

Mechanical impedance, 73. See 
also Dynamic stiffness 
Mode, 48 
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Modulus: dynamic, 73; of 

dynamic stiffness, 77, 108 ; of 
shear, 124 

Moment: of area, second, 118 ; of 
inertia, 12, 118 
Morris, viii 

Motion : damped, 20 ; decay of, 2, 
10, 24; equations of, see 

Equations of motion ; forced, 
25 ; free, 2 ; harmonic, 11 ; 
Newton’s second law of, 117 ; 
steady-state, 28 
Multi-mass systems, 57 

Natural frequency, 8 ; and gravity 
deflection, 9, 145 

Natural frequency formulae : one 
degree of freedom, 11 ; canti¬ 
lever, 89, 90; clamped-clamped 
beam, 96; compound pendu¬ 
lum, 13 ; free-free beam, 90, 
90 ; shafts, 80 ; simple pendu¬ 
lum, 13; simply-supported 
beam, 95 ; torsional pendulum, 
13; torsional two-mass system, 
15 ; two degrees of freedom, 
46 ; damped single-mass sys¬ 
tem, 25 

Newton’s second law of motion, 117 
Nodes, 47 ; and effective inertia, 68 
Normalising of numerical equations, 
133 

Numerical: Fourier analysis, 106 ; 
solution of equations, 131 

Operator; differential, 4, 88; 

versor, 38, 77, 114 

Partial inertias, 61 
Particular integral, 27 
Pendulum : compound, 13 ; simple, 
12; torsional, 13 
Period of vibration, 8 
Periodic functions, 100 
Phase angle, 7; and dynamic 
magnifier, 34 ; convention re¬ 
garding, 10; value at reson¬ 
ance, 33 


Phase velocity, 13 ; and frequency, 
14 

Potential energy of beam, 93 

Rayleigh’s energy method, 92 
Reciprocation theorem, 52 
Resistance, 20 

Resonance, 28 ; and dynamic stiff¬ 
ness, 76; and natural fre¬ 
quency, 32, 63 ; elimination of, 
29 

Response of seismic vibrograph, 39 
Rigidly coupled systems, 68 
Rolling oscillation of half-cylinder, 
119 

Runge’s method of Fourier analysis, 
106 

Scalars, 113 

Seismic vibrograph, 37 ; response 
curve, 39 

Series: exponential. 111; Fourier, 
100 ; sine and cosine, 112 
Shaft, heavy, 81 ; dynamic stiff¬ 
ness, 84 ; natural frequencies, 
86 

Shaft, light, see Torsional systems 
Shear modulus, 124 
Sine : wave, 7 ; and cosine series, 
112 

Skew-symmetry in periodic func¬ 
tions, 102 
Slugs, 9, 117 
Spectrum, frequency, 53 
Spinner, aero-propeller, 54 
Spring stiffness: linear, 3; tor¬ 
sional, 13 

Starting transient, 32, 41 ex. 6 
Static gravity deflection and natural 
frequency : formula, 9 ; table, 
145 

Steady-state motion, 28 
Stiffness : dynamic, see Dynamic 
stiffness ; linear, 3 ; torsional, 
13 

Swinging form, 05, 67 
Symmetry : in mounting systems, 
54, 130 ; in periodic functions, 
102 ; in spring characteristic, 10 
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Timoshenko, 90, 92 
Torque summation tables, 71 
Torsional-linear analogues, 12 
Torsional systems, heavy shafts, see 
Shafts 

Torsional systems, light shafts: 
one degree of freedom, 13 ; two 
degrees, 48 ; three degrees, 63 ; 
many degrees, 57 
Torsion of circular shaft, 124 
Transference formulae, dynamic 
stiffness, 74; heavy shaft, 
84 

Transient, starting, 32, 41 ex. 6 
Transmitted force, 32, 138 


“ Uncoupled ” frequencies, 53 
Uniform beams and shafts, see 
Beams and Shafts 
Units of measurement, 9,12, 73,117 

Vector, 113; rotating, 7; v. 

method, 38, 77 
Versor operator, 38, 77, 114 
Vibrograph, seismic, 37; response 
curve, 39 

Whittaker, 65 ; and Watson, 112 
Wilson, Ker, 50, 78, 107 
Work done by impressed forces, 19, 
108 






